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Executive  Summary 


This  is  the  final  report  for  the  Source  Characterization  Model  (SCM)  developed  by  Aerodyne  Research, 

Inc.  (ARI)  in  collaboration  with  the  U.S.  Army  Environmental  Center  (USAEC).  It  describes  work 
performed  under  contract  DACA72-00-C-0001  from  the  Strategic  Environmental  Research  and 
Development  Program  (SERDP)  for  the  project  CP-1 159  entitled  "A  Predictive  Capability  for  the  Source 
Terms  of  Residual  Energetic  Materials  from  Burning  and/or  Detonation  Activities".  A  portion  of  this 
report  was  previously  released  as  a  supplement  to  the  2002  Annual  Report  for  that  SERDP  project. 

Detonation  of  energetic  materials  produces  a  wide  range  of  air  and  surface  pollutants,  including  carbon 
monoxide,  nitrogen  oxides,  volatile  organic  compounds,  acid  gases,  and  particulate  matter.  These 
emissions,  including  undecomposed  or  partially  decomposed  energetic  materials,  may  lead  to  atmospheric 
pollution  or  ground  water  contamination.  The  speciation  and  amounts  of  these  emitted  pollutants  depend  on 
the  identity  and  amount  of  energetic  material  detonated,  the  detonation  order,  the  detonation  mode  (air 
burst,  surface  detonation,  buried  detonation),  and  the  munitions  type  (shell,  mine,  detonation  charge,  etc.). 

In  order  to  determine  whether  U.S.  military  training  or  munitions  disposal  activities  produce  emissions  that 
threaten  air  or  ground  water  quality,  reliable  estimates  of  emission  factors  must  be  available  for  a 
representative  fraction  of  the  thousands  of  munitions  types  in  the  inventory.  A  variety  of  recent  detonation 
experiments,  including  detonation  chamber  tests,  and  open-air  field  detonations  have  yielded  increasingly 
rich  data  sets  of  measured  gaseous  and  particulate  emissions  factors  from  a  range  of  energetic  materials 
deployed  as  both  unconfined  charges  and  standard  munitions.  These  data  have  been  used  under  CP-1 159  to 
guide  the  development  of  a  detonation  source  characterization  model  (SCM)  to  be  used  to  predict  emissions 
across  munitions  classes  and  detonation  modes.  The  SCM  predicts  the  chemical  and  physical  distribution 
of  detonation  products  in  the  initial  detonation  fireball,  calculates  the  amount,  chemical  influence,  and 
deposition  distribution  of  dirt  entrained  into  the  fireball,  and  the  effects  of  afterburning  as  oxygen  rich  air  is 
entrained  into  the  hot  rising  plume  on  both  gas  phase  and  particulate  emissions,  and  estimates  the  amount 
of  shock  ablated  energetic  material  that  escaped  the  fireball.  SCM  predictions  have  been  validated  against 
measured  emissions  factors  for  a  variety  of  munitions.  The  validated  SCM  can  be  used  to  extend  emission 
factor  predictions  to  unmeasured  munitions  types  and  to  compute  input  data  for  traditional  atmospheric 
dispersion/reaction  models  and  soil/water  transport/ transformation  models.  The  SCM  was  developed  at 
Aerodyne  Research,  Inc.  utilizing  munitions  composition  and  measured  detonation  products  databases 
developed  by  collaborators  at  the  U.S.  Army  Environmental  Center. 

SCM  Model  Summary 

As  documented  in  this  report,  the  SCM  includes  the  following  sub-models  for  the  critical  processes 
governing  the  post-detonation  plume  and  chemical  evolution  of  detonation  emission  products. 

Detonation  Thermochemistry:  predicts  thermodynamic  properties  and  chemical  speciation  of  the  high 
temperature  and  pressure  Chapman- Jouquet  (C-J)  detonation  state  and  the  adiabat  for  the  rapid  isentropic 
expansion  of  detonation  gas  towards  atmospheric  pressures. 

Crater  Dynamics:  predicts  crater  size  and  the  mass  and  particle  size  distribution  for  ejected  soil  partitioning 
into  the  buoyant  cloud. 

Buoyant  Puff  Rise:  predicts  the  chemical  speciation,  thermodynamic  properties,  height,  and  volume  of  the 
buoyant  detonation  puff  when  it  has  equilibrated  with  the  ambient  atmosphere. 

Buoyant  Puff  Structure:  treats  the  vortex  structure  of  the  detonation  puff  including  the  effect  on  particle 
deposition. 

Particle  Deposition:  treats  particle  deposition  due  to  gravitational  settling  and  the  impact  of  particle  fall-off 
on  buoyant  puff  rise. 

Meteorology:  enables  the  treatment  of  buoyant  puff  rise  in  a  one-dimensional  stratified  atmosphere. 
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Chemical  Kinetics:  treats  the  afterburning  chemistry  in  the  post-detonation  plume  with  air  entrainment. 


Collectively,  the  sub-models  listed  above  enable  the  SCM  to  support  environmental  impact  assessments  at 
DoD  test  and  training  sites  in  two  specific  ways: 

•  SCM  predictions  for  detonation  plume  properties  and  chemical  speciation  are  important  to  near  term 
(to  the  range  fence)  air  quality  issues  and  provide  the  necessary  inputs  to  longer  time  atmospheric 
dispersion  models  to  assess  air  quality  issues  over  distances  beyond  the  specific  DoD  site. 

•  SCM  predictions  for  particle  depositions  are  important  to  near  term  ground/ground  water  pollution 
issues  and  provide  necessary  inputs  to  fate  and  transport  models  for  environmental  assessments. 

Project  Highlights 

SERDP  has  sponsored  the  development  of  a  computer  model  that  predicts  potential  air  and  ground  water 
pollutants  produced  by  exploding  munitions  and  demolition  charges.  Developed  at  Aerodyne  Research, 
Inc.,  the  computer  model,  termed  a  detonation  source  characterization  model  (SCM),  draws  on  new 
databases  compiling  munitions  compositions  and  detonation  products  developed  by  collaborating  scientists 
at  the  Army  Environmental  Center. 

The  detonation  source  characterization  model  (SCM)  is  composed  of  a  number  of  coordinated  sub-models. 
These  sub-models  predict  the  chemical  content  of  the  detonation  fireball,  compute  the  mass  of  fireball 
entrained  soil  particles,  estimate  the  amount  of  shock  ablated  energetic  material  that  escapes  beyond  the 
fireball,  model  the  entrainment  of  atmospheric  oxygen  mixed  into  the  fireball  and  the  subsequent 
afterburning  effects  on  the  chemical  content  of  the  detonation  plume,  predict  the  chemical  interactions 
among  both  gas  phase  pollutants  and  soot  and  dirt  particles  in  the  plume,  calculate  the  deposition  pattern  of 
larger  plume  entrained  particles  back  to  the  surface,  and  predict  the  buoyant  rise  of  the  plume  gases  and 
entrained  fine  particles. 

The  detonation  SCM  has  been  validated  with  pollutant  emission  factor  data  obtained  for  a  range  of 
munitions  and  demolition  charges  obtained  in  an  on-going  series  of  detonation  tests  managed  by  the  Army 
Environmental  Center.  After  validation  against  a  sufficient  range  of  diverse  munitions  and  explosives,  the 
SCM  can  be  used  to  predict  emissions  factors  for  as  yet  unmeasured  types  of  munitions. 

SCM  predicted  emissions  factors  can  be  used  as  inputs  for  traditional  atmospheric  dispersion  or  soil  and 
ground  water  transport  models  that  may  be  used  by  environmental  scientists  and  engineers  to  predict  the 
impact  of  military  training,  weapons  testing,  or  munitions  disposal  activities  on  military  facilities  and 
surrounding  regions. 

Under  an  additional  task,  ARI  deployed  an  innovative  aerosol  mass  spectrometer  (AMS)  instrument  to 
support  the  SERDP  project  CP- 1 197,  consisting  of  field  measurements  by  Battelle,  Columbus,  OH  and 
managed  by  USAEC,  by  measuring  the  non-refractory  chemical  composition  of  fine  particles  from  firing 
point  muzzle  exhaust  plumes.  Muzzle  exhaust  plumes  contained  in  the  Aberdeen  Test  Center’s  (ATC) 
point-of-discharge  test  chamber  were  sampled  for  about  20  minutes  after  the  firing  of  the  each  munition. 
The  non-refractory  chemical  composition  of  the  evolving  muzzle  exhaust  plume  fine  particles  size 
distribution  was  characterized  during  these  measurements.  A  summary  of  these  results  and  a  citation  for 
the  task  report  we  presented  is  given  in  Appendix  A. 
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1.0  Introduction 


A  Source  Characterization  Model  (SCM)  has  been  developed  under  SERDP  funding  (Project  CP-1159)  to 
provide  a  predictive  numerical  tool  for  gaseous  and  particulate  emissions  that  result  from  HE  impact  and 
'blow-in-place'  detonations  at  DoD  test  and  training  sites.  Following  a  brief  review  of  the  background  for 
this  work  and  a  brief  summary  of  the  SCM  features  and  capabilities,  this  report  presents  (1) 
phenomenology  and  documentation  for  processes  treated  by  the  SCM,  (2)  a  description  of  the  USAEC 
Army  Range  and  Munitions  Emissions  Environmental  Data  (ARMED),  and  (3)  results  of  model  testing  and 
validation  studies.. 

1.1  Background 

1.1.1  Environmental  and  Public  Safety  Impact  of  Munitions  Emissions 

U.S.  Army  operations  at  test  and  training  sites  involve  the  detonation  of  High  Energy  (HE)  ordnances.  The 
chemical  emissions  from  exploding  ordnances  pose  a  potentially  dangerous  environmental  and  health 
hazard.  The  identification  of  high  levels  of  residual  energetic  material  at  the  Massachusetts  Military 
Reservation  (MMR)  significantly  heightened  public  concern  about  the  long-term  consequences  of  test  and 
training  operations  at  U.S.  Army  sites  in  general.  As  a  consequence,  the  U.S.  Army  has  identified  the  need 
to:  (1)  quantify  the  ordnance  detonation  emissions  and  (2)  assess  the  potential  health  and  environmental 
impact  of  these  emissions  as  critical  to  responding  to  regulatory  requirements  and  public  concerns,  as  well 
as  in  support  of  defining  and  possibly  modifying  test  and  training  practices. 

1.1.2  Emissions  Characterization  Test  Programs 

Several  major  emissions  testing  programs  have  been  undertaken  over  the  past  decade.  These  emission 
measurement  programs  were  largely  driven  by  concern  over  the  emissions  during  open  burn  and  open 
detonation  (OB/OD)  operations  and  were  designed  to  assist  in  meeting  regulatory  compliance  requirements 
of  the  CAAA,  EPCRA  and  RCRA.  Examples  are  the  testing  programs  undertaken  at  the  Nevada  Test  Site 
(NTS)  [Johnson,  1992]  and  Dugway  Proving  Ground  (DPG)  [Black,  1998].  The  Environmental  Protection 
Agency  (EPA)  published  data  from  DPG  testing  [Mitchell  and  Suggs,  1998]  that  documented  the  identities 
and  amounts  of  OB/OD  products  from  the  common  bulk  explosives.  The  resulting  database  greatly 
expanded  the  scarce  test  data  available  and  provided  a  proven  protocol  for  further  testing  and  modeling. 

More  recently,  emission  factors  for  a  large  number  of  gaseous  and  particle  emissions  have  been  measured 
by  researchers  at  the  Aberdeen  Test  Center  (ATC)  through  the  Exploding  Ordnance  Emission  Study,  Phase 
II  (EOES)  funded  by  the  U.S.  Army  Environmental  Center.  USAEC  has  also  sponsored  emission 
characterization  tests  at  Dugway  Proving  Ground  (DPG)  [Carson,  et  al.,  2003]  and  a  Firing  Point  Emission 
Study  at  the  ATC. 

1.1.3  Energetic  Residue  from  Munition  Detonation 

Analyses  of  the  residues  following  munition  detonations  have  confirmed  explosive  residues  that  can  result 
in  contamination  of  the  soil  and  underlying  ground  water  [Yinon  and  Zitrin,  1993;  Jenkins  et  al.,  2001; 
Thiboutot  et  al.,  1998,  Pennington  et  al.,  2001,  2002,  USACHPPM,  2000,  and  Hewitt  et  al.,  2003]. 

1.1.4  OB/OD  and  Detonation  Emissions  Modeling 

A  number  of  models  have  been  developed  for  predicting  emissions  from  OB/OD.  Examples  include 
ADORA,  PCAD,  POLU,  OBODM  and  Merlin.  These  models'  predictive  ability  is  limited  by  the 
uncertainties  in  the  measured  test  data  that  have  been  used  to  tune  the  models  and  by  the  assumption  that  a 
mass  balance  focused  on  the  major  emissions  would  adequately  represent  detonation.  For  four 
representative  scenarios,  the  predictions  of  some  of  these  models  compared  poorly  with  available 
DPG/NTS  test  data  [Eck  et  al.,  1998].  The  available  databases  and  models  still  need  to  be  improved  to 
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meet  the  public  and  regulatory  requirements  [Banta  et  al.,  1998].  None  of  the  current  models  predict  the 
particle  size,  density,  and  composition  of  condensed  particles  in  the  plume.  The  dry  deposition  patterns 
into  the  various  surrounding  environments  are  not  differentiated  in  a  single  run.  The  process  of  wet 
deposition,  for  example,  is  not  included  in  any  of  the  models.  Existing  models  fail  to  address  the  following 
problems: 

(1)  Quantification  of  afterburning  (after  detonation)  effects; 

(2)  Emission  factors  for  VOCs  or  semi-VOCs  which  are  typically  emitted  in  the  parts  per  million  range; 

(3)  Predictive  quantification  of  mass  loadings  and  size  distributions  for  emitted  particles  (smoke); 

(4)  Prediction  of  dry  and  wet  deposition  patterns  for  emitted  particles  and  acid  gases  on  various  ground 
surfaces  or  in  waters; 

(5)  Potential  interactions  of  emitted  pollutants  with  entrained  soil. 

Because  of  these  problems,  regulators  are  reluctant  to  accept  model  outputs  as  representative.  Testing  and 
sampling  continues  to  be  an  expensive  and  difficult  process,  and  it  is  unlikely  that  the  DOD  will  ever 
sample  emissions  from  all  of  the  munitions  in  its  inventory.  Regulatory  interest  in  the  impacts  from 
munitions  detonation  is  high  (for  example,  see  Hunt  and  Huntington,  1998)  and  the  DOD  needs  to  cost- 
effectively  report  emissions  and  impacts  under  the  Emergency  Planning  and  Community  Right-to-Know 
Act,  the  Safe  Drinking  Water  Act,  CERCLA,  and  RCRA.  Thus,  a  predictive  public  domain  model  that  has 
the  capability  of  addressing  all  of  these  major  concerns  is  necessary. 


1.2  Source  Characterization  Model 

1.2.1  Features  and  Capabilities 

The  SCM  addresses  a  number  of  the  technical  issues  that  are  important  in  estimating  munition  emission 

levels  and  particulate  deposition  arising  from  DoD  test  and  training  operation.  A  few  of  the  most  important 

features  of  the  SCM  are  the  following: 

•  A  source  module  for  characterizing  the  energy  release  and  chemical  speciation  with  munition 
detonation. 

•  Crater  formation  and  soil  particle  entrainment  at  ground  level  due  to  detonation. 

•  Buoyant  puff  rise  and  the  flow  structure. 

•  Explosive  residue  or  soil  particles  depletion  and  deposition  from  the  puff.  In  this  module,  the  spatial 
and  temporal  distribution  of  the  amount  and  size  of  soil  particles  is  obtained,  which  can  be  used  to 
calculate  the  soil  depletion  terms. 

•  The  effect  of  particle  depletion  and  deposition  (the  fall-out  of  particulate  matter)  on  the  puff  dynamics 
has  been  formulated.  Modifications  to  the  laws  of  mass  conservation,  momentum  conservation,  and 
energy  conservation  are  obtained  to  quantify  mass  transfer  and  heat  transfer  for  the  puff  due  to  particle 
fall-out. 

•  A  chemical  kinetics  module  for  finite-rate  afterburning  chemistry  and  a  particle  microphysics  module 
for  characterizing  the  dynamics  of  particulate  emissions  and  gas-to-particle  conversion. 

•  A  data-driven  model  for  energetic  material  residue  emissions. 

Each  of  these  features  is  discussed  in  more  detail  in  the  subsequent  sections  of  this  report. 
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1.3  Report  Outline 


The  remainder  of  this  report  is  organized  as  follows.  Section  2  presents  a  broad  overview  of  the  SCM  from 
a  phenomenology  perspective  including  a  discussion  of  the  specific  physiochemical  processes  that  are 
treated,  the  computational  processing  flow,  and  model  inputs.  The  SCM  is  described  in  detail  in  Section  3. 
This  includes  a  discussion  of  the  modules  for:  (1)  detonation  thermochemistry,  (2)  crater  dynamics,  (3) 
buoyant  puff  rise,  (4)  detonation  plume  vortex  structure,  (5)  particle  deposition,  (6)  ambient  meteorology, 
and  (7)  afterburning  chemical  kinetics,  (8)  heterogeneous  update  kinetics,  and  (9)  energetic  material 
emissions.  Section  4  describes  the  USAEC  ARMED  database  for  munition  emissions.  Section  5  discusses 
the  initial  data  analysis  and  model  validation  studies  that  have  been  performed.  A  summary  is  presented  in 
Section  6.  Technical  references  are  presented  in  Section  7. 


3 


2.0  SCM  Model  Phenomenology 


2.1  Model  Overview 

Figure  1  presents  a  pictorial  representation  of  the  phenomenology  that  is  the  basis  for  the  SCM  and  the  key 
processes  governing  the  release,  dispersion,  and  deposition  of  emissions  from  the  detonation  of  HE 
explosives.  As  shown  in  Figure  1,  three  distinct  stages  can  be  identified  that  are  differentiated  by  the 
physiochemical  processes  governing  the  important  emissions,  as  well  as  the  characteristic  time  scales. 


KEY  PROCESSES 


Detonation 


Buoyant  Puff  Rise 


Soil  Interactions 


Detonation 
Energetics  & 
Chemistry 


atmospheric  stratification 

Air  Entrainment  & 
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Passive  Dispersion 
Equilibration 


Figure  1.  Schematic  of  the  Source  Characterization  (SCM)  model  showing  the  key  microphysical 
processes  governing  the  release,  dispersion,  and  deposition  of  gaseous  and  particulate  emissions  from  HE 
detonations. 

2.1.1  Stage  1:  Detonation  Energetics  and  Chemistry 

Stage  1  is  taken  to  encompass  the  initial  detonation,  a  very  rapid  isentropic  expansion  of  the  detonation 
plume  to  near  ambient  conditions,  formation  of  a  detonation  crater,  and  the  initial  entrainment  of  any  soil 
particles.  In  addition,  energetic  material  is  ejected  due  to  ablation  and/or  ablation  with  some  pyrolysis. 

The  physiochemical  processes  associated  with  Stage  1  are  highly  complex  and  collectively  involve 
characteristic  time  scales  ranging  from  a  few  picoseconds  to  a  few  microseconds.  Consequently,  a  rigorous 
treatment  of  Stage  1  would  be  exceedingly  difficult  to  implement,  time  consuming  to  execute,  and  subject 
to  a  great  deal  of  uncertainty  due  to  the  lack  of  experimental  data  and  limited  knowledge  of  critical  physical 
and  chemical  parameters.  Consequently,  the  SCM  treats  Stage  1  using  non-ideal  gas  chemical  equilibrium 
in  combination  with  empirical  data  for  NOx  emission  and  post-blast  energetic  residues  to  provide  a  data- 
driven  treatment  of  the  Stage  1  processes. 
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Figure  2.  Important  physical  and  chemical  processes  of  pollutants  emission  and  the  soil  scavenging  effects. 


Despite  their  complexity,  Stage  1  processes  play  a  fundamental  role  in  determining  the  HE  detonation 
emissions.  Figure  2  shows  a  diagram  on  the  important  processes  in  obtaining  the  pollutant  source  terms 
coming  out  of  a  detonation  crater.  Most  of  the  explosive  energy  is  converted  to  thermal  energy  that  drives 
the  fireball  as  it  moves  upwards.  A  considerable  amount  of  the  energy  is  used  to  compress  the  soil  and  lift 
the  soil  particles  to  form  a  crater.  Some  energy  is  dissipated  through  sound  waves  propagated  away  in  the 
air  as  well  as  in  the  ground.  Immediately  after  the  detonation,  some  small  fraction  of  detonation  products 
and  thermal  energy  in  the  fireball  is  transferred  to  the  ground.  When  the  soil  particles  rise  and  fall  through 
the  hot  fireball,  they  scavenge  the  pollutants  as  well  as  the  thermal  energy.  The  eventual  emitted  pollutants 
will  experience  further  chemical  reaction  as  more  humid  air  is  entrained  into  the  hot  fireball  (Stage  2). 


2.1.2  Stage  2:  Buoyant  Puff  Rise 


Stage  2  is  taken  to  include  rise  of  the  detonation  plume  in  a  stratified  atmosphere  due  to  buoyant  forces  as 
well  as  downstream  dispersion  due  to  advection.  Initial  conditions  for  Stage  2  processes  are  based  on  the 
energy  release,  chemical  transformations,  and  soil  interactions  that  comprise  Stage  1 .  The  processes 
governing  the  buoyant  puffs  thermodynamic  properties  and  chemical  composition  include  the  finite  rate 
chemical  kinetics  associated  with  afterburning  processes  in  the  hot  plume  as  oxygen  is  entrained  with 
ambient  air,  as  well  as  the  loss  of  heavier  particles,  including  any  gas  phase  species  adsorbed  on  those 
particles,  due  to  deposition  processes.  Stage  2  is  taken  to  continue  until  the  detonation  puff  has  equilibrated 
with  the  ambient  atmosphere  where  upon  the  subsequent  dispersion  is  governed  by  passive  atmospheric 
dispersion  (Stage  3). 

2.1.3  Stage  3:  Atmospheric  Dispersion 


Stage  3  refers  to  the  atmospheric  dispersion  of  released  gas  phase  and  smaller  particulate  emissions  over 
longer  time  scales  and  spatial  dimensions.  Stage  3  is  not  explicitly  treated  by  the  SCM.  Rather,  the  SCM 
provides  initial  conditions  for  existing  state-of-the-art  atmospheric  codes  that  currently  exist.  The  results  of 
the  SCM  include:  (1)  the  altitude,  downstream  distance,  and  volume  of  the  detonation  plume,  (2) 
thermodynamic  properties  of  the  detonation  plume,  (3)  the  chemical  composition  of  the  detonation  plume 
in  the  form  of  source  emission  factors,  and  (4)  the  near  term  particle  fall-out  pattern  due  to  deposition. 
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2.2  Computational  Processing  Flow 


Figure  3  outlines  the  basic  system  architecture  for  the  SCM  in  terms  of  the  phenomenological  description 
that  is  shown  in  Figure  1.  This  provides  a  more  detailed  description  of  the  processing  flow  adopted  in  the 
SCM,  as  well  as  a  broad  overview  of  the  component  sub-models  that  collectively  are  the  basis  for  the  SCM. 


Source  Characterization  Model 
System  Architecture 
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Dispersion  Models 

Photochemistry 
Dry  an d_W_et _Dep_o_sjt]on 
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Figure  3.  Schematic  of  the  Source  Characterization  Model  (SCM)  system  architecture  showing  the  general 
processing  flow  and  relationship  between  the  key  processes  governing  FIE  detonation  emissions. 


As  indicated  in  Figure  3,  the  basic  inputs  to  the  SCM  fall  into  three  broad  groups  corresponding  to  the  key 
parameters  and  properties  for  the  projectile  or  explosive  charge,  ambient  air,  and  soil.  Figure  4  shows  a 
schematic  diagram  of  model  parameters.  They  include  the  following: 

Group  1:  projectile/charge  parameters:  shell  mass,  shell  volume,  and  explosive  type  and  charge. 

Group  2:  ambient  parameters:  wind  speed,  temperature,  stratification,  ambient  turbulence,  and  humidity. 

Group  3:  Soil  parameters:  soil  density,  compressibility,  soil  strength,  temperature,  particle  size,  and  water 
content. 
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Figure  4.  Model  parameters  for  the  quantification  of  chemical  pollutants  emitted  during  an  explosive 
detonation  in  military  test  and  training  ranges. 

As  an  example,  for  a  bare  charge,  such  as  a  simple  TNT  demolition  charge,  the  input  would  include  the 
chemical  composition  and  total  charge  mass.  For  a  projectile,  the  inputs  would  also  include  the 
composition  of  the  shell  casing.  In  addition,  it  is  necessary  to  know  the  chemical  composition  and  mass  of 
any  secondary,  supplemental,  or  detonating  charges  for  blow-in-place  detonations  as  part  of  field 
measurement  programs  or  the  chemical  contribution  from  the  projectile  propellant  for  impact  detonations. 
These  inputs  provide  the  necessary  data  for  an  analysis  of  the  explosive  gas  phase  chemistry  that  forms  the 
first  part  of  Stage  1  in  Figure  3.  That  analysis  provides  estimates  for  both  the  energy  release  and  detonation 
species  emission  factors. 

The  analysis  of  crater  formation  that  forms  the  second  part  of  the  Stage  1  modeling  requires  the  energy 
release  in  the  detonation  coupled  with  additional  inputs  in  order  to  specify  soil  properties  and  detonation 
depth.  The  result  of  that  analysis  is  an  estimate  for  crater  properties  and,  more  importantly,  the  mass  of  soil 
entrained  in  the  detonation  plume  within  the  first  few  milliseconds  after  detonation. 

The  net  result  of  the  Stage  1  model  is:  (1)  the  energy  release  to  the  plume  which,  in  turn,  determines  its 
initial  temperature  and  size,  (2)  the  initial  chemical  composition  of  the  detonation  plume,  and  (3)  the  mass 
and  size  distribution  of  entrained  soil  particles.  Collectively,  these  data  provide  the  initial  conditions  for 
the  treatment  of  buoyant  plume  rise,  afterburning  chemical  kinetics,  and  particle  deposition  in  Stage  2. 

With  respect  to  modeling  the  Stage  2  processes,  Figure  2  shows  that  while,  at  one  level,  it  is  possible  to 
distinguish  between  physical  processes  governing  the  plume  buoyant  rise  and  chemical  processes 
governing  gas  phase  afterburning  kinetics  and  gas  phase  species  interactions  with  small  entrained  soil 
particles  or  release  particulates  from  projectile  casings,  there  is  a  coupling  between  these  physical  and 
chemical  processes.  Enthalpy  changes  due  to  afterburning  kinetics,  for  example,  impact  the  plume 
temperature  and,  consequently,  the  mechanisms  influencing  the  puffs  rise  and  downstream  dispersion. 

The  basic  outputs  from  Stage  2  modeling  and,  in  essence,  the  SCM  model  as  a  whole,  are  the  chemical 
composition,  altitude,  downstream  distance,  and  size  of  the  plume  once  it  has  equilibrated  with  the 
atmosphere.  These  outputs  provide  both  a  basis  for  analysis  of  the  local  (fence  line)  environmental  impact 
of  detonation  emissions  as  well  as  inputs  to  further  processing  by  long  time  atmospheric  dispersion  (air 
quality  issues)  or  fate/transport  models  (ground/water  pollution  impact  issues). 
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3.0  SCM  Model  Description 


3.1  Detonation  Thermochemistry  and  Emission  Factors 

3.1.1  Phenomenology 

The  treatment  of  explosive  detonation  focuses  on  predicting  the  energy  release  and  other  thermodynamic 
properties  of  the  explosive  state  and  the  chemical  composition  of  the  mixture  of  detonation  products.  The 
thermodynamic  properties  of  the  detonation  state  play  a  fundamental  role  in  determining  the  temperature 
and  size  of  the  detonation  plume  while  the  chemical  composition  or  emission  factors  for  detonation 
products  provides  the  first  best  estimate  for  pollutants  released  in  the  blast.  Empirical  evidence  indicates 
that  the  detonation  energy  is  expended  in  a  volume  determined  by  the  equivalent  mass  of  TNT  and  that  this 
initial  fireball  quickly  expands  isentropically  to  near  the  ambient  pressure  [Boughton  and  DeLaurentis, 
1987].  Consequently,  the  initial  detonation  plume  temperature  and  volume  can  be  estimated 
thermochemically  given  the  chemical  composition  of  the  explosive  charge.  This  general  process  is 
illustrated  schematically  in  the  top  panel  of  Figure  5. 

However,  if  the  chemical  speciation  is  also  estimated  thermochemically,  there  are  significant  differences 
with  experimental  data.  Bomb  calorimetric  data  and  analyses  reported  by  Ornellas  [1982]  suggest  that  the 
chemistry  is  frozen  during  the  isentropic  expansion  resulting  in  a  detonation  plume  whose  chemical 
composition  is  not  in  thermochemical  equilibrium  at  the  temperature  and  pressure  of  the  isentropically 
expanded  detonation  plume.  The  study  by  Ornellas  [1982]  suggests  that  the  chemistry  freezes  at 
temperatures  around  2000  K  which  for  an  isentropic  expansion  corresponds  to  pressures  on  the  order  of 
10,000  atm.  Consequently,  the  SCM  model  estimates  the  initial  chemical  composition  of  the  detonation 
plume  from  the  high  detonation  pressure  mixture  by  assuming  that  the  expansion  is  significantly  more  rapid 
than  the  time  to  reach  chemical  equilibrium  and  that  chemical  speciation  freezes  out  at  temperatures  in  the 
range  reported  by  Ornellas  [1982].  The  panel  in  the  right  bottom  of  Figure  5  illustrates  this  for  an 
explosive  charge  consisting  of  Composition  B. 
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Thermochemical  Model  for  Detonation 


Chapman-Jouguet 
Detonation 

Explosive  Fireball 


Expansion 


Initial  Puff 


C-J  Detonation  Products  as  a 
Function  of  Explosive  Mass  Density 


CO 

co2 

ch4 

C2H4 

ch2o2 

NO 

nh3 


C-J  State  and  Adiabatic  Expansion 
Composition  B,  density  =1.7  g/cc 


pressure  (atm) 


■o. 

<t> 

3 


Figure  5.  Overview  of  the  thermochemical  detonation  model  showing:  (top)  the  evolution  from  the  initial 
explosive,  to  the  initial  fireball  following  detonation,  and  then  to  the  initial  puff  following  isentropic 
expansion  to  near  ambient  conditions;  (left  bottom)  an  example  of  detonation  emission  factors  for  selective 
species  and  their  variation  with  the  density  of  the  initial  explosive  charge;  and  (right  bottom)  profiles  for 
detonation  puff  temperature  and  mass  density  as  a  function  of  pressure  along  the  adiabat. 


3.1.2  Chemical  Equilibrium 


Detonation  pressures  are  on  the  order  of  several  tens  of  thousands  of  atmospheres  and,  as  noted,  finite  rate 
chemical  kinetics  in  the  expanding  detonation  puff  freezes  at  pressure  of  approximately  10,000  atm. 
Consequently,  non-ideal  gas  effects  may  play  an  important  role  in  determining  the  initial  speciation  of  a  HE 
detonation.  As  part  of  the  research  in  developing  the  SCM,  both  ideal  gas  and  non-ideal  gas  equilibrium 
codes  have  been  used  to  estimate  chemical  speciation.  Specifically,  the  Chemical  Equilibrium  with 
Applications  (CEA)  code  written  by  McBride  and  Gordon  [1996]  was  used  as  a  prototypical  ideal  gas 
equilibrium  code.  CHEETAH  developed  by  Fried  [1996]  and  based  on  the  earlier  TIGER  equilibrium  code 
[Cowperthwaite  and  Zwisler,  1973]  was  selected  as  the  non-ideal  gas  equilibrium  code  for  this  analysis. 
Although  CHEETAH  comes  bundled  with  four  product  libraries,  the  remainder  of  this  discussion  is  based 
on  the  BKW  library  which  is  a  computer-optimized  parameterization  of  the  BKW  equation  of  state  which 
was  noted  by  Fried  [1996]  to  give  good  results  for  explosives  containing  C,  H,  O,  N,  F,  and  Cl.  The  CEA 
codes  was  integrated  into  the  SCM  while  CHEETAH  functions  as  a  stand-alone  application  whose  outputs 
are  read  into  SCM. 


These  codes  were  used  to  explore  the  various  options  for  a  thermochemical  equilibrium  treatment  of  both 
the  C-J  state  properties  and  adiabatic  expansion.  Three  features  were  examined:  (1)  initial  explosive 
density,  (2)  C-J  state  detonation  products,  and  (3)  chemical  species  evolution  along  the  adiabat.  The  results 
of  that  study  are  discussed  below.  However,  in  the  discussion  that  follows  it  should  be  noted  that  energy 
losses  to  crater  formation,  soil  entrainment,  and  shell  fragmentation  are  ignored. 
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3. 1.2.1  Initial  Explosive  Mass  Density 


Within  the  context  of  a  non-ideal  gas  thermodynamic  equilibrium  treatment  the  detonation  products  depend 
on  the  initial  density  of  the  explosive.  The  left  bottom  panel  in  Figure  5  illustrates  this  dependence  through 
the  results  of  simulations  for  Composition  B  by  showing  the  emission  level  of  selected  species  in  units  of 
mole/kg-explosive  as  a  function  of  the  density  for  the  initial  explosive. 

3. 1.2. 2  C-J  State  Properties 

Table  1  lists  C-J  state  properties  and  mixing  ratios  for  the  major  detonation  products.  The  first  four 
columns  of  data  illustrate  the  variation  in  properties  and  composition  with  explosive  by  showing  data  for 
RDX,  TNT,  HMX,  and  Comp  B  (60%  RDX  /  39%  TNT  /  1%  wax).  This  data  was  generated  using  the 
CHEETAH  nonideal  gas  equilibrium  code  [Fried,  1996].  The  last  column  shows  data  for  Comp  B 
generated  using  the  CEA  ideal  gas  equilibrium  code  [McBride  and  Gordon,  1999].  Part  of  the  differences 
in  chemical  composition  for  the  last  (ideal  gas  equation  of  state)  set  of  data  is  due  to  the  lack  of  a  solid 
carbon  phase  (C  ).  The  resulting  increased  level  of  carbon  in  the  CEA  gas  phase  species  increases  the 
levels  of  CO,  CH4,  and  other  species  with  carbon  in  a  lower  oxidation  state  and  decreases  the  level  of  CO2. 


Table  1:  C-J  State  Properties  and  Major  Detonation  Products 


C-J  State  Properties  and  Composition 

CHEETAH  Nonideal  Thermochemical  Eq 

uilibrium 

CEA  Ideal 
Equilibrium 

RDX 

TNT 

HMX 

Comp  B 

Comp  B 

Density  (g/cc) 

1.8 

1.654 

1.9 

1.74581 

na 

Po  (atm) 

1 

1 

1 

1 

1 

Vo  (cc/g) 

0.55369 

0.60449 

0.52491 

0.57198 

na 

Eo  (cal/g) 

66.21027 

-78.4265 

60.52725 

14.14103 

14.16 

Pressure  (atm) 

341763.7 

191869 

388620.4 

282121.4 

290188.4 

Specific  Volume  (cc/g) 

0.4206 

0.454 

0.3995 

0.4301 

0.04256 

Temperature  (K) 

4191.9 

3685 

4113.6 

4035.6 

3834.9 

n2 

0.3332 

0.1365 

0.3332 

0.2542 

0.2588 

H20 

0.3331 

0.2266 

0.3332 

0.2903 

0.1468 

C02 

0.1525 

0.1303 

0.1584 

0.1444 

0.0738 

CO 

0.0280 

0.05797 

0.0162 

0.0382 

0.3579 

NO 

2.09E-04 

4.35E-05 

1.77E-04 

1.32E-04 

1.76E-04 

ch4 

1.65E-05 

1.35E-04 

4.35E-06 

6.03E-05 

1.460E-02 

c2h4 

1.24E-05 

5.75E-05 

3.02E-06 

2.58E-05 

9.01  E-04 

h2 

5.46E-06 

5.89E-05 

1.14E-06 

1.29E-05 

8.17E-02 

CH2OH 

3.76E-06 

1.51E-05 

1.16E-06 

7.1  IE-05 

2.15E-04 

02 

3.98E-05 

3.28E-06 

4.79E-05 

1.76E-05 

1.79E-06 

nh3 

1.57E-06 

1.13E-05 

3.99E-07 

3.62E-07 

1.55E-02 

ch2o 

7.61  E-08 

1.2E-06 

1.27E-08 

2.26E-10 

<  10-15 

c* 

0.1527 

0.4480 

0.1586 

0.2723 

<  10-15 
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3. 1.2.3  Equilibrium  Speciation  Along  C-J  Adiabat 


Figure  6  shows  species  mixing  ratios  for  expansion  along  the  C-J  adiabat.  In  Figure  6b,  the  equilibrium 
mixture  composition  was  obtained  using  the  ideal  gas  code  CEA.  Notable  for  the  nonideal  gas  equilibrium 
case  is  the  peak  level  for  organics,  methane,  ethene,  and  formaldehyde,  at  a  pressure  of  about  10,000  atm. 
This  pressure  corresponds  to  temperatures  between  1800  K  and  2200  K  depending  on  the  particular 
explosive.  Ornellas  [1982]  reports  that  chemical  kinetics  freeze  during  expansion  at  temperatures  between 
1800  K  and  2000  K.  Consequently,  the  estimated  levels  of  organics  based  on  chemical  equilibrium  with  a 
pre-specified  pressure  or  temperature  for  which  the  equilibrium  products  are  kinetically  frozen  is  sensitive 
to  the  point  at  which  the  equilibrium  products  are  frozen.  Additionally,  this  can  further  impact  other 
oxidation  products  from  afterburning  as  ambient  air  is  entrained  in  the  explosive  gas.  A  higher  level  of 
methane  for  example  increases  fuel  content  of  the  post-detonation  gas  resulting  in  a  higher  afterburning 
temperature. 


Nonideal  Gase  Expansion  Species  Profile 


o 
-+— > 
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c 
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Figure  6:  Species  mixing  ratios  for  expansion  along  the  C-J  adiabat  using  the  (a)  nonideal  equilibrium  code 
CHEETAFl  and  (b)  the  ideal  gas  equilibrium  code  CEA. 

3. 1.2. 4  Ideal  Gas  Versus  Non-ideal  Gas  Equations-of-State 

Based  on  the  discussion  in  Sections  3.1.2. 1  through  3. 1.2. 3,  it  is  clear  that  a  non-ideal  gas  thermochemical 
equilibrium  treatment  of  the  detonation  state  and  adiabatic  expansion  is  appropriate.  At  the  same  time, 
however,  the  differences  are  not  as  significant  as  might  be  expected  given  the  high  detonation  pressures. 
Rather,  the  largest  differences  may  be  associated  with  the  release  of  carbon  particulates  that  has  a  strong 
impact  on  the  gas  phase  levels  of  CO,  C02,  and  trace  organics.  Flowever,  this  is  assuming  as  reported  by 
Ornellas  [1982]  that  chemical  kinetics  freezes  at  temperatures  of  about  2000  K.  If  the  chemistry  for 
selected  groups  of  chemical  species  freezes  earlier  (higher  pressure  and  temperature)  along  the  adiabat, 
then  the  differences  between  an  ideal  gas  versus  non-ideal  gas  equation  of  state  equilibrium  solutions 
would  be  substantially  larger  as  shown  in  Figure  6. 

3.1.3  Initial  Cloud  Radius 

The  thermodynamic  properties  discussed  above  are  used  in  the  SCM  to  specify  the  initial  temperature, 
pressure,  and  volume  of  the  detonation  puff.  The  temperature  and  pressure  follow  directly  from  the 
adiabatic  expansion.  The  approach  adopted  to  specify  the  initial  puff  radius  is  given  below. 
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We  define  the  detonation  puff  as  the  initial  condition  when  this  fast  adiabatic  expansion  finishes  and  before 
ambient  air  and  soil  is  entrained  into  the  fireball.  The  kinetic  energy  of  the  projectile  is  negligible.  The 
fragmented  shell  parts  are  ejected  away  far  from  the  fireball.  Thus,  the  mass  inside  the  fireball  is  the  same 
as  the  amount  of  explosives.  (Note,  subsequent  treatment  of  the  dust  cloud  resulting  from  crater  formation 
in  Section  3.2  modifies  the  treatment  below  by  incorporating  the  mass  of  entrained  soil.) 


M0  =  PoV0 


(3.1.1) 


where  po  and  VQ  are  the  density  and  volume  of  the  initial  cloud,  respectively.  We  further  assume  that  all 
the  heat  released  from  the  detonation  is  used  in  increasing  the  fireball  enthalpy: 


Hr  M0CpT0 


(3.  1.2) 


where  Cp  and  T0  are  specific  heat  at  constant  pressure  and  initial  cloud  temperature,  respectively.  For 
simplicity,  we  assume  the  gases  in  the  fireball  follow  the  perfect  gas  law,  which  is: 

p  =  pRT  (3.1.3) 

where  p,  T,  p,  and  R  are  pressure,  temperature,  density  and  the  gas  constant,  respectively.  Using 
Equations  (3.1.1)  - 

fireball  volume  is 

y-ilHR 


Y 

Equations  (3.1.1)  -  (3.1.3)  and  relation,  c  =  — - — rR.  where  y  is  the  ratio  of  specific  heats,  the  initial 

p  y  —  1 


V  = 


Y 


(3.  1.4) 


The  initial  radius  of  the  fireball  (assumed  spherical)  can  be  readily  obtained  as: 


3  1 

Y-l 

|hr] 

4ti 

l  Y  , 

1  Pa  . 

(3.  1.5) 


If  we  use  pa  =  paRTa,  which  should  hold  for  ambient  air,  we  have: 


r 

o 


3  hr 

47i  p  c  T 
a  p  a 


1/3 


(3.  1.6) 


CL  is  a  coefficient  (smaller  than  1)  that  represents  all  the  neglected  energy  losses.  The  value  of  CL  should 
be  obtained  by  field  test  measurements  for  puff  size. 

As  an  example,  consider  the  81-mm  mortar  and  60-mm  mortar  with  explosive  formulations  given  by: 

81 -mm  mortar 

0.95  kg  of  Composition  B,  60%  RDX  and  39%  TNT 
0.57  kg  of  C4  charge,  91%  RDX 

60-mm  mortar 

0.43  kg  of  Composition  B,  60%  RDX  and  39%  TNT 
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In  this  calculation,  only  the  major  charges  are  included.  The  minor  charges  in  the  primer,  ignition 
cartridge,  propellant  charge,  booster,  and  fuse,  and  the  effect  of  wax  and  plasticizers  are  not  taken  into 
account.  Energy  loss  to  detonation,  crater  formation,  soil  entrainment,  and  shell  fragmentation  is  ignored. 
We  assume,  for  simplicity,  that  all  the  heat  released  from  the  detonation  goes  into  the  fireball.  The  heat 
generation  due  to  the  afterburning  is  also  neglected. 

The  heat  release  for  the  80-mm  and  60-mm  rounds  is  771 1  kJ  and  2232  kJ,  respectively.  Letting  the 
ambient  air  density,  air  temperature  and  heat  capacity  be  1.2  kg/m3,  270  °K,  and  1.15  kJ/kg/K,  respectively, 
then  based  on  the  approach  presented  above,  the  initial  cloud  radii  for  the  81 -mm  round  and  60-mm  round 
are  respectively,  1.6  m  and  1.1m.  The  initial  temperatures  of  the  two  rounds  are  obtained  as  441 1  °K  and 
4513  °K  respectively.  The  initial  densities  of  the  fireballs  are:  0.089  kg/m3  and  0.077  kg/m3,  respectively. 

3.2  Crater  Dynamics 

In  order  to  have  an  understanding  of  the  amount  of  soil  particles  eventually  entrained  into  the  puff,  we  need 
to  take  account  the  bomb  crater  dynamics  and  the  soil  characteristics.  The  crater  model  is  a  module  of  the 
Source  Characterization  Model  (SCM)  that  predicts  the  (1)  crater  mass,  (2)  dust  cloud  mass,  and  (3)  dust 
cloud  particle  size  distribution  resulting  from  an  explosion.  The  crater  module,  in  conjunction  with  the 
detonation  thermochemistry  module  discussed  in  Section  3.1,  is  used  to  establish  initial  conditions  for  the 
detonation  puff.  Figure  7  depicts  the  parameters  used  in  the  discussion  that  follows  to  describe  the  basic 
modeling  approach  adopted  in  SCM  to  describe  crater  formation  and  the  associated  dust  cloud. 


Figure  7.  Schematic  of  the  model  for  the  crater  formed  as  a  the  result  of  HE  detonation. 


3.2.1  Crater  Volume 


The  physics  of  a  ground  level  explosion  includes  the  chemistry  of  the  explosive  and  physical  interactions 
with  the  soil  and  atmosphere.  The  governing  differential  equations  are  well  understood  [Roddy  et  al., 

1977]  and  have  been  solved  numerically  to  give  reasonable  agreement  with  experimental  data  [Knowles 
and  Brode,  1977].  A  less  detailed  approach  has  been  developed,  based  partially  on  dimensional  analysis,  in 
which  empirical  correlations  were  developed  using  experimental  data  [U.S.  Army  Engineering,  1961].  The 
key  mathematical  relations  are  presented  below. 

The  crater  radius  rc  (units  in  m)  and  crater  depth  dc  (units  in  m)  are  given  by  the  following  relations: 


—  CIq  +  Cl^Ac  +  Q^AC  +  d2Ac 

(3.2.1) 

=  b0  +b\Ac  +b2Arc  +  b2A]  +b4A4c 

(3.2.2) 

13 


where  rc,  d,  .  and  Xc  are  the  scaled  crater  radius,  scaled  crater  depth,  and  scaled  explosion  depth 
respectively  given  by: 


W3 


d  c  — 


d 


/L  = 


W 

d 


0.3 


w- 


(3.2.3) 


(3.2.4) 


(3.2.5) 


Here  W  is  the  energy  release  of  the  explosive  charge  (units  in  equivalent  kg  of  TNT)  and  dh  is  the  distance 
above  ground  of  the  charge  (units  in  m)  (dh  is  negative  for  charge  below  ground).  The  coefficients  a,  and 
bh  where  i  denotes  vector  elements,  are  uniquely  determined  for  a  given  soil  type.  Soil  was  categorized 
into  the  seven  different  types  listed  in  Table  2.  The  corresponding  coefficient  parameters  a,-  and  h,  are 
given  in  Tables  3  and  4,  respectively. 


Table  2.  Description  of  Soil  Categories 


Category 

Description 

I 

Rock 

II 

Dry  cohesive  soils  (lightly  cemented) 

III 

Dry  sandy  soils 

IV 

Dry  to  moist  sand  and  frozen  ground 

V 

Soft  rock 

VI 

Wet  sand,  moist  cohesive  soils,  and  ice 

VII 

Wet  cohesive  soils  (not  saturated)  and  snow 

Table  3.  Coefficients  for  Scaled  Crater  Radius 


Soil 

a0 

ai 

a3 

I 

0.271 

-0.684 

0.390 

0.886 

II 

0.271 

-0.684 

0.390 

0.886 

III 

0.386 

-0.849 

0.367 

0.993 

IV 

0.503 

-0.954 

0.450 

1.190 

V 

0.629 

-1.080 

0.264 

1.120 

VI 

0.629 

-1.080 

0.264 

1.120 

VII 

0.806 

-1.280 

-0.178 

0.852 

Table  4.  Coefficients  for  Scaled  Crater  Depth 


Soil 

b0 

bi 

b2 

b3 

b4 

I 

0.113 

-0.477 

0.270 

1.84 

1.05 

II 

-0.571 

0.343 

2.24 

1.31 

III 

0.189 

-0.840 

0.447 

3.30 

2.10 

IV 

0.251 

-1.17 

0.494 

4.72 

3.34 

V 

0.189 

-0.840 

0.447 

3.30 

2.10 

VI 

0.331 

-1.49 

0.579 

4.92 

3.13 

VII 

0.449 

-1.82 

0.322 

4.11 

2.02 
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Scaled  crater  radius  and  scaled  crater  depth  are  plotted  in  Figures  8  and  9,  respectively.  All  curves  follow 
similar  patterns.  For  the  range  oflc  shown,  the  curves  lie  approximately  within  an  order  of  magnitude,  and 
they  increase  with  decreasing  scaled  burst  depth  up  to  a  critical  scaled  burst  depth  where  scaled  crater 
radius  and  depth  reach  a  maximum  and  decrease  thereafter.  For  a  given  scaled  burst  depth  the  scaled  crater 
radius  and  depth  increase  with  increasing  soil  category  number.  Recall  from  Table  2  that  the  soil  category 
numbers  increase  with  decreasing  soil  “hardness.”  Thus,  for  given  energy  release  and  charge  depth,  the 
crater  formed  will  be  larger  for  a  “softer”  soil. 


0.2  0.0  -0.2  -0.4  -0.6  -0.8 

Scaled  burst  depth  (m*kg  ) 


Figure  8.  Scaled  crater  radius  vs  scaled  burst  depth. 
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Scaled  Crater  Depth  vs  Scaled  Burst  Depth 


The  crater  formed  is  assumed  to  be  ellipsoidal  in  shape.  One  soil  characteristic  not  included  up  to  this 
point  is  vegetation  cover.  It  was  determined  in  a  previous  report  [Ebersole  et  al.,  1980]  that  the  crater 
volume  is  reduced  by  the  presence  of  vegetation  cover.  The  crater  volume,  including  the  effect  of 
vegetation  cover,  is  given  by: 


l 

3 


v 


(3.2.6) 


where  Vc  is  the  crater  volume,  which  contains  soil  only  (no  sod),  rc  and  d,  are  the  crater  radius  and  depth 
mentioned  earlier,  and  ds  is  the  sod  depth.  Equation  (3.2.6)  is  plotted  in  Figure  10.  In  the  figure,  the 
equation  has  been  normalized  such  that: 

K,o  =  |  ^cdc  (3.2.7) 

The  curve  is  parabolic  in  shape  with  a  minimum  at  djdc  =  1. 
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Volume  Reduction  Due  To  Vegetation 


Figure  10.  Volume  reduction  due  to  vegetation  cover.  Vc  0  is  crater  volume  without  vegetation  effect. 

3.2.2  Dust  Cloud  Mass 

The  primary  system  of  interest  for  the  SCM  model  is  the  dust  cloud  produced  by  the  explosion.  Of  the 
mass  associated  with  the  crater  volume  only  a  fraction  ends  up  in  the  initial  hot  air  mass  (or  puff)  resulting 
from  the  explosion.  In  this  section,  relations  for  predicting  the  total  mass  of  soil  particles  in  the  puff,  which 
were  developed  in  a  previous  report  [Pinnick  et  al.,  1983],  are  presented. 

The  mass  of  soil  particles  in  the  puff  Mdmt  is  given  by: 

dust  —  P soil  rl  case  V  sod  V  silt  *1  soil  *1  parittion  (3.2.8) 


where: 


Vsoil 


^1  compaction  ^1  lip^l  agglomeration 


(3.2.9) 


where  the  if  s  are  efficiency  factors  which  effectively  deplete  the  crater  volume  mass.  These  efficiencies 
are  presented  in  Table  5.  As  stated  in  the  report  by  Pinnick  et  al.  [1983],  “Values  indicated  with  a  question 
mark  show  uncertainties  due  to  insufficient  experimental  data.” 
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Table  5.  Crater  Volume  Depletion  Efficiencies. 


Symbol 

Description 

Value 

Given* 

Value 

Used 

n  casing 

Effect  of  artillery  shell  fragment 
kinetic  energy  on  available 
explosive  energy 

0.6  for  projectile 

1  for  bare  charge 

0.6  for  projectile 

1  for  bare  charge 

Hsod 

Dust  depletion  due  to  sod  cover 

1 

1 

Hsilt 

Fraction  of  small  (<  75mm) 
soil  particles  subject  to 
aerodynamic  arrest 

0.25-0.35 

0.3 

n  compaction 

Dust  depletion  due  to  compaction 
of  crater  soil 

0.6  (?) 

0.6 

hiip 

Dust  depletion  due  to  crater  ejecta 
deposition  on  crater  lip 

0.4  (?) 

0.4 

n  agglomeration 

Dust  depletion  due  to 
agglomeration  of  ejecta  under 
dynamic  (explosive)  loading 

0.1  -  0.5  (?) 

0.3 

n  partition 

Fraction  of  initial  dust  deposited  in 
buoyant  cloud 

0.2 -0.8  (?) 

0.5 

1— n  partition 

Fraction  of  dust  deposited  in 
non-buoyant  dust  skirt 

0.8  -  0.2  (?) 

0.5 

*  Question  marks  show  uncertainty  due  to  insufficient  data. 


Consider  a  sample  calculation  of  bare  charge.  The  mass  of  soil  particles  in  the  buoyant  cloud  is  calculated 
by: 


( 


V 


M 


dust 


M 


crater  J 


^casing  sod  silt  compaction  lip  agglomeration  partition 


(3.2.10) 


M 


dust 


crater  J 


1.08% 


Thus,  for  this  case  only  approximately  1%  of  the  crater  soil  particles  end  up  in  the  buoyant  cloud. 

3.2.3  Dust  Cloud  Particle  Distribution 


In  order  to  quantify  the  particle  size  distribution  experimental  measurements  are  needed.  We  have  found 
only  one  instance  where  a  size  distribution  was  measured  directly  on  an  explosion-produced  dust  cloud, 
reported  by  Pinnick  et  al.  [1983].  The  results  from  this  study,  which  are  presented  below,  are  used  to 
quantify  the  soil  particle  distribution  entrained  in  the  cloud. 
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This  group  set  off  explosions  on  four  occasions,  in  three  different  types  of  soils  and  with  varying  charge 
type  and  size.  They  measured  the  soil  distribution  approximately  30m  downwind  of  the  blast  (as  the  cloud 
passed  over).  The  three  soil  types  were  “clay”,  “sandy  clay”,  and  “sand  to  silty  sand”.  Their 
instrumentation  was  sensitive  to  soil  sample  particles  in  the  range  of  2  to  100  pm,  and  to  airborne  dust 
particles  in  the  range  of  0.2  to  100  pm.  Their  tests  showed  variations  in  the  soil  distributions  from  location 
to  location,  but  the  cloud  distributions  were  amazingly  similar  at  each  site.  All  measured  cloud  size 
distributions  could  be  roughly  fit  by  a  single  bimodal  lognormal  distribution. 

The  single  lognormal  number  distribution  is  given  by,  [following  Seinfeld  and  Pandis,  1998]: 
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N 


d N  _  N 

dD/;  42ttDp  In  (crj 


f 


exp 

V 


?  )  -  ln(z)  PgN ) 


|2  A 


21n2K) 


(3.2.11) 


where  N  is  the  total  number  of  particles  per  unit  volume,  Dp  is  particle  diameter,  D  pgN  is  the  geometric 
mean  or  median  of  the  particle  number  distribution,  and  ag  is  the  geometric  standard  deviation.  Pinnick’s 
results  are  summarized  in  Table  6. 

Table  6.  Particle  Size  Distribution  Parameters  Measured  by  Pinnick  et  al.  [1983], _ 


Small  mode 

Large  Mode 

N  (particles/cm  ) 

200 

0.07 

D  pgN  (microns) 

1.0 

45.0 

ffg 

2.6 

1.87 

A  bimodal  lognormal  distribution  if 

nN  iPp  )  =  nA r, small  ^P  p  ) 

>  simply  the  sum  of  two  lognormal  distributions. : 

+  '%,/arg  e{Dp)  (3.2.12) 

The  number  distribution,  using  the  parameters  in  Table  6,  is  plotted  in  Figure  11.  It  is  evident  that  the 
small  mode  dominates  the  particle  number  distribution.  This  result  can  be  somewhat  misleading.  The 
particle  number  distribution  is  useful  when  dealing  with  the  kinetics  of  the  system.  However  when  dealing 
with  the  mechanics  of  the  system,  the  mass  distribution  is  the  primary  distribution  of  interest. 
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Particle  Number  Distribution 


The  particle  mass  distribution  is  related  to  the  particle  number  distribution  by: 


dM  =  p—DldN 


•«(*>,)= 


(3.2.13) 


Combining  equations  (3.2.13)  and  (3.2.11)  gives: 


ttd;n 


(n  ) \_AM_ 

‘“K)X=6M®exp 


■[ln(Dj-ln(D«At)] 

21n2K) 


(3.2.14) 


As  it  stands  the  above  equation  is  not  in  the  usual  lognormal  form.  Rearranging  gives: 


•u(pPh 
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4lxDp  ln(cT) 


exp 


-1 

[lnl 

I-: 
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[p  pgM  )]2 

2  In2 1 

J 

1 

(3.2.15) 


where  M  is  the  total  mass  of  particles  per  unit  volume,  and  D  pgM  is  the  median  particle  diameter  for  the 
mass  distribution.  The  mass  distribution  has  the  same  geometric  standard  deviation  as  the  number 
distribution.  Mis  given  by: 
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(3.2.16) 


It  — 3 

M  =  pN  —  DpgM  exp 
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and  D pgM  is  given  by: 
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l  s)\ 

(3.2.17) 


Similar  to  the  number  distribution,  the  bimodal  lognormal  mass  distribution  is  the  sum  the  lognormal  mass 
distributions. 

nM^Pp^)=  nM,sma\\[Pp  )+  nM,  large  iP  p  )  (3.2.18) 

The  mass  distribution,  using  the  parameters  in  Table  6,  is  plotted  in  Figure  12.  The  bimodal  mass 
distribution  is  more  evenly  distributed  between  the  small  diameter  mode  and  the  large  diameter  mode  than 
was  seen  in  the  number  distribution. 


Particle  Mass  Distribution 


Pinnick  also  presented  values  of  total  particle  mass.  Here  the  results  are  reproduced  in  an  effort  to  provide 
a  self-consistency  check.  A  soil  density  of  p  =  2.5  g  cm"3  was  assumed.  Using  Equation  (3.2.16)  gives  M  = 
1.59E-5  kg  m"3  for  the  small  mode  and  M=  4.89E-5  kg  m"3  for  the  large  mode.  These  values  are  the  same 
as  those  reported  by  Pinnick.  The  total  number  and  mass  of  the  particle  distributions  are  summarized  in 
Table  7.  The  table  shows  that  Nsmau  »  Niarge  while  Miarge  >  Msmau.  This  result  reflects  on  the  relatively 
large  difference  in  median  diameters  for  the  two  modes. 
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Table  7.  Total  Particle  Number  and  Mass  Reported  by  Pinnick  [1983]. 


Small  mode 

Large  mode 

N  (particles  m'  ) 

2.00E+6 

0.07E+6 

M  (kg  in3) 

1.59E"5 

4.89E"5 

3.3  Puff  Fluid  Dynamics 

A  phenomenological  description  of  the  puff  rise  model  is  shown  in  Figure  2.  In  this  relatively  simple 
model,  the  puff  is  treated  as  a  spatially  homogeneous  sphere  of  mass  M  in  a  stratified  atmosphere.  The 
ambient  properties  are  uniform  in  a  horizontal  plane  for  a  given  altitude,  z,  with  the  wind  direction  dictating 
a  unique  horizontal,  x,  direction.  The  puff  then  rises  due  to  buoyancy  relative  to  the  atmosphere  and  moves 
horizontally  due  to  advection.  In  addition,  the  total  mass  of  the  puff  changes  as  ambient  air  is  entrained. 
Lastly,  the  enthalpy  of  the  evolving  puff  is  changing  due  to  a  number  of  factors  including  radiative  cooling, 
air  entrainment,  gas  phase  chemical  reactions,  and  gas-particle  phase  changes. 


Processes  Governing  Puff  Rise 


Figure  13.  A  conceptualized  view  of  the  model  for  buoyant  puff  rise. 


3.3.1  Analytical  Results 

A  buoyant  puff  rise  model  for  an  instantaneous  release  of  hot  materials  in  a  stable  atmosphere  with  cross 
wind  has  been  developed  [Zhang  et  al,  1996].  The  major  results  are  summarized  below.  The  model  applies 
the  conservation  laws  of  mass,  momentum,  and  energy  to  calculate  the  puff  motion.  The  entrained  air  into 
the  puff  is  assumed  to  be  proportional  to  the  puffs  outer  surface  area  and  the  velocity  difference  between 
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the  puff  and  the  surrounding  air.  Simple  analytical  formulas  for  the  puff  trajectory  and  final  equilibrium 
height  in  a  crosswind  and  stratified  atmosphere  are  obtained.  In  the  case  where  the  ambient  wind  is  absent, 
the  model  results  are  fully  validated  in  terms  of  entrainment  relation,  puff  rise  trajectory,  and  final  puff  rise. 

When  a  positively  buoyant  puff  is  released  into  the  atmosphere,  it  rises  due  to  its  buoyancy  and  its  size 
increases  due  to  entrainment.  In  a  stably  stratified  atmosphere,  the  puff  eventually  equilibrates  at  a  certain 
height. 

(a)  No  Ambient  Wind 


Figure  14.  A  schematic  diagram  of  puff  rise  in  the  absence  of  ambient  wind. 

In  the  absence  of  an  ambient  wind  the  puff  rises  vertically  as  shown  in  Figure  14.  We  assume  that  the  puff 
is  of  spherical  shape  with  radius  r(t)  and  height  z(t)  at  time  t.  The  puff  final  rise  is: 


zf  = 


3V0g' 


7i:a3N‘ 


1  /4 


(3.3.1) 


where  a  is  the  entrainment  coefficient;  V0  is  the  initially  released  pollutant  volume;  g'  =  g— — -  is  the 

ra 

reduced  gravity;  g  is  the  gravitational  acceleration;  pa  and  p0  are  the  ambient  air  density  and  the  puff 

density  at  its  releasing  time,  respectively.  N=  \  —  ^  ^Pa  is  the  so-called  Brunt- Vaisaila  frequency 

V  Pa  QZ 

characterizing  the  strength  of  the  ambient  stratification.  Its  value  can  be  calculated  by  using  : 


g  dpa_  g 
Pa  dz  Ta 


f<K+  S} 

dz  CP 


,  where  =  9.76  x  10  deg  K  /  m. 


(3.3.2) 
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In  the  initial  stage,  Nt  «1,  the  puff  radius  and  rise  trajectory  are: 


4  1 

t2  (3.3.3) 

I 

4  l 

t2  (3.3.4) 

When  t  is  small,  Nt  «  1,  the  trajectory  for  a  puff  in  stratified  atmosphere,  Equation  (3.3.4),  does  not 
involve  the  stratification  parameter.  This  is  because  the  stratification  effect  on  the  puff  rise  is  small 
initially. 

The  entrainment  relation  between  puff  radius  and  puff  rise  is  simply: 

r  =  az  (3.3.5) 

It  is  found  that  the  entrainment  coefficient  a  has  a  small  range  of  variation  from  0.26-0.30  (Zhang  et  al, 
1998). 

(b)  Large  Ambient  Wind,  P  =  uaN/g'  »  1 


3aV0g' 

271 

3V0g' 
2n  a3 


Figure  15.  A  schematic  diagram  of  puff  rise  in  the  presence  of  ambient  wind. 


This  case  occurs  when  the  ambient  wind  is  higher  or  the  stratification  is  strong  and  the  puff  buoyancy  is 
weak.  We  first  examine  the  puff  expansion  and  trajectory  when  Nt  «  n/2  and  the  dimensionless  parameter 
P  =  uaN/g'  »  1,  where  ua  is  the  horizontal  ambient  wind  velocity. 
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(3.5.6) 


(3.5.7) 
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The  maximum  puff  rise  in  the  presence  of  large  wind,  P  =  uaN/g'  »  1,  can  be  obtained  as: 


Ze=4- 

e  5 


32a3u3N5 


(3.5.8) 


Using  the  dimensionless  parameter  P  =  uaN/g'  and  denoting  the  maximum  puff  rise,  see  Equation  (3.3.1),  in 
the  absence  of  wind  as  zf,  we  can  rewrite  Equation  (3.3.8)  as  simply: 


|*  =  0.3  P“3/4 


(3.3.9) 


The  relation  between  puff  radius  and  downwind  distance  x  is: 
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(3.3.10) 


The  puff  rise  trajectory  formula  expressed  as  a  function  of  downwind  distance  x  is: 
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(3.3.11) 


This  relation  can  be  used  in  designing  the  proper  allocation  of  measurement  devices  in  performing  field 
experiments. 

(c)  Wind  Effect 

The  ratio  of  puff  rise  with  and  without  wind  is: 


Ze 


{ 


0  3p-3/4  P  >  i 

e- LVP  P<1 


(3.3.12) 


3.3.2  Governing  Equations 

Boughton  and  DeLaurentis  [1987]  developed  an  integral  model  of  puff  rise  from  HE  detonations.  This 
model  is  based  on  an  integration  of  the  partial  differential  conservation  equations  across  the  puff  cross 
section  assuming  a  hilly  axisymmetric,  spherical  puff.  It  includes  atmospheric  stratification  and  cross  flow, 
as  well  as  moisture  thermodynamics  assuming  an  equilibrium  mixture  of  dry  air  and  water.  Quantitative 
analysis  of  model  results  based  on  the  predicted  top  cloud  height  found  that  model  predictions  were  within 
25%  of  measurements  for  high  explosive  tests  using  TNT  and  Cyclotol  that  were  preformed  at  the  Tonopah 
Test  Range. 

The  SCM  puff  rise  model  was  developed  starting  with  the  treatment  given  by  Boughton  and  Delaurentis 
[1987].  The  model  was  then  extended  by  including  enthalpy  terms  for  gas  phase  chemistry  (Section  3.7) 
and  changes  in  total  mass  through  particle  deposition  (Section  3.5).  The  resulting  system  of  ODEs  are  then 
given  by: 
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(3.3.13a) 
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(3.3.13c) 
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(3.3. 13f) 


In  the  above  equations,  t  is  time,  M  is  the  puff  mass,  r  is  the  puff  radius,  g  is  the  gravitational  constant,  ( x , 
z)  are  the  horizontal  and  vertical  positions  relative  to  the  source,  (pr,  pz)  are  horizontal  and  vertical 
momenta,  and  pambient  and  PPuff& re  the  densities  of  the  ambient  atmosphere  and  puff,  respectively.  The 
change  in  momentum  in  the  horizontal  direction  depends  on  the  wind  speed  (yv),  while  the  change  in 
momentum  in  the  vertical  direction  depends  on  the  difference  in  density  between  the  puff  and  ambient 
environment.  The  change  in  mass  depends  on  the  entrainment  rate,  E,  and  ambient  density,  pambient-  hi 
eqaution  (3.3. 13f)  pressure  changes  are  not  explicitly  shown  since  they  are  negligible  relative  to  the  other 
terms. 


The  entrainment  rate  is  given  by: 


E 


4  7tr^  a 


V 


(3.3.14) 


where  r  is  the  puff  radius  and  |  Vr\  is  the  relative  velocity,  |  V  -  Vamb\ .  The  entrainment  coefficient,  a,  is  a 
semi-empirical  factor  that  is  determined  by  model  calibration  with  field  data. 

The  set  of  ODEs  for  puff  rise  depend  on  the  puff  temperature  and,  as  indicated  in  the  above  equations,  the 
enthalpy  of  the  puff  is  changing  in  response  to  a  number  of  factors,  including  adiabatic  cooling,  radiative 
heat  losses,  and  enthalpy  changes  induced  by  chemical  reactions. 
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To  solve  the  ODEs,  an  interface  was  developed  and  specifically  applied  to  the  Livermore  Solver  for 
Ordinary  Differential  Equations  (lsode),  which  can  be  used  to  solve  initial-value  problems  for  stiff,  and 
non-stiff  systems  for  first  order  ODEs. 

The  SCM  puff  rise  module  has  been  used  to  compute  the  puff  height  and  radius  as  a  function  of  time  for  an 
above  ground  detonation  of  a  Composition  B  mixture.  A  comparison  between  model  predictions  and 
available  observations  is  shown  in  Figure  16.  In  general,  the  model  predictions  where  found  to  agree  well 
with  the  reported  data. 

To  better  evaluate  the  results  shown  in  Figure  16,  the  sensitivity  of  model  predictions  for  the  detonation 
puff  height  and  radius  to  the  mass  of  the  explosive  charge  is  shown  in  Figure  17.  Model  predictions  for 
puff  height  were  found  to  be  more  sensitive  than  puff  radius. 

Observed  and  Calculated  Fireball  Height  and  Radius 
Above  Ground  Explosion:  Comp.  B  Mixture 
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Figure  16.  Observed  and  calculated  fireball  height  and  radius  for  a  composition  B  mixture. 


Sensitivity  to  Explosive  Mass 
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3.4  Buoyant  Puff  Flow  Structure 


3.4.1  Vortex  Method  Modeling 

To  understand  the  particle  deposition  process  from  a  buoyant  puff  (Section  3.5),  it  is  important  to  know  the 
flow  structure  inside  and  outside  of  the  buoyant  puff,  since  the  uprising  puff  exerts  an  upward  drag  force  on 
falling  particles.  Since  not  much  field  experimental  work  is  available  in  providing  the  information,  we  will 
rely  on  numerical  results  obtained  by  solving  Navier-Stokes  fluid  dynamics  equations. 

In  the  following,  we  will  show  some  numerical  results  to  illustrate  the  turbulent  flow  structure  in  a  buoyant 
puff.  For  details  regarding  the  modeling  techniques,  readers  are  referred  to  Zhang  and  Ghoniem  (1993)  for 
planar  two-dimensional  calculations  and  Issam  (2000)  for  axisymmetric  two-dimensional  calculations. 

Most  of  the  figures  in  this  section  are  from  the  planar  2-d  case. 

Zhang  and  Ghoniem  (1993)  developed  a  computational  model  for  the  rise  and  dispersion  of  a  wind-blown, 
buoyancy-driven  puff  in  a  stratified  atmosphere.  Lagrangian  numerical  techniques,  based  on  the  extension 
of  the  vortex  method  to  variable  density  flows,  are  used  to  solve  governing  equations.  The  puff  rise 
trajectory  and  the  dispersion  of  its  materials  in  the  crosswind  plane  are  obtained.  It  is  found  that  the 
computed  trajectory  agrees  well  with  the  two-thirds  power  law  of  a  buoyancy-dominated  puff.  The 
crosswind  section  of  a  wind  blown  buoyant  puff  is  similar  to  the  meridional  section  of  a  rising  puff  in  the 
absence  of  wind.  The  puff  cross-section  exhibits  a  kidney-shaped  pattern,  as  observed  in  many  laboratory 
and  field  experiments.  This  pattern  is  due  to  the  formation  of  two  counter-rotating  vortices  that  develop  as 
baroclinically  generated  vorticity  rolls  up  on  the  side  of  the  puff.  Figure  18  shows  the  evolution  of  the 
cross-section  of  a  puff  at  different  downwind  locations  (only  half  of  the  plane  is  shown  here  due  to 
symmetry). 


Figure  18.  The  evolution  of  the  cross-section  of  a  puff  shown  at  different  downwind  locations. 

The  flow  structure  is  dominated  by  a  large-scale  vortex.  Figure  19  shows  the  flow  structure  inside  and 
outside  of  the  puff  as  well  as  on  the  puff  external  surface.  Small-scale  eddy  structures  can  also  be  seen. 
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Figure  19.  (a)  The  buoyancy-induced  velocity  field  and  the  puff  outer  surface,  and  (b)  positions  and 
velocities  of  the  puff  outer  surface. 


z 


Figure  20.  A  schematic  diagram  showing  the  mechanism  of  (a)  large-eddy  generation  and  (b)  small-eddy 
generation. 


The  flow  field  and  the  vorticity  generation  pattern  is  illustrated  in  Figure  20,  which  is  a  schematic  diagram 
showing  the  baroclinic-vorticity  generation.  The  vorticity  generation  term  is  proportional  to  the  horizontal 
density  gradient.  For  a  puff  of  initially  elliptical  shape,  or  circular  shape,  the  large  scale,  counter-rotating, 
streamwise  vortices  are  formed  due  to  the  opposite  horizontal  density  gradient,  leading  to  a  large  scale 
rollup  which  changes  the  puff  shape  significantly  from  that  described  by  the  traditional  Gaussian  type 
concentration  distribution.  The  vortex  structure  induces  a  velocity  field  that  moves  the  puff  upwards. 

After  this  large-scale  rollup,  the  puff  external  surface  is  stretched  and  distorted,  secondary  vortices,  with 
the  same  or  the  opposite  sign  of  vorticity,  are  generated  and  the  puff  motion  becomes  turbulent  as  shown  in 
Figure  20  (b).  After  the  large-scale  rollup,  both  positive  and  negative  circulations  are  generated,  but  the 
total  circulation  remains  constant. 

Figure  21  shows  the  effect  of  Reynolds  number  on  the  puff  shape.  For  large  puff  with  strong  buoyancy,  the 
puff  motion  is  turbulent;  for  small  puff  with  weak  buoyancy,  the  motion  is  laminar.  However,  the  puff 
sectional  shape  is  always  dominated  by  the  large-scale  vorticity. 


Figure  21.  Puff  cross-sections  shown  for  different  buoyant  Reynolds  numbers. 

In  summary,  the  numerical  simulation  shows  the  following  characteristics  of  the  puff  motion,  which  will  be 
used  in  the  following  sections: 

•  A  large-scale  vortex  ring  dominates  the  puff  motion. 

•  The  circulation  of  the  vortex  ring  keeps  constant  as  the  puff  rises. 
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The  value  of  this  circulation  is  proportional  to  the  horizontal  density  gradient  integrated  over  the  vertical 
extent  of  the  puff. 


3.4.2  Extending  Hill’s  Vortex  Modeling 


3.4.2. 1  Hill’s  Spherical  Vortex 


For  detailed  derivation  of  Hill’s  vortex  flow  structure,  the  readers  are  referred  to  Article  165  of  Lamb 
[1945].  Denoting  r,  z  as  the  horizontal  and  vertical  distance  from  the  vortex  center,  the  stream  function  of 
the  flow  field  of  a  Hill’s  vortex  can  be  expressed  as: 
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where  a  is  the  radius  of  the  spherical  vortex,  R  =  »J r2  +  Z2,  and  W  is  the  upward  self-induced  vertical 
velocity. 

Spherical  Vortex  -  Streamlines 


Figure  22.  Streamlines  of  Hill’s  spherical  vortex. 

The  velocity  field  induced  by  this  vortex,  in  terms  of  vertical  component  Uz  and  horizontal  component  Ur 
can  be  described  by: 

Vertical  component: 


u  =  -  w 
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(3.4.2b) 


Horizontal  component: 
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The  velocity  field  is  shown  in  Figure  23. 


(3.4.3a) 


(3.4.3b) 


Figure  23.  Velocity  field  of  a  Hill’s  spherical  vortex. 

The  strength  of  the  vortex  ring  is  related  to  vertical  velocity  of  the  spherical  vortex  as: 

r  =  -jTtUa  (3.4.4) 

3. 4.2. 2  Extended  Hill’s  Vortex  Model 

Comparing  Figure  19  and  Figure  22,  we  find  that  the  flow  field  of  a  buoyant  puff  may  be  simplified  by 
using  a  Hill’s  spherical  vortex.  The  buoyant  puff  generates  a  vortex  ring  in  the  air  and  its  rising  motion  can 
be  viewed  as  driven  by  a  vortex  ring.  The  vorticity  field  in  the  meridional  section  is  generated  due  to 
barolicinc  vorticity  generation  mechanism.  Since  the  gravity  is  vertical,  the  chief  generation  term  is  the 
horizontal  density  difference  integrated  over  the  vertical  puff  extent.  Scaling  analysis  shows  that  the 
strength  (circulation)  of  the  vortex  ring,  F,  follows  the  relation  below: 


r  oc  Buoyant  Velocity  x  Vertical  extent  oc  V  (PapaPo)  g  2  Ro  2Ro  V^2  g'  v0  (3.4.5) 
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where:  pa  and  p0  are  the  ambient  air  density  and  the  initial  puff  density,  respectively,  R0  is  the  initial  puff 


radius,  g'  =  g- 


is  the  reduced  gravity,  and  V0  is  the  volume  of  the  initial  puff.  There  are  many 


small  eddies  generated  as  the  puff  move  upwards,  but  the  motion  of  the  bulk  part  can  be  modeled  by  this 
large  scale  vortex  ring. 


The  value  of  this  circulation  remains  the  same  even  when  ambient  air  is  entrained  into  the  puff. 
Consequently,  Hill’s  spherical  vortex  is  used  to  model  the  flow  structure  inside  and  outside  of  the  puff  by 
maintaining  the  circulation.  Since  the  size  of  a  Hill’s  vortex  does  not  change,  in  order  to  incorporate  the 
entrainment  process  for  an  expanding  buoyant  puff,  the  Hill’s  vortex  concept  is  extended  by  assuming  the 
vortex  size  increases  as  the  puff  rises.  As  an  example,  a  simple  but  frequently  used  relation,  that  the 
entrainment  velocity  is  linearly  proportional  to  the  puff  rising  velocity,  gives: 


dR  w 

—  =  aw 

dt 


(3.4.6) 


where  a  is  a  constant.  Denoting  R0  and  z0  as  the  initial  puff  radius  and  initial  puff  height,  the  integration  of 
equation  (3.4.6)  yields: 

R  =  R0  +  a(z  -  z0)  (3.4.7) 


We  assume  the  size  of  the  extended  spherical  vortex  increases  with  height  as  equation  (3.4.7)  indicates. 

In  this  work,  we  have  assumed  that  the  puff  is  spherical  and  the  initial  puff  height  is  the  same  as  the  initial 
puff  radius.  We  assume  the  relation  given  by  Equation  (3.4.7)  still  holds  so  that  Equation  (3.4.6)  gives: 
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(3.4.8) 


Noticing  that  the  puff  radius  becomes  much  larger  than  the  initial  puff  radius  very  quickly,  equation  (3.4.8) 
simplifies  to: 
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(3.4.9) 


Therefore,  as  the  puff  expands,  the  vortex  circulation  keeps  the  same  and  vertical  velocity  decreases  as  the 
puff  rises. 


3. 4. 2. 3  The  rise  of  extended  Hill’s  vortex 


Very  close  to  the  source  location,  the  effect  of  ambient  wind  and  stratification  can  be  ignored.  Then,  the 
height  of  a  buoyant  puff  is  proportional  to  the  square  root  of  the  time  elapsed. 
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In  terms  of  the  puff  height,  the  vertical  velocity  can  be  obtained  as: 
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Comparing  the  (3.4.9)  and  (3.4.10),  we  have: 
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3.5  Particle  Deposition 


The  characterization  of  particle  deposition  from  a  buoyant  puff  generated  from  explosive  detonations  is  one 
of  the  key  elements  to  successful  modeling  of  the  environmental  impact  of  munitions  training  activities. 
There  are  very  few  field  tests  that  address  particle  deposition.  This  is  mainly  because  the  initial  puff 
temperature  is  very  high,  the  time  duration  is  short,  the  variability  is  high,  and  very  little  understanding  of 
the  physical  and  chemical  processes  is  available.  In  this  work,  we  will  approach  this  problem  via 
mathematical  modeling  based  on  physical  understanding  as  well  as  available  field  and  laboratory 
measurements.  The  quantified  results  can  be  used  to  design  future  tests  and  model  the  emitted  pollutants 
for  air  quality  and  ground  water  quality  analysis. 

The  particulate  pollutants  emitted  from  a  detonation  can  be  divided  into  two  broad  categories  based  on 
particle  size.  The  first  category  includes  particles  with  diameters  less  than  about  50  pm  that  disperse  with 
gaseous  emissions  as  the  buoyant  puff  evolves  in  the  atmosphere.  The  second  category  includes  larger 
particles  (>5()liiti)  that  fall  to  ground  due  to  gravity.  The  temporal  and  spatial  deposition  pattern  of 
particles  is  important  in  estimating  how  much  pollutants  in  the  form  of  particles  are  deposited  to  the 
ground.  Since  the  particle  dropout  may  change  the  puff  mass  loading  and  the  thermal  energy  content 
significantly,  the  deposition  process  is  also  important  to  the  estimation  of  the  final  puff  rise  that  is  a  critical 
parameter  in  the  air  quality  dispersion  modeling, 

The  larger  particles  in  a  buoyant  puff  are  composed  of  soil  particles  entrained  during  detonation, 
unexploded  explosives,  and  also  fragmented  shell  parts.  During  the  deposition  process  of  these  particles, 
they  scavenge  pollutants  in  the  puff,  absorb  heat  energy,  and  may  even  locally  quench  the  afterburning 
combustion  in  the  puff.  It  is  important  to  understand  the  scavenging  effect  of  the  entrained  soil.  To  model 
these  scavenging  processes,  it  is  important  to  know  the  residence  time  of  a  particle  inside  the  hot  buoyant 
puff. 

The  particle  residence  time  depends  on:  the  particle  size  and  density  that  determines  the  gravitational  force; 
the  puff  buoyancy  that  determines  the  uplifting  drag  force  acting  on  the  particles;  and  the  local  atmospheric 
conditions  that  affect  diffusion  and  also  the  afterburning  chemistry. 

3.5.1  Deposition 

We  will  first  consider  the  deposition  in  the  absence  of  puff  buoyancy.  A  puff  produced  by  the  release  of 
detonated  materials  may  consist  of  a  mixture  of  gases  and  particles.  Larger  particles  may  fall  to  the  ground 
due  to  gravitational  settling.  Smaller  particles  and  gases  may  deposit  on  the  ground  surface  due  to  dry 
deposition,  which  is  caused  by  a  combination  of  chemical  reaction  and  physical  interception  by  the  ground 
and  vegetation.  In  the  presence  of  rain,  fog,  or  snow,  the  pollutant  may  be  removed  from  the  puff  and 
deposited  on  the  ground  either  by  absorption  or  collection  by  the  water  droplets  or  snowflakes.  Sehmel 
[1984]  and  Slinn  [1984]  give  thorough  reviews  of  the  topics  of  dry  deposition  and  wet  scavenging, 
respectively.  For  this  project,  we  considered  deposition  due  solely  to  gravitational  settling. 

In  the  SCM,  Stoke's  law  is  used  to  define  the  gravitational  settling  velocity,  Vg;  as: 


(ps  -  pJgDpC 

18r| 


(3.5.1) 


where,  ps  is  the  particle  density  (g/m2);  pajr  is  the  air  density  (g/m2);  g  is  the  gravitational  acceleration 

(9.81  m/s2);  Dp  is  the  particle  diameter  (m);  and  q  is  air  viscosity  (1.81  x  10'^  kg/m/s).  Since  particles 
entrained  in  the  puff  (e.g.  soil  particles)  span  a  wide  size  range,  Equation  3.5.1  includes  the  Cunningham 
correction,  C,  for  small  particles  given  by: 
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C  =  i  +  2k- 

Dn 


(3.5.2) 


a!  +  a2  exp(-^j 


where  X  is  the  mean  free  path  of  air  molecules  (6.53  x  10'^  cm),  and  a\,  a2,  ay  are  constants  (1.257,  0.40, 
0.55,  respectively). 

3.5.2  Particle  Motion  in  Buoyant  Puff 

Figure  24  shows  a  schematic  diagram  on  buoyant  puff  rise  and  the  particle-laden  stem  beneath.  It  is 
important  to  know,  given  particle  size  and  puff  buoyancy,  when  the  particles  will  fall  off  from  the  puff. 


Figure  24.  Buoyant  puff  rise  and  particle-laden  stem  formation. 

The  motion  of  a  particle  is  governed  by  its  gravity  as  well  as  by  the  drag  force  due  to  the  ambient  air 
motion.  The  two  components  of  the  governing  equations  are: 


=67nirp(uf-u) 

(3.5.3) 

=  67rnrp(wf  -  w)  -  4^(ps  -  pa)g 

(3.5.4) 

where  ps  is  the  particle  average  density,  rp  is  the  radius  of  the  particle,  u,  w  are  the  horizontal  and  vertical 
components  of  the  particle  velocity,  q  is  the  air  viscosity,  and  uf  and  wt-  are  the  horizontal  and  vertical 
velocity  components  of  the  ambient  air  flow.  Equations  (4.2)  and  (4.3)  can  be  used  for  Uf  and  Wf  for 
buoyancy  generated  large-scale  puff  motion.  These  two  equations  can  be  simplified  to: 
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where  dp  is  the  particle  diameter  and  g" 


(Ps  -  Pa)  g 

P^  ' 

As  an  example,  we  applied  Equations  (3.5.5)  and  (3.5.6)  for  the  deposition  of  particles  from  a  buoyant  puff 
generated  from  detonation.  The  initial  puff  is  released  at  ground  with  radius  1  m  and  density  0.8  kg/m3. 
The  ambient  air  density  is  1.2  kg/m3  with  viscosity  1.81  x  105  kg/m/s.  The  density  of  the  soil  particles  is 
2650kg/m3.  We  assume  the  entrainment  rate  has  a  value  of  0.3.  The  following  discussion  specifically 
focuses  on  larger  particles  although,  as  noted  above,  SCM  treats  a  wide  size  range  including  small  particles 
that  remain  in  the  puff  and  are  carried  downstream  without  significant  fall-out. 

Figure  25  shows  the  rise  and  fall  of  particles  with  diameter  100  pm  in  the  buoyant  puff.  The  initial  particle 
cloud  has  height  2  m  and  width  10m. 

The  particle  cloud  is  generated  due  to  the  ejection  of  crater  material  and  to  the  uplifting  of  soil  particles 
during  the  shaking  of  the  ground  by  the  detonation.  The  whole  rise  and  fall  process  can  be  divided  into  six 
stages  as  analyzed  below. 

Stage  1,  around  t  =  1  second,  the  buoyant  puff  rises  from  the  dust  cloud.  Some  particles  are  entrained  into 
the  puff  moving  upward  and  the  rest  fall  down  due  to  gravity. 

Stage  2,  around  t  =  3  seconds,  the  buoyant  puff  separates  from  the  dust  cloud.  The  number  of  particles 
moving  up  with  the  puff  is  more  than  the  number  of  particles  initially  inside  the  puff.  This  is  due  to  the 
effect  of  entrainment.  When  the  puff  entrains  ambient  air,  it  also  entrains  the  soil  particles.  At  this  time, 
the  dust  cloud  falls  to  the  ground  almost  completely. 

Stage  3,  around  t  =  10  seconds,  the  particle-laden  buoyant  puff  rises.  The  distribution  of  particles  in  the 
puff  is  not  uniform.  Most  of  the  particles  are  associated  with  the  vortex  ring.  The  kidney  shaped  puff 
structure  is  apparent.  This  is  consistent  with  the  numerical  simulation  results  discussed  in  Section  3.4. 

Stage  4,  around  t  =  20  seconds,  the  particles  in  the  puff  start  falling.  This  is  when  the  buoyant  puff  entrains 
much  ambient  air  and  the  flow  structure  is  not  strong  enough  to  hold  the  particles  moving  with  the  puff. 

The  long  maintained  coherent  kidney  shaped  structure  starts  to  change. 

Stage  5,  around  t  =  30  seconds,  a  particle-laden  stem  is  formed  below  the  bottom  of  the  buoyant  puff  and  is 
stretching  downwards.  At  this  time,  the  top  part  of  the  stem  still  rises  since  it  is  closer  to  the  vortex  flow 
structure.  The  lower  part  falls  down  due  to  gravity.  We  can  see  that  the  particles  at  the  vortex  ring  region 
at  t=20  seconds  are  first  stretched  forming  the  lower  part  of  the  stem.  The  particles  outside  of  the  vortex 
ring  region  at  t=20  seconds  are  stretched  later  forming  the  top  part  of  the  stem. 

Stage  6,  all  the  particles  have  fallen  outside  of  the  buoyant  puff.  The  top  part  of  the  stem  may  or  may  not 
be  attached  to  the  puff.  The  lower  part  of  the  stem  is  reaching  to  the  ground. 
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Figure  25.  The  rise  and  fall  of  particles  with  diameter  100  pm. 
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Figure  26  shows  the  rise  and  fall  of  particles  with  diameters  ranging  from  20pm  to  120  pm.  The  particles 
with  diameter  20pm  near  the  ground  are  the  ones  associated  with  the  initial  dust  cloud.  They  remain  near 
the  ground  because  of  their  smaller  size  and  lower  gravitational  settling  velocity  as  expected  from  Equation 
(3.5.1).  This  has  been  observed  in  field  observations. 

Around  t  =  15  seconds,  the  main  particles  in  the  stem  are  of  size  120pm.  The  smaller  particles  are  either  in 
the  buoyant  puff  or  near  the  ground. 

Around  t  =  30  seconds,  the  main  particles  in  the  stem  are  of  size  100pm.  The  smaller  particles  are  still  in 
the  puff  and  the  particles  with  size  120  pm  have  deposited  to  the  ground  completely. 

Around  t  =  45  seconds,  the  puff  s  size  keeps  increasing  because  of  the  entrainment  of  the  ambient  air.  The 
particles  with  size  100  pm  have  detached  from  the  buoyant  puff.  The  kidney  shaped  pattern  of  particles 
with  size  80pm  start  to  deform. 

Around  t  =  60  seconds,  the  particles  with  size  100  pm  have  deposited  to  the  ground  completely.  The 
particles  with  size  80pm  form  the  top  part  of  the  stem. 

Figure  26  indicates  that  heavier  particles  fall  down  quickly  and  are  deposited  close  to  the  source  location  in 
the  absence  of  wind.  The  stem  is  composed  of  particles  with  different  sizes  at  different  times,  heavier 
particles  earlier  and  lighter  particles  later. 

Figure  27  shows  the  number  of  particles  in  the  buoyant  puff  as  a  function  of  time  for  seven  cases  with 
different  particle  diameter.  The  number  of  particles  in  the  initial  buoyant  puff  is  about  760  for  all  cases  see 
Figure  27b.  As  the  buoyant  puff  rises,  it  entrains  particles  into  the  puff.  More  light  particles  are  entrained 
into  the  puff  than  heavier  particles.  The  maximum  number  of  20  pm  particle  is  about  1000  at  about  t  =  1.5 
seconds  and  the  maximum  number  of  140  pm  particles  is  about  860  at  about  t  =  0.5  seconds.  See  the 
description  of  stage  2  on  Figure  25.  Immediately  after  the  number  of  particles  in  the  buoyant  puff  reaches 
maximum,  there  is  a  drop  off  of  the  particles  from  the  puff.  This  is  due  to  the  separation  of  the  buoyant 
puff  from  the  ground  based  dust  cloud. 
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Figure  26.  The  rise  and  fall  of  particles  of  various  diameters  in  a  buoyant  puff:  red-20pm,  green-40pm, 
blue-60|im,  cyan-80  m,  yellow-100  m,  and  puiple-120  m. 
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Entrainment  of  Particles  into  Buoyant  Puff 


Figure  27.  Number  of  particles  in  the  buoyant  puff  at  different  times:  (a)  for  duration  up  to  100  seconds, 
(b)  for  the  first  10  seconds. 

After  the  separation,  the  number  of  particles  in  the  puff  is  about  540  for  20pm  particles,  480  for  40pm 
particles,  and  430  for  60pm  particles,  etc.  These  numbers  of  particles  can  be  sustained  in  the  puff  for  a 
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period  of  time  before  they  fall  off:  we  call  this  time  the  particle  residence  time.  The  residence  times  for 
smaller  particles  are  longer  than  100  seconds.  The  residence  times  for  80pm,  100pm,  and  120pm  particles 
are  about  80  seconds,  32  seconds  and  15  seconds,  respectively.  The  140  pm  particles  fall  off  quickly  after 
the  separation  process.  There  is  not  a  significant  residence  time  period. 

It  should  be  pointed  out  that  the  results  presented  in  Figure  25  through  27  are  for  the  example  case.  The 
timing  of  particle  fall  off  depends  on  the  buoyant  puff  conditions  and  ambient  atmospheric  conditions.  For 
the  ease  of  application,  we  derive  the  characteristics  particle  residence  time  and  particle  fall  off  height 
below.  Their  applicability  awaits  the  validation  by  field  experimental  measurements. 

When  the  buoyant  puff  rises,  the  puff  material  and  ambient  air  exerts  an  upward  drag  force  to  the  particle. 
As  the  puff  rises,  it  entrains  more  air  and  the  upward  motion  slows  down.  Eventually,  the  gravitational 
force  on  the  particle  exceeds  the  drag  force  and  the  particle  falls  down. 

At  the  height  that  the  particle  begins  to  fall,  the  right  hand  side  of  equation  (3.5.6)  vanishes,  and  w  =  0  and 
wf=  W: 


1 8r)8 
's  p 


g  =  g 


(3.5.7) 


Substituting  equation  (3.4.9)  and  (3.4.12)  into  (3.5.7),  the  particle  fallout  height  is  then  given  by: 
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and  the  particle  residence  time,  which  characterizes  how  long  a  particle  stays  in  the  puff,  can  be  obtained 
from  equations  (3.4.10)  and  (3.5.8): 
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«  g,  equation  (3.5.8)  and  (3.5.9)  can  be  further  simplified  to: 

zd  = 

9r| 

Psd2gV 

/  3Vog' 

2na3 

(3.5.10) 

t  H  = 

81r|2 

/  3v-g; 

(3.5.11) 

d 

P2dpg2  ^ 

/  2na 

3.5.3  Effect  of  Particle  Fall-off  on  Buoyant  Puff  Rise 

When  heavier  particles  fall  from  the  buoyant  puff  (Figure  28),  the  total  mass  and  heat  content  of  the  puff 
are  decreased.  These  effects  may  significantly  affect  the  final  puff  rise.  We  will  formulate  a  model  for  the 
buoyant  puff  rise  with  significant  particle  fall  off  in  this  chapter.  We  assume: 

•  Atmospheric  stratification  is  not  important. 

•  Ambient  wind  is  not  important. 

•  The  buoyant  cloud  is  spherical. 

•  Soil  particles  are  of  same  size  (to  be  extended). 

•  Initial  condition:  t=0,  z=0,  w=0 


41 


Gas  puff  with 
small  particles 


Figure  28.  A  schematic  diagram  on  the  buoyant  puff  rise  and  particles  fall  off. 

The  mass  conservation  says  that  the  rate  of  mass  change  in  the  puff  is  the  difference  between  the  rate  of  air 
entrained  into  the  air  and  the  rate  of  particle  mass  fall  out  of  the  puff: 

^-(Pp^fr3)  =4m-2awpa  -  pd7tr2Vd  (3.5.12) 

where  p/;  is  the  averaged  puff  density,  r  is  the  puff  radius,  a  is  the  entrainment  rate,  w  the  puff  rising 
velocity,  pa  is  the  ambient  air  density,  pd  is  the  mass  of  particles  divided  by  the  puff  volume,  and  Vd  is  the 
particle  fall  off  velocity.  Assuming  there  are  Nd  particles  in  the  puff,  denoting  ps  and  rs  as  the  particle 
material  density  and  particle  radius  respectively,  by  definition,  we  have: 

Pd^r3  =  NdPs^r2  (3.5.13) 

Newton’s  second  law  in  the  vertical  direction  states  that  the  rate  of  change  in  vertical  momentum  of  the 
puff  is  the  same  as  the  net  force  in  the  vertical  direction  acting  on  the  puff. 

-|-(PpXr3w)  =  ^r3g(Pa-pp)-Pd™2(w  +  Vd)w  -  37n)rs(w  +  Vd)Nd  /  C  (3.5.14) 

On  the  right  hand  side  of  equation  (3.5.3),  the  first  term  is  the  buoyant  force  owing  to  the  fact  that  the  puff 
density  is  lower  than  the  ambient  air  density,  the  second  term  is  the  rate  of  vertical  momentum  loss  due  to 
the  fall  out  of  particles  from  the  puff,  the  third  term  is  the  viscous  frictional  force  acting  on  the  puff  due  to 
the  downward  motion  of  particles  in  the  puff.  C  is  Cunningham  correction  factor  see  equation  (3.5. 13). 
Applying  Newton’s  second  law  to  the  group  of  particle,  we  have: 

-  -^(Ndps^r3Vd)  =  -  Ndps4^r3g  +  37n-|rs(w  +  Vd)Nd  /  C  (3.5.15) 

The  rate  of  change  in  vertical  momentum  of  particles  is  due  to  the  gravitational  force  acting  on  the  particles 
and  the  frictional  drag  force.  In  equations  (3.5.3)  and  (3.5.4),  the  relative  velocity  between  the  upward  puff 
motion  and  the  downward  particle  motion  should  be  calculated  separately  for  cases  particles  are  inside  and 
outside  of  the  puff  as  below: 
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(3.5.16) 


w  +  Vd 


(  w  +  Vd,  if  z  -  zd  <  r 
t  Vd,  if  z  -  zd  >  r 


where  z  is  the  height  of  the  buoyant  puff  center,  zd  is  the  center  height  of  the  falling  particles,  and  r  is  the 
radius  of  the  buoyant  puff.  Kinematics  relations  also  hold: 

#  =  w  (3.5.17) 


dzd 

dt 


-vd 


(3.5.18) 


The  last  equation  is  obtained  from  the  law  of  energy  conservation.  A  similar  equation  has  been  derived  for 
the  rise  and  dispersion  of  reacting  puff  (Zhang  et  al,  1997): 

■^:(Pp^r3Hp)  =  HEntr  +  HReact  +  HRad  +  Hxhermal  +  HGround  +  HDepo  (3.5.19) 


where: 


Hp  =  CpTp  (3.5.20) 

is  the  specific  heat  of  the  puff.  Equation  (3.5.19)  states:  The  rate  of  enthalpy  change  in  the  puff  equals  to 
the  sum  of  the  following  terms: 

HEntr  =  47tr2awcpTa  (3.5.21) 

is  the  entrainment  rate  of  the  ambient  enthalpy. 

HRCact  is  the  heat  generation  rate  due  to  chemical  reactions.  HRad  is  the  rate  of  heat  loss  due  to  thermal 
radiation  to  the  surroundings.  HThermal  is  the  rate  of  heat  change  due  to  phase  changes.  HGround  is  the  heat 
transfer  from  the  hot  puff  to  the  cold  ground  when  the  puff  touches  the  ground. 

HDepo  =  NdPd^d3  CpTs  (3.5.22) 

is  the  gravitational  fall  out  rate  of  enthalpy  associated  with  particles,  where  Ts  is  the  mean  particle 
temperature.  This  temperature  can  be  obtained  by  considering  heat  transfer  from  the  hot  puff  to  the  cold 
soil  particles. 


3.6  Meteorology 

Reliable  predictions  for  the  puff  downwind  distance  and  vertical  height  as  a  function  of  time  require 
solving  the  ODEs  for  puff  motion  in  a  stratified  atmosphere.  In  the  SCM,  a  meteorology  module  has  been 
implemented  which  extracts  ambient  temperature,  pressure,  wind  speed,  wind  direction,  and  relative 
humidity  from  user-supplied  profiles.  The  meteorology  module  also  includes  algorithms  for  calculating 
other  measures  of  atmospheric  stability  and  turbulence  including  the  friction  velocity,  Monin-Obukhov 
length,  sensible  heat  flux,  and  adiabatic  lapse  rate. 

3.6.1  Atmospheric  Stratification 

The  puff  rise  dynamics  depends  on  the  variation  of  atmospheric  properties  with  altitude.  This  includes  the 
wind  speed,  temperature,  and  relative  humidity. 
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Letting  z  represent  the  vertical  height  in  the  atmosphere,  then  the  pressure  profile  is  given  by:: 


*& =-/*)* 

dz 


(3.6.1) 


where  g  is  the  acceleration  due  to  gravity  and  the  ideal  gas  density  is: 


P(z) 


p{z)M 

RT{z) 


(3.6.2) 


where  p(z)  is  the  pressure,  T(z)  is  the  temperature,  and  M  is  the  molecular  weight.  Assuming  the  pressure 
at  the  ground  is  known  and  that  the  temperature  as  a  function  of  height  is  also  know  for  a  few  points,  then 
the  pressure  profile  can  be  obtained  by  combining  Equations  3.6.1  and  3.6.2  and  integrating: 


p{z)=  p{0)e 


(3.6.3) 


For  an  air  parcel  that  is  expanding  adiabatically,  the  first  law  of  thermodynamics  can  be  used  to  give: 


dT(z) 

dp(z) 


(3.6.4) 


where  C  is  the  heat  capacity  at  constant  pressure  per  unit  mass.  Combining  Equations  3.6.1  and  3.6.4 
gives  the  adiabatic  lapse  rate: 


adiabatic 


dT_^ 

V  dz  J  adiabatic 


(3.6.5) 


where: 


c 


p 


(3.6.6) 


and  ttl  and  /II  iV  are  the  mass  of  water  vapor  and  air  and  c  and  c  are  the  constant  pressure 

heat  capacities  of  water  and  air,  respectively.  The  adiabatic  lapse  rate  is  constant  for  a  dry  atmosphere  and 
is  given  by: 


dry ,  adiabatic 


—  = - — —  =  -0.00976  O' 


dz 


dry ,  adiabatic 


m 


(3.6.7) 


3.6.2  Atmospheric  Stability 

The  SCM  enables  the  user  to  specify  the  atmospheric  stability  based  on  the  Pasquill  stability  classes  given 
in  Table  8.  The  approach  is  useful  because  the  stability  classes  are  a  convenient  way  to  define  atmospheric 
stability  based  on  observation.  Once  specified,  the  stability  parameter  than  influences  the  calculation  of 
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other  atmospheric  properties  such  as  the  Friction  velocity  and  Monin-Obukhov  parameters  discussed  in 
Section  2.4.4. 

Table  8.  Pasquill  Atmospheric  Stability  Classes 


Stability 

Parameter 

Description 

Pasquill  Stability 

Class 

1 

Extremely  unstable 

A 

2 

Moderately  unstable 

B 

3 

Slightly  unstable 

C 

4 

Neutral 

D 

5 

Slightly  stable 

E 

6 

Moderately  stable 

F 

3.6.3  Sensible  Heat  Flux 

The  sensible  heat  flux  for  daytime  conditions  is  given  by: 


H  = 


1  -a  + 


(' Qrad  Q ground  )  Pf 


(3.6.8) 


where  Qrad  and  Qground  are  the  radiative  heat  flux  and  ground  heat  flux,  respectively.  The  parameter  a  is  a 
surface  moisture  coefficient.  The  parameters  y  and  s  are  given  by: 


r  =  c,x 

s  =  q/T 


(3.6.9) 


where  /,  is  the  latent  heat  of  vaporization,  Cp  is  the  constant  pressure  heat  capacity,  and  q  is  the  saturation 
specific  humidity. 

=  1 19.56  -  0.7843J  + 1. 2887x1  (T3J2  (3.6.10) 


The  radiation  heat  flux  is  given  by: 


Q 


rad 


(l  -  r)g„  +  CtT‘-oT‘  +  C2NrhlJ 

1  +  c3 


(3.6.11) 


where  T  is  the  temperature,  r  is  the  surface  albedo,  a  is  the  Stefan-Boltzman  constant,  and  Rsoi  is  the  solar 
radiation  at  ground  level.  The  empirical  constants  Q  and  C2  are  5.31x10-13  W  m"2  K6  and  60  W  m"2, 
repsectively.  The  constant  C3  is  a  coefficient  for  surface  and  depends  on  whether  the  surface  is  dry  or  wet. 
For  a  dry  surface  C3  =  0.25  and  for  a  wet  surface  C3  =  0. 12. 

The  solar  radiation  at  ground  level  in  the  absence  of  cloud  cover  can  be  expressed  in  terms  of  the  solar 
elevation  angle  <|>  by: 

Rso,=Rso,=a  2  +alsin</)  (3.6.12) 
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where  the  turbidity  coefficients  ai  and  a2  are  990  W  m"2  and  -30  W  m'2,  respectively.  When  cloud  cover  is 
present  then  the  radiation  at  ground  level  is  reduced  from  that  given  in  Eq.  2.4.12  and  can  be  expressed  as: 


(3.6.13) 


In  Equation  (3.6.13),  Tcioud  is  the  cloud  transmittance,  Ncioud  is  the  fraction  of  cloud  cover,  and  the 
coefficient  b2  =  3.4.  Tcioud  varies  by  the  type  of  cloud.  A  simple  set  of  data  is  Tcioud  =  0.61  for  high  altitude 
cirrus  type  clouds,  Tcioud  =  0.26  for  medium  altitude  cumulus  type  clouds,  and  Tdoud  =  0.17  for  low  altitude 
stratus  or  nimbostratus  clouds. 


3.6.4  Friction  Velocity  and  Monin-Obukhov  Length 

The  friction  velocity,  u*,  and  the  Monin-Obukhov  length,  LMO,  is  are  related  by: 
*  3 

PC pT  u 

L MO  =  ~ 

KgHsens 


(3.6.14) 


where  T  is  the  ambient  temperature,  p  is  the  ambient  density,  Cp  is  the  specific  heat,  g  is  the  gravitational 
constant,  k  is  the  von  Karman  constant,  and  H  is  the  sensible  heat  flux.  Additionally,  the  following 
relations  define  the  friction  velocity  for  neutral,  stable,  and  unstable  atmospheric  conditions  [Seinfeld  and 
Pandis,  1998]. 
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(3.6.15) 


where  u  is  the  measured  wind  speed,  z  is  the  wind  measurement  height,  z0  is  the  roughness  scale.  The 
parameter  a  is  given  by  a  =  z/L  when  z  <  L  and  by  a  =  1  when  z  >  L. 


3.7  Chemical  Kinetics 


3.7.1  Phenomenology 

Test  data  indicate  that,  after  detonating  or  burning  energetic  materials,  chemical  reactions  continue  as  the 
high  temperature  puff  cools  with  air  entrainment.  These  afterburning  reactions  influence  the  eventual 
chemical  composition  in  the  puff  and  also  release  heat  that  enhances  the  puff  buoyancy.  The  latter 
increases  the  final  puff  rise.  The  afterburning  is  a  coupled  process  between  chemistry  (thermochemical 
transformation)  and  physics  (turbulent  entrainment/mixing)  that  can  be  modeled  by  an  operator  splitting 
algorithm. 
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Semi-volatile  organic  compounds,  polychlorinated  dibenzo-p-dioxins  and  furans,  toxic  metals,  volatile 
organic  compounds,  and  inorganic  gases  have  been  measured  in  Bang  Box  tests  [Mitchell  and  Suggs,  1998] 
and  contained  detonation  and  exploding  ordnance  tests  [METDC,  2000;  2001].  Existing  models  applied 
thermodynamic  equilibrium  codes  throughout  the  entrainment  process  to  calculate  the  emission  factors. 
VOCs  were  not  predicted.  The  modeling  of  these  products  requires  the  application  of  a  chemical  kinetics 
solver  with  relevant  rate  parameter  data  retrieved  from  literature,  estimated  based  on  available  information, 
or  measured  in  laboratories. 

The  SCM  includes  a  module  to  treat  gas  phase  after-burning  kinetics  and  the  coupling  to  the  puff  flow 
dynamics  by  changes  in  the  enthalpy  due  to  chemical  reactions.  This  is  a  unique  aspect  of  the  current 
model  that  has  been  largely  lacking  in  other  models  for  explosive  events. 

3.7.2  Chemical  Kinetic  Mechanisms 

For  HE  munitions  detonation,  the  chemical  kinetics  of  particular  interest  include  the  gas  phase  afterburning 
kinetics  of  small,  unburned  hydrocarbons,  nitrogen  and  sulfur  oxides  and  acids,  and  volatile  (VOC)  and 
semi-volatile  (SVOC)  organic  compounds.  In  addition,  the  emission  of  carbonaceous  and  metal 
particulates  introduces  the  need  to  include  the  heterogeneous  uptake  kinetics.  Lastly,  decomposition  and 
oxidation  reactions  for  energetic  material  such  as  RDX  and  HMX  have  also  been  included  in  the  database 
due  to  the  detection  of  residual  energetic  material  at  several  DoD  sites  and  the  current  uncertainty  regarding 
their  source. 

Chemical  kinetic  mechanisms  were  developed  by  first  formulating  the  reactions  needed  to  describe  the 
chemistry  of  specific  species  or  sets  of  specific  species  and  then  combining  these  smaller  reaction  subsets 
to  formulate  the  kinetic  mechanisms  for  the  more  complex  systems.  The  reaction  subsets  incorporated  are 
shown  in  Figure  29. 

For  the  evolution  of  CO,  NOx,  and  relatively  small  VOCs  rate  parameters  were  adopted  based  on  evaluated 
kinetic  recommendations  from  the  combustion  and  atmospheric  modeling  communities  [Tsang  and  Herron, 
1991;  Baulch  et  al.,  1994;  DeMore  et  ah,  1997;  Atkinson  et  ah,  1997;  Sander  et  ah,  2003]  and  the  electronic 
NIST  chemical  kinetic  database  [Milliard,  1994].  For  more  complex  VOCs,  the  monographs  by  Atkinson 
[1989;  1994],  which  have  compiled  key  experimental  reaction  rate  constants  and  mechanisms  for  the 
oxidation  of  a  wide  range  of  VOCs  and  also  present  methods  for  estimating  key  oxidation  rate  parameters 
when  experimental  data  are  lacking  were  used.  An  illustrative  example  of  the  kinetic  mechanism  for  COx, 
HOx,  NOy,  NHX,  and  small  hydrocarbons  and  is  given  in  Section  7. 
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Figure  29.  Overview  of  SCM  chemical  kinetic  reaction  sets. 

3.7.3  Numerical  Model 

Through  previous  and  ongoing  projects,  Aerodyne  has  developed  a  chemical  kinetics  driver,  FLOWKIN, 
based  on  the  CHEMKIN-I1  [Kee  et  al.,  1989]  library  for  treating  finite -rate  chemical  kinetics  in  simple 
flows.  This  driver  is  a  modified  version  of  SENK1N  [Lutz  et  al.,  1991]  based  kinetic  integration  routine  for 
gas  phase  kinetics.  FLOWKIN  has  also  been  interfaced  with  AR1  computational  codes  for  treating  aerosol 
microphysics  [Brown  et  al,  1996a]  enabling  a  unified  treatment  of  the  chemistry  for  both  vapor  and 
condensed  phases.  The  entire  chemical  package  has  been  integrated  with  the  SCM  module  for  buoyant 
puff  rise.  This  provides  the  numerical  tools  required  to  treat  after  burning  chemistry  in  the  buoyant  puff,  as 
well  as  include  a  source  term  for  enthalpy  changes  due  to  chemistry  that  can  impact  the  puffs  dynamics. 

3.7.4  RDX  Half-life 

As  an  example,  the  SCM  chemical  kinetic  model  was  used  to  characterize  the  half-life  of  gas  phase  RDX 
during  the  initial  puff  expansion.  The  kinetic  mechanism  is  given  in  Appendix  5.  Figure  30  shows  selected 
decomposition  and  oxidation  products.  Figure  31  shows  the  RDX  half-life.  In  general,  the  half-life  is  short 
indicating  very  little,  if  any,  gas  phase  RDX  is  expected  to  survive  the  initial  adiabatic  expansion.  These 
gas  phase  kinetic  calculations,  however,  do  not  address  the  extent  to  which  RDX  is  released  as  a 
component  of  particulates  from  low  order  detonations  or  quenched  out  of  the  fireball  by  entrained  soil 
particles. 
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RDX  Decomposition  and  Oxidation  in  Post-Detonation  Gases 
(P0  =  100  atm,  T0  =  930.6  K,  V0  =  24.2  cc  g'1) 


RDX 
-©-  H2CN 

h2cnno2 

--V--  HCN 
N02 
N20 
—A—  HNO 
-A-  HONO 


time  (ps) 


Figure  30.  RDX  decomposition  and  oxidation  in  post-detonation  gases 


RDX  Half-Life  in  Post-Detonation  Gases 
as  a  Function  of  Temperature  During  Expansion 
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Figure  31.  RDX  half-life  and  C-J  adiabat  pressure  as  a  function  of  puff  temperature. 
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3.8  Heterogeneous  Update  Kinetics 


3.8.1  Overview 

The  reacting  detonation  puff  will,  in  general,  contain  two  types  of  particles,  entrained  soil/dust  and 
detonation  produced  soot,  that  may  interact  chemically  with  the  trace  gaseous  pollutants  produced  in  the 
detonation  of  energetic  materials  and/or  in  their  subsequent  afterburning  reactions. 

Soil  and  dust  particles  that  are  entrained  into  the  detonation  fire  ball  will  be  processed  at  very  high 
temperature,  with  organic  soil  components  largely  oxidized  or  pyrolyzed  to  soot.  The  remaining  particles 
will  be  primarily  inorganic  oxide,  hydroxide  and/or  carbonate  minerals.  Typical  soil  minerals  include: 

Si02,  AI2O3,  Fe203,  CaO,  CaC03,  MgO,  Ti02  and  various  clays  (typically  mixtures  of  Ca,  Mg,  Al,  Fe,  Na 
and  K  oxides,  hydroxides,  and/or  carbonates.) 

Soot  is  a  partially  hydrogenated  graphitic  carbon  material  with  a  nominal  C/H  ratio  of  ~1.  It  may  also 
contain  minor  amounts  of  oxygen,  nitrogen,  sulfur,  and  metallic  atoms,  depending  on  the  composition  of 
the  energetic  material  (fuel)  and  the  availability  of  oxidizer  (oxygen,  nitrogen  oxides)  during  the  initial 
detonation  process.  Soot  may  also  be  coated  with  a  variety  of  aromatic  hydrocarbon  compounds,  including 
polycylic  aromatic  hydrocarbon  (PAH)  compounds,  especially  when  it  is  generated  under  fuel  rich 
conditions  (Sander  et  al.,  2003.) 

Gaseous  species  can  be  physisorbed  or  chemisorbed  on  airborne  particles;  they  can  also,  in  some  cases, 
react  directly  with  the  particle’  surface  material  (first  order  reaction)  or  an  adsorbed  species  can  react  with 
another  gaseous  species  arriving  at  the  particle’s  surface  (second  order  reaction).  The  first  three  processes 
can  be  characterized  by  a  simple  uptake  parameter  (y)  that  is  defined  as  the  probability  that  a  gaseous 
species  striking  the  particle’s  surface  will  either  stay  adsorbed  on  the  surface  or  be  lost  via  reaction;  with  a 
range  of  0<y<  1  (Sander  et  al.,  2003;  Usher  et  al.,  2003a).  Under  some  circumstances  even  a  bimolecular 
surface  process  can  be  characterized  by  a  single  uptake  coefficient  if  one  reactant  is  much  more  abundant 
than  the  other,  its  uptake  may  be  controlled  by  the  y  value  for  the  less  abundant  species.  In  general,  the 
heterogeneous  uptake  processes  relevant  for  the  short  duration  of  the  reactive  detonation  puff  will  be 
physisorption  ,  chemisorption,  or  first  order  reactions,  so  specification  of  an  uptake  coefficient,  y,  will  be 
sufficient  to  describe  the  heterogeneous  kinetics. 

3.8.2  Soil  Particle  Uptake  Kinetics 

Heterogenous  uptake  kinetics  on  soil  dust  minerals  and  mineral  compositions  has  recently  been  reviewed 
by  Usher  et  al.  (2003).  While  relatively  few  laboratory  studies  aimed  at  providing  mineral  aerosol  and/or 
surface  uptake  parameters  have  been  published,  data  are  available  for  selected  nitrogen  oxides,  ozone, 
sulfur  dioxide,  and  selected  partially  oxygenated  organic  compounds  on  a  range  of  individual  mineral 
compounds  as  well  as  selected  types  of  soil  dust  particles.  Most  studies  are  based  on  either  Knudsen  cell  or 
fast  flow  reactor  techniques.  Many  studies  provide  time  resolved  uptakes  that  show  relatively  fast  uptake 
until  surface  saturation,  followed  by  slower  uptake  kinetics.  Since  soil  particles  entrained  through  the 
detonation  fireball  will  tend  to  be  “baked”  free  of  adsorbed  surface  species,  the  initial  uptake  coefficient,  y°, 
provided  by  the  time  dependent  experiments  is  deemed  to  be  most  representative  for  detonation  puff 
modeling.  Data  for  individual  reactive  detonation  species  is  presented  below.  In  general,  uptake  on  very 
unreactive  surfaces  like  Si02  is  much  less  efficient  than  uptake  on  more  labile  surfaces  like  CaO,  CaC03, 
MgO,  or  A1203.  Since  the  SCM  model  will,  in  general,  not  have  a  detailed  soil  chemical  analysis  to  input, 
a  single  uptake  coefficient,  y°,  representing  a  weighted  average  of  available  data  is  recommended  for  each 
trace  gaseous  species.  Where  data  are  available  for  actual  soil  dust  samples  the  uptake  coefficients  for 
those  samples  have  been  weighted  heavily.  Nearly  all  studies  where  performed  at  room  temperature  (295- 
298  K),  so  the  results  are  most  appropriate  for  detonation  puff  conditions  after  the  initial  expanding  fireball 
is  cooled  and  diluted  by  entrained  air.  For  warmer  puffs  the  room  temperature  values  will  likely  be  too 
high  for  physisorption  and  chemisorption  processes,  but  may  underestimate  the  efficiency  of  unimolecular 
surface  reactions.  Since  air  entrainment  and  puff  cooling  is  rapid  the  recommended  y°  values  should  be 
reasonably  accurate  for  most  of  the  time  scale  modeled  by  the  SCM. 
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NO 


Mineral  surface  uptake  studies  for  nitric  oxide,  NO,  have  been  reported  by  Hanish  and  Crowley  (2003a)  on 
Saharan  dust  and  Judeikis  and  Wren  (1978)  sandy  loam,  adobe  clay  and  cement.  The  Hanisch  and  Crowley 
study  demonstrated  that  while  surface  reaction  for  NO  did  take  place  with  the  co-deposition  of  significant 
amounts  of  ozone,  no  measurable  uptake  took  place  in  the  absence  of  ozone.  From  their  data  a  y°  of  <1  x 
10"5  can  be  estimated.  The  Judeikis  and  Wren  study  measured  deposition  velocities  for  NO  and  NO2  rather 
than  uptake  coefficients,  their  data  are  consistent  with  uptake  coefficients  on  soil  for  NO  in  the  range  of  3 
to  8  less  than  those  for  N02.  Based  on  the  y°N02  discussed  below,  a  y°N0  of  5  x  10"6  is  recommended. 

N02 

Relevant  nitrogen  dioxide  uptake  kinetics  have  been  reported  by  Underwood  et  al.  (2001a)  and  Ullerstam  et 
al.  (2003).  Knudsen  cell  studies  were  performed  by  Underwood  et  al.  a-Al203,  y-Al203,  a-Fe203,  y-Fe203, 
Ti02,  MgO,  CaO,  Saharan  sand,  and  Chinese  loess.  Lowest  y°  values  were  recorded  for  y-Al203  (2.2  x  10" 

8)  and  Saharan  sand,  which  is  largely  Si02  (1  x  10'6).  Significantly  larger  values,  in  the  range  of  1.1  to  8.5  x 
10"5  where  measured  for  a-Al203,  a-Fe203,  MgO,  CaO  and  Chinese  loess,  while  Ti02  at  9  x  10"6  and  y- 
Fe203  at  4  x  10"6  had  intermediate  uptake  coefficients.  Ullerstam  et  al.  measured  a  small  value  of  6.3  x  10"7 
for  Cape  Verde  dust,  which  is  mostly  Si02.  Based  on  this  wide  range  of  values  an  estimated  weighted 
mean  y°N02  of  2  x  1 0"  is  recommended,  which  is  about  half  of  the  value  measured  for  China  loess  by 
Underwood,  et  al. 

HNO3 

Nitric  acid  vapor  uptake  studies  on  mineral  surfaces  have  recently  been  reported  by  Underwood  et  al. 
(2001a;  2001b)  and  by  Flanisch  and  Crowley  (2001a;  2001b).  While  their  actual  laboratory  results  are  in 
reasonable  accord,  their  data  analyses  lead  to  significantly  different  reported  values  for  y°hno3-  The  point  of 
controversy  is  whether  the  proper  surface  area  accessed  by  the  experiments  is  the  geometrical  top  surface 
area  of  the  mineral  sample  of  a  much  larger  surface  area  estimated  with  a  “pore  diffusion”  model.  This 
controversy  is  not  new;  it  has  surfaced  previously  in  uptake  studies  on  ice  and  alkali  salt  particles  (Sander 
et  al.,  2003).  In  general  for  small  molecules  with  relatively  low  uptake  coefficients,  using  the  top 
“geometrical”  surface  area  is  not  correct  and  corrections  for  pore  diffusion  are  required.  Flowever,  HN03 
vapor  tends  to  have  a  large  uptake  coefficient  on  many  surfaces  and  until  the  top  surface  is  saturated,  very 
few  of  its  vapor  molecules  may  make  it  past  the  top  surface. 

Based  on  a  pore  diffusion  analysis,  Underwood  et  al.  (2001a)  report  y°HN03  values  that  range  from  1.1  x  10"3 
to  1.6  x  10  2  for  on  a-Al203,  CaC03,  MgO,  and  China  loess  with  the  value  for  basic  CaO  exceeding  the 
others  by  roughly  and  order  of  magnitude.  Underwood  et  al.  (2001b)  also  report  results  for  several  other 
mineral  oxides,  Gobi  dust  and  Saharan  sand.  Values  ofy°HN03  reported  for  these  samples  range  from  3.7  x 
10"3  to  2.0  x  10"\  Both  Hanisch  and  Crowley  papers  analyze  their  results  using  geometric  surface  areas,  but 
demonstrated  that  y°HN03  values  measured  on  particle  of  significantly  different  sizes  were  consistent  using 
this  assumption.  They  report  results  for  a-Al203,  CaC03,  Saharan  dust,  Arizona  dust,  Chinese  dust,  and 
eight  clays  that  range  from  0.06  to  0.196.  For  nitric  acid  vapor  uptake,  the  analyses  of  Hanisch  and 
Crowley  appear  to  be  more  compelling  and  based  on  their  results  a  y°hno3  of  0.13  is  recommended. 

o3 

Ozone  uptake  coefficients  on  mineral  surfaces  have  been  reported  by  Hanisch  and  Crowley  (2003b), 

Michel  et  al.  (2003),  and  Usher  et  al.  (2003b).  Hanisch  and  Crowley  report  a  mean  value  of  3.5  x  10"4for 
Saharan  dust.  Michel  et  al.  measured  y°03  values  for  a-Al203,  a-Fe203,  Si02,  kaolite,  China  loess,  and 
Saharan  sand  ranging  from  2  x  10"4to  2.7  x  10"5;  while  Usher  et  al.  of  the  same  lab  showed  that  pre-treating 
alumina  with  nitric  acid  decreased  y°03  by  72%  while  treatment  with  sulfur  dioxide  increased  y°03  by  33%. 
They  also  treated  Si02  with  saturated  and  unsaturated  hydrocarbons,  causing  its  y°03  to  decrease  40%  and 
increase  40%,  respectively.  Based  on  these  data  a  mean  y°03  value  of  5  x  10  is  recommended. 


51 


so2 

Sulfur  dioxide  uptake  on  a  variety  of  relevant  surfaces  have  been  studied  by  Ullerstam  et  al.  (2002;  2003), 
Goodman  et  al.  (2001)  and  Usher  et  al.  (2002).  Ullerstram  et  al.  measured  y°so2  values  ranging  from  4.1  x 
10"7  to  4.6  x  10"6  for  various  Cape  Verde  dust  samples.  Usher  et  al.  and  Goodman  et  al.  measured  y0S02 
values  for  a-Al203,  CaC03,  a-Fe203,  MgO,  Ti02,  Si02,  and  China  loess  ranging  from  5.1  x  10"4  (MgO)  to 
<1  x  10'7  (Si02).  China  loess  yielded  a  y°S02  value  of  3  (±  1)  x  10"5,  matching  a  recommended  mean  value 
of  3  x  10"5. 

VOCs 

Li  et  al.  (2001)  and  Carlos-Cueller  et  al.  (2003)  measured  uptake  coefficients  for  a  variety  of  partially 
oxidized  volatile  organic  compounds  (VOCs)  for  a  range  of  representative  inorganic  mineral  oxides.  Based 
on  these  measurements  y°voc  values  for  formaldehyde  (8  x  10"5),  acetaldehyde  (4  x  1  O'5),  propionaldehyde 
(6  x  10"5),  acetone  (8  x  10"5),  acetic  acid  (2  x  10"3),  and  methanol  (1  x  10"4)  can  be  recommended.  The  y°s 
for  acetic  acid  and  methanol  can  be  used  for  other  light  organic  acids  and  alcohols,  respectively. 

3.8.3  Soot  Uptake  Kinetics 

Heterogeneous  uptake  kinetics  on  soot  have  been  recently  reviewed  by  Sander  et  al.  (2003).  Uptake  on 
soot  usually  shows  a  relatively  high  initial  uptake  coefficient  followed  by  a  drop  off  as  reactive  or 
absorptive  sites  saturate.  These  initial  uptake  coefficients  are  most  appropriate  for  fresh  detonation  soot 
processed  through  the  fireball. 

no2 

Initial  uptake  of  N02  for  a  variety  of  fresh  soots  where  measured  by  Stadler  and  Rossi  (2000)  and  Alcala- 
Jomod  et  al.  (2000)  using  a  novel  cylindrical  reactor  technique  operating  at  room  temperature.  Initial 
uptake  coefficients  on  hexane,  decane,  toluene,  and  acetylene  soot  where  all  measured  to  be  0. 1  or  larger. 
Based  on  these  measurements  a  y°N02  for  soot  of  0. 1  is  recommended. 

hno3 

The  uptake  of  nitric  acid  vapor  was  studied  by  Rogaski  et  al.  (1997)  and  Choi  and  Leu  (1998)  at  room 
temperature.  Rogaski  et  al.  measured  a  y°HN03  of  0.038  for  Degussa  FW2  soot,  while  Choi  and  Leu 
obtained  y°uN03  values  of  0.060,  0.023,  and  0.067  on  kerosene,  hexane,  and  Degussa  FW2  soots, 
respectively.  Based  on  these  results  a  mean  y°HN03  of  0.05  is  recommended 

o3 

The  uptake  of  ozone  on  soot  displays  a  high  initial  uptake  rate,  followed  by  a  rapid  drop  off  to  a  much 
lower  steady  state  rate  (Longfellow  et  al.,  2000;  Sanders  et  al.,  2003).  Knudsen  cell  studies  with  good  time 
resolution  have  shown  that  the  initial  y°o3  can  be  as  large  as  1  x  10~3  (Rogaski  et  al.,  1997,  Stephens  et  al., 
1986).  A  y°03  of  1  x  10'3  is  recommended. 


so2 

The  uptake  of  sulfur  dioxide  on  soot  surfaces  has  been  reported  by  Rogaski  et  al.  (1997  and  Koehler  et  al. 
(1999).  The  former  reported  a  y°so2  of  3  x  1 0"3,  while  the  latter  report  a  value  of  2  x  10°.  A  value  of  2.5  x 
10"3  is  recommended. 
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nh3 

Muenter  and  Koehler  (2000)  studied  the  uptake  of  gaseous  ammonia  on  n-hexane  soot  as  a  function  of 
temperature  and  saw  no  significant  uptake  at  temperatures  above  173  K.  A  y°nh3  of  0  is  recommended. 

VOCs 

No  studies  of  the  uptake  kinetics  of  gaseous  VOC  compounds  on  soot  have  been  identified.  However,  soot 
and  related  black  carbon  materials  are  known  to  have  a  significant  absorptive  capability  for  both  partially 
oxidized  and  aromatic  VOCs.  In  the  absence  of  experimental  data,  a  y°voc  of  1  x  10"2  is  recommended  for 
organic  acids,  aldehydes,  ketones,  alcohols  and  aromatic  compounds. 

3.8.4  Summary 

The  available  scientific  literature  published  through  2003  on  the  heterogeneous  kinetics  of  gaseous 
pollutants  with  soil  mineral  and  soot  particulates  has  been  reviewed.  Published  uptake  kinetics  data  has 
been  assessed  to  develop  recommended  initial  absorptive  and/or  reactive  uptake  coefficients,  y°s,  for 
detonation  pollutant  products  of  environmental  interest.  The  recommended  uptake  parameters  should  be 
updated  as  additional  laboratory  data  becomes  available. 

3.9  Energetic  Material  Model 

Measurements  have  demonstrated  the  presence  of  energetic  residue  due  to  the  detonations  of  HE  munitions. 
For  high  order  detonations,  little  un-reacted  energetic  material  is  expected  since  essentially  all  of  the 
energetic  material  is  consumed  in  the  detonation  and  any  energetics  entrained  in  the  fireball  will  undergo 
thermal  decomposition,  c.f.  Section  3.7.  Low-order  detonations  can  yield  significantly  higher  levels  of 
residual  energetic  materials  so  that  in  addition  to  knowledge  about  the  levels  of  released  energetic  material, 
estimates  for  the  probability  or  percentage  of  misfired  munitions  resulting  in  a  low-order  detonation  are 
required. 

Mechanisms  for  the  Ejection  and  Deposition 
of  Energetic  Material  in  HE  Detonations 

Ablation  from  Energetic  Material  Ejection  Local  Combustion  Extinction 

Unconfined  Charges  for  Confined  Charges  and  Scavenging 

Expanding  Plume 


•**  Energetic  Material 

Particle 

Ablation 

-  |  Explosive  Charge  |  * 

Shock  Wave 


Mechanisms  for  ejection  and  deposition  of  energetic  material  formulated 
based  on  field  measurements  and  data  published  in  archival  literature. 


Figure  32.  Graphical  Representation  of  Mechanisms  for  the  Ejection  and  Deposition  of  Energetic  Material 
in  HE  Detonations. 
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There  are  several  potential  mechanisms  for  the  release  of  small  amounts  of  residue  emissions  including: 
ablation,  ablation  and  local  combustion  extinction  on  casing  fragments,  and  combustion  extinction 
associated  with  particle  scavenging  of  entrained  material.  These  mechanisms  are  depicted  graphically  in 
Figure  32.  Surface  ablation  of  energetic  molecules  and  molecular  clusters  has  been  demonstrated  in 
laboratory  studies  of  shock  detonation  of  common  energetic  materials  [Blais  et  ah,  1997;  Engelke  et  ah, 
2001], 

However,  modeling  these  energetic  material  ejection  and/or  combustion  extinction  processes  reliably  is 
exceedingly  difficult.  Consequently,  the  SCM  treats  energetic  residue  emissions  using  a  data-driven  model 
based  on  the  experimental  data  from  recent  field  measurements  [Pennington,  2002;  Hewitt,  2003].  Many 
of  these  tests  utilize  the  systematic  approach  developed  by  Jenkins  and  co-workers[2000a,  2000b,  2001] 
which  relies  on  a  fresh  snow  surface  as  a  collection  medium.  Available  field  test  data  are  listed  in  Tables  9 
and  10  for  live-fire  and  blow-in-place  detonations  respectively. 


Table  9.  Average  RDX  and  TNT  mass  deposited  from  live-fire  detonations 


Munitions 

Detonation 

Main  Charge 

Average  Mass  Deposited  (pg) 

RDX 

TNT 

60-mm  mortar 

live-fire 

Comp  B 

74 

18 

40-mm  rifle 
grenade 

live-fire 

Comp  B 

1608 

5.2 

M67  hand  grenade 

live-fire 

Comp  B 

25 

ND 

81 -mm  mortar 

live-fire 

Comp  B 

8500 

1100 

120-mm  mortar 

live -fire 

Comp  B 

4200 

320 

105-mm  howitzer 

live-fire 

Comp  B 

95 

170 

Table  10.  Average  RDX  and  TNT  mass  deposited  from  blow-in-place  detonations 


Munitions 

Detonation 

Main  Charge 

Average  Mass  Deposited  (pg) 

RDX 

TNT 

C4  demolition 
charge 

blow-in-place 

C4 

15 

not  reported 

81 -mm  mortar/C4 

blow-in-place 

Comp  B 

14 

0.081 

M15  anti-tank/C4 

blow-in-place 

Comp  B 

40 

0.076 

M19  anti-tank/C4 

blow-in-place 

Comp  B 

2.7 

ND 

155-mm 

howitzer/C4 

blow-in-place 

TNT 

not  reported 

7889 

PMA- 1 A  mine 

blow-in-place 

TNT 

not  reported 

690 

PPM2  mine 

blow-in-place 

TNT 

1.2 

276 

PPM2  mine  w/C4 

blow-in-place 

TNT 

48.5 

4500 

VS50  mine  w/C4 

blow-in-place 

RDX 

135 

not  reported 

As  noted,  released  energetic  material  from  misfired  munitions  is  an  important  source  of  range 
contamination.  Table  10  above  shows  the  estimated  mean  mass  of  RDX  and  TNT  released  with  blow-in¬ 
place  detonations  from  recent  field  measurements.  These  data  were  generally  acquired  by  attaching  a  small 
charge  to  a  test  munition  and  blowing-it-in  place  and  are  not  from  misfires.  Further,  some  of  the  measured 
energetic  material  residues  may  be  due  to  the  detonation  charge,  and  not  the  test  article.  Consequently,  the 
results  should  be  taken  to  as  indicative,  not  completely  descriptive,  of  the  energetic  material  residues 
created  by  low  order  detonations. 

To  implement  the  energetic  material  model,  mean  mass  deposited  per  mass  of  explosive  for  given 
munitions  were  collected  from  reported  data  for  various  munitions.  The  deposited  mass  is  associated  by 
the  experiments  with  given  munitions  and  type  of  detonation,  e.g.  impact,  blow-in-place,  or  proximity. 
Additionally,  the  average  values  for  different  munitions  with  the  same  primary  explosive  charge  were 
adopted  to  provide  one  method  of  estimation  for  munitions  that  contain  the  same  primary  charge  but  whose 
energetic  residue  emissions  have  not  been  explicitly  studied. 
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In  practice  for  high-order  detonations,  the  reported  mass  deposited  is  scaled  by  the  tested  explosive  mass  to 
provide  mass  deposited/explosive  mass.  The  latter  is  then  multiplied  by  the  explosive  mass  for  the 
munitions  under  consideration  to  estimate  the  mass  of  energetic  residual  material  resulting  from  a 
detonation.  For  multiple  test  items,  the  explosive  mass  per  item  and  the  total  number  of  items  are  utilized 
to  estimate  the  deposited  mass  of  residual  energetic  material: 


deposited 


UNitemsM 


exp  losive,i 


reported ,i 

M  TestExplsovie,i 


(3.9.1) 


The  summation  in  Equation  (3.9.1)  is  over  different  types  of  munitions.  Recent  field  studies  clearly 
indicate  that  high-order  detonations  do  not  create  enough  energetic  material  residues  to  explain  the  degree 
of  contamination  found  on  military  firing  ranges.  Recent  studies  have  further  shown  that  low-order 
detonations  caused  by  misfired  munitions  create  many  orders  of  magnitude  higher  energetic  material 
residues. 

Even  though  only  a  few  percent  or  less  of  munitions  firings  will  be  subject  to  low  order  detonations,  much 
higher  residues  from  these  munitions  will  probably  dominate  the  rates  of  range  contamination.  Thus, 
Equation  (3.9.1)  must  be  modified  to  represent  both  the  lower  level  of  contamination  per  shot  from  high- 
order  detonations  and  the  higher  level  of  contamination  from  the  smaller  number  of  low  order  detonations: 


items, i  exp,/ 

(1-  p 

\  1  low-order,  if 

high, reported, i 

-P low-order  ,i 

low ,  reported  ,i 

TestExplsovie,i 

TestExplsovie,i 

i 

(3.9.2) 

In  order  to  use  Equation  (3.9.2)  to  predict  range  energetic  material  contamination  for  a  series  of  firings,  a 
knowledge  of  the  mean  residue  levels  both  of  the  mean  high  and  low-order  detonations  for  each  energetic 
material  in  munitions  type  is  required,  along  with  an  estimate  of  the  fraction  of  munitions  expected  to 
produce  a  low -order  detonation.  The  usefulness  of  Equation  (3.9.2)  relies  on  continued  efforts  to  develop 
data  of  the  sort  displayed  in  Tables  9  and  10  for  as  wide  range  of  munitions  as  possible  and  access  to  that 
data  through  available  databases  such  as  the  USAEC  database  discussed  below  in  Section  4.0. 

4.0  Army  Range  and  Munitions  Environmental  Data  (ARMED) 

ARMED  is  an  Oracle  database  system  developed  to  support  the  collection,  accessibility,  and  management 
of  emissions  data  gathered  from  the  munitions  testing.  ARMED  currently  stores  emissions  data  provided 
from  munitions  testing  conducted  at  Aberdeen  Test  Center  (ATC),  Aberdeen  Proving  Ground,  MD  and 
West  Desert  Test  Center,  Dugway  Proving  Ground  (DPG),  UT. 

The  ARMED  project  consists  of  two  phases.  Phase  I  consists  of  a  Microsoft  Access  script  used  to 
automatically  upload  data  from  Microsoft  Excel  spreadsheets  into  a  central  Oracle  database.  Munitions 
data  is  documented  on  a  standardized  Microsoft  Excel  template  by  ATC  and  DPG.  Personnel  at  the  U.S. 
Army  Environmental  Center  (USAEC)  Acquisition  and  Technology  Division  (AT)  maintain  the  template. 
Once  the  data  goes  through  an  initial  quality  assurance/quality  control  (QA/QC)  process,  it  is  uploaded  into 
the  Oracle  Database  by  way  of  the  Microsoft  Access  script.  Any  recognized  errors  are  identified  by  the 
automated  script  and  relayed  back  to  the  upload  administrator.  These  errors  are  then  fixed  and  the  upload 
process  begins  again  until  all  the  data  is  uploaded  successfully.  The  requirements  identified  for  Phase  I  of 
ARMED  were  identified  through  interviews,  requirements  gathering  use  case  analysis,  and  meetings  with 
functional  and  technical  leads.  There  is  no  set  practice  as  to  how  data  is  input  into  the  template.  The  test 
center  users  have  been  given  the  ability  to  add  new  information,  such  as  property  values,  weapon 
detonators,  etc.,  as  needed.  This  decision  was  made  to  expedite  the  data  collection  process  for  the  test 
center  so  that  the  user  would  be  able  to  add  information  on  an  as  needed  basis  without  first  having  to 
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contact  a  system  administrator.  It  is  expected  that  the  systems  analyst  will  make  changes  that  would  be 
needed  to  accommodate  test  center  user  additions  for  the  upload  tool  to  function  and  succeed. 


Phase  II  will  allow  ARMED  stakeholders  to  access  the  data  maintained  in  the  system  using  a  web  based 
application. 

4.1  Phase  I  Business  Practice 


The  current  ARMED  business  process  for  Phase  I  is  illustrated  below  in  workflow/ swim  lane  format. 


Figure  33.  ARMED  Phase  I  Business  Process 


The  ARMED  system  functionality  is  illustrated  below  in  context  diagram  format. 
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CONTEXT  DIAGRAM 
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Enter  Test  Results 
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Project  Office 


Return  Template 


Figure  34.  The  ARMED  system  functionality. 

4.2  User  Characteristics 


The  user  characteristics  for  this  system,  as  displayed  in  the  above  diagram,  shall  adhere  to  the  following: 

•  Test  Center  -  Users  shall  have  the  ability  to  read,  write,  and  submit  data  via  a  standardized 
template. 

•  Project  Officer  -  Users  shall  have  the  ability  to  read  and  write  data  in  the  standardized  template. 

•  Systems  Analyst  -  Users  shall  have  the  ability  to  read  and  write  data  in  the  standardized  template 
as  well  as  make  modifications  to  template  elements  and  drop-down  lists. 

•  System  Administrator  -  Users  shall  have  the  ability  to  read  and  write  data  in  the  standardized 
template  as  well  as  the  ability  to  modify  the  database  tables  to  accommodate  new  elements. 
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In  addition  to  the  specific  user  characteristics  listed  above,  there  will  be  a  stakeholder  end  user.  This  user 
will  be  able  to  view  the  data  stored  in  the  database  and  generate  pre-defined  reports  for  convenience  in 
analyzing  data. 


4.3  Data  Collection 

Phase  I  is  an  ongoing  process  to  successfully  upload  emissions  data  from  over  200  tested  munitions.  To 
date,  the  Army  has  successfully  tested  and  collected  data  for  over  140  different  types  of  munitions 
commonly  used  in  training  activities.  Munition  items  are  tactically  fired  in  a  closed  facility  to  capture  the 
airborne  emissions  generated.  The  emissions  data  is  collected  using  state  of  the  art  sampling  techniques. 

In  order  to  ensure  defensible  data  for  the  Army’s  future  use,  munition  emissions  are  collected  in  accordance 
with  U.S.  Environmental  Protection  Agency  (EPA)  approved  Detailed  Test  Plans.  During  each  test,  over 
280  chemical  species  and  particulates  are  monitored  to  determine  if  they  are  present,  and  in  what  quantities, 
when  weapons  are  functioned.  All  of  this  information  will  be  accessible  to  users/stakeholders  during  Phase 
II  of  the  project,  which  is  scheduled  for  testing  during  the  spring  of  2004. 

4.4  Data  Searching  and  Viewing 

Upon  logging  into  ARMED,  users  will  be  able  to  search  for  emission  data  using  four  filters.  The  user  can 
search  for  data  based  on: 

1)  the  Department  of  Defense  Identification  Code  (DoDIC); 

2)  the  weapon; 

3)  a  specific  chemical  compound;  and/or 

4)  the  nomenclature. 

Also,  the  user  has  the  ability  to  run  pre-defined  reports  and  export  data  to  a  Microsoft  Excel  spreadsheet. 
ARMED  currently  contains  three  pre-defined  reports,  from  which  the  user  can  apply  runtime  parameters. 
These  parameters  are  listed  in  the  footer  of  the  report.  The  reports  are  listed  below: 

•  Emissions  Factors/Compound  Report  -  This  report  lists  the  emission  factors  for  the  user- 
selected  compound.  The  report  shows  the  DODIC,  nomenclature  and  emission  factors 
as  the  Average  Adjusted  Emission  Factor,  Ib/item  and  Average  Adjusted  Emission 
Factor,  Ib/lb  Net  Explosive  Weight  (NEW).  Information  is  derived  from  data  stored  in  the 
database. 

•  Sampling  and  Analysis  Methods  Report  -  This  report  lists  the  sampling  and  analysis 
methods  used  to  obtain  emission  factor  data  for  the  user-selected  test.  These  methods 
have  been  submitted  to  the  U.S.  Environmental  Protection  Agency  (EPA)  prior  to  each 
test  for  approval.  A  more  detailed  description  for  each  of  these  methods  is  available  upon 
request  from  USAEC. 

•  Parameter/Compound  List  Report  -  This  report  lists  all  of  the  compounds  that  each 
munition  has  been  tested  for  as  potential  emission  factors  to  date.  The  report  groups  all 
of  the  compounds  by  parameter  type. 

4.5  Conclusion 

The  data  contained  in  ARMED  has  undergone  an  extensive  QA/QC  process  to  ensure  complete  data 
integrity.  USAEC  anticipates  going  live  with  Phase  II  of  the  database  for  the  summer  of  2004,  so  that  the 
emission  factor  information  will  be  accessible  to  users  and  stakeholders  via  the  web.  The  munition  items 
that  have  not  been  tested  at  the  Test  Centers  will  have  relevant  data  accessible  using  the  Source 
Characterization  Model. 
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5.0  Model  Validation 


5.1  Overview  of  Measurement  Data 

5.1.1  Nevada  Test  Site  Emissions  Characterization 

This  dataset  set  include  results  from  four  experiments  that  were  performed  at  the  Department  of  Energy's 
Nevada  Test  Site  (NTS)  [Watkins  et  al.  1998;  Velsko  et  al.  1999].  Tests  at  the  NTS  were  performed  for  the 
M107,  high  explosive  (Composition  B)  155-mm  projectiles.  Four  detonations  were  conducted  between 
December  1996  and  March  1997.  The  four  tests  included  the  simultaneous  detonation  of  6,  24,  24,  and  60 
projectiles.  An  additional  3  detonations  were  performed  between  July  and  August  1999.  Each  of  these 
tests  was  for  the  detonation  of  24  projectiles.  Two  of  these  latter  tests  consisted  of  the  24  projectile 
detonations  in  1997  to  determine  experiment  reproducibility.  The  third  test  changed  the  orientation  of  the 
projectiles  from  the  horizontal  to  the  vertical.  Measurements  inside  the  chamber  included  pressures, 
temperatures,  relative  humidity,  and  gas  concentrations.  In  addition,  grab  samples  were  collected  for  gas, 
organic,  metal  and  particulate  analyses.  A  general  description  of  the  results  of  these  measurements  has 
been  presented  by  Velsko  et  al.  [1999]  and  is  summarized  below. 

•  CO2  and  NO  showed  peaks  shortly  after  detonation  and  slowly  decreased  in  concentration 
thereafter  while  CO  changes  concentration  in  the  chamber  due  to  detonation  but  then  remains  at  a 
constant  volume  fraction  of  the  chamber  gases.  The  C02  concentrations  indicated  incomplete 
combustion.  For  the  four  detonations  in  1997,  the  overall  conversion  of  the  explosive  carbon  to 
gas  products  ranged  from  85-94%. 

•  For  volatile  organic  compounds  (VOCs),  the  emissions  of  EPA  regulated  materials  were  found  to 
be  present  in  trace  levels  of  <  10  ppmv. 

•  The  total  suspended  particulate  (TSP)  levels  were  in  the  range  of  hundreds  of  mg/m3.  The  size 
fraction  <10  microns  (PMi0)  was  in  the  range  of  tens  of  mg/m3. 

•  The  only  energetic  material  detected  in  these  tests  was  4 1  mg  RDX  /  kg  residue  for  the  test  in 
1997  using  6  projectiles. 

5.1.2  Aberdeen  Test  Center 

USAEC  funded  the  U.S.  Army  Aberdeen  Test  Center  (ATC),  Aberdeen  Proving  Ground  (APG)  Maryland, 
to  develop  the  Exploding  Ordnance  Emission  Study,  Phase  II  [METDC,  2000].  The  goal  of  this  study  is  to 
identify  and  quantify  the  emissions  resulting  from  detonation  of  ordnance  currently  fielded  by  the  U.S. 
Army.  The  Series  I  tests  have  been  completed  and  the  results  have  been  summarized  in  a  draft  report 
[METDC,  2001].  The  test  items  included  (1)  a  1/4-pound  demolition  charge  consisting  of  trinitrotoluene 
(TNT),  (2)  a  demolition  block  of  military  dynamite  (M591),  (3)  a  block  demolition  charge  consisting  of 
Composition  C4  (Ml  12),  (4)  a  fragmentation  hand  grenade  (M67),  and  (5)  a  66-mm  antitank  rocket 
(M72A2).  An  electric  blasting  cap  (M6)  is  used  to  initiate  the  high  explosives  although  the  blasting  cap 
itself  is  not  a  test  item. 

The  list  of  emission  measurement  target  analytes  for  the  test  ordnances  in  Table  1  is  extensive.  In  addition 
to  the  gases  carbon  monoxide  (CO),  carbon  dioxide  (CO?),  sulfur  dioxide  (S02),  ammonia  (NH3),  and 
nitrogen  oxides  (NO,  N02),  it  includes  inorganic  acids,  hydrocarbons,  aldehydes,  polyaromatic 
hydrocarbons  (PAHs),  volatile  organic  compounds  (VOCs),  semivolatile  organic  compounds  (SVOCs), 
dioxins/furans,  energetic  materials,  and  metals.  A  complete  list  of  the  emission  measurement  targets  is 
given  in  Table  2.3.3  of  the  test  plan  [METDC,  2000]. 
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5.1.3  Dugway  Proving  Ground 

USAEC  sponsored  a  6  phase  series  of  tests  utilizing  the  BangBox  at  DPG  to  quantify  the  emission  products 
produced  when  training  munitions  are  used  in  their  tactical  configurations  [Radian  International,  2003]. 
Phase  I  characterized  airborne  emissions  for  eight  simulator  and  pyrotechnic  munitions.  Phase  II 
considered  6  additional  ordnance  items.  Phase  III  studied  smoke  grenades.  Phase  IV  characterized  the 
airborne  emissions  of  4  exploding  ordnance  items. 

The  analysis  of  these  tests  using  the  SCM  treated  the  munitions  studied  in  Phase  IV  and  included  the 
following: 


M031 

Demolition  block 

1/2  lb  TNT 

M456 

Detonating  cord 

PETN  type  1 

ML09 

Flexible  Linear  Shaped  Charge 

20  g/ft 

ML  15 

Flexible  linear  Shaped  Charge 

225  g/ft 

In  summary,  CO  and  C02  were  the  primary  gas  phase  compounds  emitted.  TSP  was  the  primary  regulated 
compound  emitted  from  all  items.  The  highest  emissions  were  from  the  linear  shaped  charges.  Higher 
levels  of  RDX,  HCN,  NH3,  formaldehyde,  and  CO  indicated  less  than  complete  combustion  for  the  linear 
shaped  charges.  Airborne  metals  were  also  emitted.  There  was  no  indication  of  dioxins/furans  for  the  TNT 
demolition  block  or  the  detonation  cord.  Low  levels  of  dioxins/furans  were  reported  for  some  samples 
from  the  linear  shaped  charges.  Most  of  the  organic  compounds  measured  where  light  hydrocarbons 
(ethylene  and  acetylene).  Most  of  the  organic  regulated  compounds  were  VOCs,  reported  at  lto  2  orders  of 
magnitude  higher  than  the  SVOC. 

5.2  Aberdeen  Test  Center  Exploding  Ordnance  Emissions  Study  (EOES)-  Series  I,  Phase  II 
5.2.1  EOES  Test  Data 

Table  9  summarizes  the  test  items.  The  NTS  experiments  treated  the  M107  155-mm  rocket  with 
Composition  B  and  included  the  detonation  of  6,  24,  24,  and  60  projectiles.  The  USAEC  test  items 
included  (1)  a  1/4-pound  demolition  charge  consisting  of  trinitrotoluene  (TNT),  (2)  a  demolition  block  of 
military  dynamite  (M591),  (3)  a  block  demolition  charge  consisting  of  Composition  C4  (Ml  12),  (4)  a 
fragmentation  hand  grenade  (M67),  and  (5)  a  66-mm  antitank  rocket  (M72A3).  An  electric  blasting  cap 
(M6)  was  used  to  initiate  the  high  explosives  although  the  blasting  cap  itself  is  not  a  test  item. 
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Table  9:  Summary  of  USAEC/ATC  Exploding  Ordnance  Emission  Studies  -  Series  1  Test  Items 


USAEC 

NTS 

Ordnance 

M030 

M591 

Ml  12 

M72A2 

M67a 

M107 

Demolition 

Demolition 

Demolition 

Type 

charge 

charge 

charge 

66-mm  rocket 

Grenade 

155-mm  rocket 

PEP 

TNT 

Dynamite 

Composition  C4 

Octal 

Composition  B 

Composition  B 

PEP  weight  (lbs) 

0.25 

0.37 

1.25 

0.67 

0.41 

15.4 

Rounds 

1 

1 

1 

1 

1 

6,  24,  and  60 

PEP  composition 

TNT 

100.00% 

15.00% 

0.00% 

30.00% 

39.00% 

36.00% 

RDX 

0.00% 

73.00% 

91.00% 

0.00% 

60.00% 

63.00% 

HMX 

0.00% 

0.00% 

0.00% 

70.00% 

0.00% 

0.00% 

polyisobutylene 

0.15% 

2.25% 

mineral  oil 

1.68% 

1.31% 

Lubrizol 

0.16% 

0.13% 

dioctyl  sebacate 

5.31% 

guar  gum 

5.00% 

other 

0.01% 

Detonation 

M6  blasting  cap 

M6  blasting  cap 

M6  blasting  cap 

M6  blasting  cap 

M2 13  fuse 

RP1 

a.  Detonator  assembly  constituents  not  listed 


The  test  for  the  M67  grenade  and  M107  155-mm  rocket  are  of  particular  interest  because  Composition  B 
was  the  explosive  in  both  cases.  Consequently,  two  data  sets  for  Composition  B  are  available  for 
evaluating  the  results  of  the  Source  Characterization  Model  (SCM). 

Table  10  lists  a  portion  of  the  emission  factors  determined  under  the  Exploding  Ordnance  Emission  Study, 
Phase  II.  It  is  largely  limited  to  the  major  emissions  listed  in  Table  1.1-1  of  the  preliminary  data  report 
[METDC,  2001]  with  a  few  minor  exceptions.  The  emission  factors  are  reported  as  the  average  adjusted 
emission  factor  in  units  of  lb  of  analyte  per  lb  of  net  explosive  weight  (NEW).  As  noted  in  the  summary 
for  the  test  plan,  there  were  two  trials  with  one  detonation  per  trial  for  each  test  item.  The  emission  factors 
in  Table  2  are  the  average  for  the  two  trials  and  have  been  adjusted  for  background  levels  and  losses  due  to 
the  lack  of  containment.  The  water,  nitrogen,  and  hydrogen  emissions  were  not  directly  measured.  In 
addition,  no  data  was  reported  for  graphitic  carbon.  Also  notable  are  the  reported  levels  for  energetic 
material  including  RDX  for  all  test  items. 
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Table  10:  Average  Adjusted  Emission  Factor  from  Exploding  Ordnance  Emission  Study,  Phase  II 


M030 

M591 

Ml  12 

M67 

M72A3 

Permanent  Gases  and  Cyanide 

co2 

1.390 

0.464 

0.604 

0.517 

0.535 

CO 

2.03E-02 

7.96E-03 

2.00E-02 

4.07E-02 

7.41E-02 

NOx 

5.83E-03 

2.30E-03 

4.03E-03 

1.57E-03 

4.71E-03 

hno3 

1.46E-04 

1.25E-04 

1.06E-03 

6.12E-05 

ND 

HCN 

5.50E-04 

2.18E-04 

4.64E-04 

1.08E-03 

1.38E-04 

Particulates 

TSP 

6.92E-02 

2.33E-02 

2.80E-02 

9.65E-02 

9.43E-02 

Pio 

5.31E-02 

1.42E-02 

2.27E-02 

8.02E-02 

8.27E-02 

PM2.5 

1.98E-02 

5.49E-03 

1.67E-02 

4.55E-02 

4.77E-02 

Selected  Metals 

Aluminum 

1.14E-03 

2.26E-04 

1.86E-04 

7.11E-04 

1.67E-03 

Barium 

6.33E-04 

1.67E-05 

3.85E-06 

1.84E-04 

5.42E-04 

Chromium 

8.04E-04 

3.35E-04 

4.86E-06 

3.74E-04 

2.04E-04 

Copper 

2.31E-05 

3.71E-05 

1.91E-04 

2.95E-05 

2.24E-05 

Lead 

4.97E-04 

4.36E-05 

8.28E-06 

1.54E-04 

4.77E-03 

Magnesium 

5.93E-04 

2.39E-04 

1.04E-04 

1.30E-03 

1.85E-03 

Manganese 

6.94E-05 

5.34E-05 

1.52E-04 

6.25E-05 

8.40E-04 

Selected  Energetics 

2,4-dintrotoluene 

1.55E-06 

2.13E-07 

ND 

ND 

2.36E-07 

1,3,5-trinitrobenzene 

9.90E-6 

1.42E-06 

1.46E-06 

ND 

4.01E-07 

2,4,6-trinitrotoluene 

8.08E-05 

1.44E-05 

5.38E-07 

1.55E-06 

1.89E-06 

RDX 

2.41E-06 

2.15E-06 

1.58E-06 

2.31E-06 

2.99E-05 

Selected  Hydrocarbons 

methane 

6.14E-05 

3.35E-04 

1.25E-03 

6.80E-04 

8.64E-04 

ethylene 

2.31E-05 

9.84E-03 

1.04E-04 

5.56E-05 

6.85E-05 

acetylene 

3.63E-05 

9.84E-03 

4.27E-05 

1.87E-05 

9.58E-05 

Selected  Carbonyls 

formaldehyde 

7.25E-05 

3.12E-05 

1.21E-04 

8.04E-05 

6.36E-05 

acetaldehyde 

2.82E-05 

1.07E-05 

7.96E-06 

ND 

1.73E-05 

butyraldehyde 

1.21E-04 

ND 

5.91E-05 

1.16E-04 

3.36E-05 

ND  =  not  detected 


5.2.2  Mass  Balance 

A  mass  balance  analysis  was  performed  for  the  emissions  measured  under  the  EOES.  The  analysis  was 
limited  by  the  fact  that  experimental  data  for  water,  molecular  hydrogen,  molecular  nitrogen,  and 
particulate  carbon  were  not  reported.  Consequently,  for  the  present  analysis  empirical  carbon  mass  balance 
is  reported  based  on  the  gas  phase  emissions  that  were  measured  and  assuming  no  particulate  carbon, 
assuming  that  all  PM2.5  emissions  were  carbon  particulates,  and  assuming  all  total  suspended  particles 
(TSP)  were  carbon  particulates.  An  independent  hydrogen  mass  balance  was  not  performed  since  water 
vapor  and  molecular  hydrogen  emissions  were  not  reported.  Instead,  a  water  emission  index  was  calculated 
assuming  all  hydrogen  from  the  explosive  was  converted  to  water.  This  emission  factor  was  then  adjusted 
to  give  a  corrected  water  emission  index  by  taking  into  account  the  hydrogen  in  HCN  and  the  organics  for 
which  emission  levels  were  measured.  The  corrected  emission  index  for  water  was  then  used  in  the  oxygen 
mass  balance.  The  corrected  emission  index  for  water  is  expected  to  be  high  since  it  assumes  no  gas  phase 
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molecular  hydrogen.  An  independent  nitrogen  mass  balance  was  not  performed  because  no  data  was 
reported  for  molecular  nitrogen. 


Table  11.  Carbon  and  Oxygen  Mass  Balance  for  the  Exploding  Ordnance  Emission  Study,  Phase  II. 


M030* 

M591 

M112 

theor. 

empir. 

% 

theor. 

empir. 

% 

theor. 

empir. 

% 

carbon 

7 

7.34 

104.9 

2.84 

2.4 

85.0  % 

3.00 

3.7 

107.7 

% 

% 

with 

7.72 

110.3 

2.5 

88.2  % 

3.5 

118.0 

PM2.s 

% 

% 

with  TSP 

8.65 

123.6 

2.8 

98.4  % 

3.2 

125.0 

% 

% 

oxygen 

6 

17.06 

284.3 

4.64 

5.8 

124.2 

6.00 

9.3 

154.9 

% 

% 

M72 

M67 

Ml 07  (24  projectiles) 

theor. 

empir. 

% 

% 

% 

% 

% 

% 

% 

carbon 

5.1 

4.1 

79.5  % 

4.63 

2.9 

63.8  % 

4.67 

2.68 

62.9% 

with 

5.2 

100.8 

3.8 

82.0  % 

PM2.5 

% 

with  TSP 

6.1 

121.5 

4.7 

102.4 

% 

% 

oxygen 

7.4 

10.9 

148.2 

5.89 

8.4 

143.4 

5.86 

7.04 

120.1 

% 

% 

*  Excess  carbon  may  be  due  to  high  pre -detonation  chamber  CO2  from  a  large  number  of  observers  who 
inspected  the  chamber  prior  to  the  test  as  this  was  one  of  the  first  tests  in  the  series. 


Table  1 1  summarizes  the  results  of  the  mass  balance  analysis.  With  the  exception  of  TNT,  the  carbon  mass 
balance,  assuming  all  carbon  from  the  explosive  is  emitted  in  gas  phase  species,  ranges  from  about  80%  to 
100%  (a  little  over  100%  for  the  Ml  12).  This  suggests  that  depending  on  the  test  item  up  to  approximately 
20%  of  the  carbon  may  be  in  particulates.  The  block  demolition  charge  of  TNT  (M030)  is  the  most  oxygen 
deficient  of  the  explosives  tested  and  the  carbon  mass  balance  falls  well  below  100%  reaching  only  about 
50%  assuming  that  the  TSP  are  composed  of  carbon.  The  oxygen  mass  balance  ranges  from  1 12%  to 
155%.  These  numbers  partly  reflect  that  over  estimate  for  water  as  discussed.  For  the  most  part,  however, 
they  are  interpreted  to  reflect  afterburning  in  the  test  chamber. 

5.2.3  Detonation  Product  Emissions:  Kinetic  Results 

Post-detonation  chemical  kinetic  calculations  were  performed  to  compare  predicted  emission  levels  for 
selected  species  to  those  reported  as  part  of  the  EOES  data.  To  make  the  comparison  the  measured 
emission  factors  were  converted  to  mixing  ratios.  To  make  that  conversion,  as  well  as  to  do  the  kinetic 
calculations,  nominal  industrial  background  levels  for  trace  species  were  adopted:  C02  (400  ppmv),  CO  (4 
ppmv),  CH4  (2  ppmv). 

A  simple  treatment  of  the  expansion  and  entrainment  was  developed  in  which  the  radius  of  the  expanding 
puff  increases  proportionally  to  the  time  squared.  The  explosive  gas/air  mixing  ratio  was  taken  to  vary 
linearly  between  1  (all  explosive  gas)  to  fully  mixed  out  conditions.  Various  mixing  profiles  were  tested. 
These  results  are  for  a  profile  in  which  the  change  in  the  mass  fraction  of  explosive  varies  from  1  at  time  tQ 
fully  mixing  at  time  f  =  0.7  s  linearly.  In  addition,  a  few  parametric  calculations  were  performed  for 
chamber  volumes  25%  smaller  and  larger  then  reported.  All  these  parametric  studies  sometimes  resulted  in 
improvement  for  some  species,  but  resulted  in  a  poorer  comparison  to  measurement  for  others  than  the 
baseline  case. 
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However,  there  was  one  particular  important  aspect  of  the  calculation  that  played  an  important  part  in 
improving  the  correlation  between  SCM  predictions  and  measurements  for  NOx  emissions.  Specifically, 
the  SCM  predictions  for  NOx  levels  were  significantly  lower  than  measurements  for  each  of  the  five  test 
cases.  Analysis  of  the  reported  NOx  levels  related  to  the  amount  of  nitrogen  in  a  given  explosive  and  the 
NOx  levels  predicted  by  non-ideal  gas  equilibrium  at  detonation  suggested  that  the  measured  NOx  emission 
levels  might  be  indicative  of  an  alternate  kinetic  mechanism.  Specifically,  reported  emission  factors  for 
energetics  and  energetic  material  byproducts  suggest  that  some  energetic  material  is  not  consumed  during 
detonation.  Energetic  material  that  is  ejected,  while  not  fully  reduced  to  primary  detonation  products,  may 
still  experience  pressures  and  temperatures  that  induce  partial  decomposition  through  pyrolysis 
mechanisms.  The  pyrolysis  of  energetic  material  for  the  test  explosives  of  interest  here  would  result  in 
higher  NOx  levels  than  would  be  predicted  assuming  all  of  the  explosive  mass  is  consumed  in  the 
detonation. 

Consequently,  a  pyrolysis  model  for  enhanced  NOx  emission  was  developed  that  used  measured  levels  for 
residual  energetic  material  as  an  additional  NOx  source.  It  was  found  that  this  data-driven  pyrolysis  model 
had  no  impact  on  the  other  predictions  for  species  emission  factors  because  the  reported  levels  of  energetics 
were  relatively  small.  However,  it  did  increase  the  levels  of  NOx  in  the  isentropically  expanded  detonation 
puff  prior  to  further  afterburning  chemistry.  Afterburning  chemistry  was  treated  as  described  previously 
with  initial  conditions  derived  from  the  non-ideal  gas  equilibrium  with  chemical  kinetics  frozen  for 
temperatures  <  2000  K,  except  that  initial  NOx  levels  were  increased  based  on  the  proposed  pyrolysis 
mechanism  for  NOx  enhancement. 

Table  12  summarizes  the  SCM  model  predictions  for  each  of  the  test  cases  relative  to  reported  emission 
factors  from  measurement  by  listing  the  ratio  of  calculated-to-measured  emission  factors  for  several  key 
species.  The  predicted  to  measured  ratios  of  detonation  product  mixing  ratios  for  the  five  disparate 
munitions  listed  in  Table  12  are  very  encouraging.  Ratios  for  the  primary  detonation  products  (CO,  C02, 
CH20,  and  NOx  fall  within  an  order  of  magnitude  and  most  within  a  factor  of  three.  Predicted  mixing  ratio 
for  small  alcohols  such  as  ethanol  and  butanol,  as  well  as  5-6  carbon  alkanes  are  within  a  factor  of  3  to  10 
of  the  measured  values.  These  results  support  the  non-ideal  gas  equilibrium  freeze-out  of  carbon 
containing  species  presented  in  Section  3.1  as  modified  by  the  limited  afterburning  effects  as  discussed  in 
Section  3.7.  However,  calculated  mixing  ratios  for  benzene,  toluene,  napthalene,  xylenes,  and  alkenes 
such  as  1 -butene  and  isobutene  are  up  to  a  factor  of  100  times  smaller  than  measured  levels. 
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Table  12.  Ratio  of  SCM  Predicted  to  Measured  Mixing  Ratios  for  Selected  Detonation  Products 


Species\Test 

M030 

M591 

M112 

M72 

M67 

CO 

3.5 

2.0 

2.9 

0.9 

2.8 

C02 

0.4* 

0.9 

1.2 

1.1 

0.8 

ch2o 

1.3 

1.6 

0.3 

0.5 

0.2 

ch4 

6.8 

2.6 

2.3 

5.8 

0.3 

HCN 

3.5 

4.3 

0.4 

1.1 

9.5 

NOx 

0.8 

0.8 

0.9 

1.2 

0.9 

c2h4 

9.4 

3.1 

2.6 

4.2 

4.1 

c2h2 

1.5 

0.4 

0.3 

0.7 

0.2 

ch3hco 

0.2 

0.1 

NA 

NA 

2.0 

Propene 

4.2 

2.2 

2.4 

1.8 

2.5 

Ethanol 

0.4 

0.2 

0.3 

0.2 

0.2 

Butanol 

0.35 

0.3 

0.3 

0.1 

0.4 

Napthalene 

0.08 

0.05 

0.07 

0.05 

0.09 

Pentane 

0.4 

0.3 

0.5 

0.4 

0.6 

Hexane 

0.1 

0.1 

0.3 

0.2 

0.3 

1  butene 

0.08 

0.06 

0.07 

0.08 

0.06 

Isobutene 

0.08 

0.08 

0.05 

0.06 

0.08 

Benzene 

0.08 

0.06 

0.06 

0.07 

0.07 

Toluene 

0.01 

0.01 

0.02 

0.01 

0.03 

m/p/o-xylene 

0.01 

0.01 

0.03 

0.02 

0.02 

5.3  Dugway  Proving  Ground  Ordnance  Emissions  Characterization  -  Phase  IV 


In  June  2001  Phase  IV-A  testing  was  performed  using  the  BangBox  at  DPG  under  a  USAEC  program  to 
measure  the  emissions  from  weapons  firing,  pyrotechnic  devices,  and  exploding  ordnance  used  as  part  of 
range  operations.  URS  provided  air  sampling  and  analytical  support  in  this  program  and  the  data  given  in 
their  report  [Carson  et  al,  2003]  were  used  to  specify  SCM  initial  conditions  and  to  compare  calculated 
emission  factors. 

The  Phase  IV-A  test  items  included  a  demolition  block  charge  of 'A  lb  TNT,  a  PETN  detonating  cord,  and 
two  flexible,  linear  demolition  charges.  These  test  items  and  their  chemical  compositions  are  listed  in 
Table  13  along  with  some  other  general  data.  For  each  item,  the  composition  given  in  Table  13  is  from 
Carson  et  al.  [2003]  and  was  calculated  based  on  data  listed  in  the  MIDAS  database.  The  composition 
includes  the  test  item  and  energetic  portions  of  the  blasting  cap.  These  compositions  were  used  to  specify 
initial  conditions  for  the  SCM  calculations. 
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Table  13:  Summary  Items  Tested  at  PPG 


Test  Item  NEW  and  Composition 

M031 

M456 

ML09 

ML15 

Charge, 

Charge, 

Charge, 

Demolition 

Cord, 

Demolition, 

Demolition, 

Type 

Block 

Detonating 

Linear 

Linear 

Explosive 

TNT 

PETN 

RDX 

RDX 

NEW  (lbs) 

0.49 

0.49 

0.481 

0.482 

No.  items  /length 

1  item 

70  ft 

168  ft 

15  ft 

No.  trails 

2 

2 

2 

2 

Composition 

lb/lb  NEW 

C 

0.371 

0.190 

0.180 

0.180 

RDX 

4.32E-03 

4.23E-03 

0.982 

9.84 

Pb 

9.00E-04 

9.00E-04 

1.85E-03 

1.84E03 

HC1 

3.80E-05 

3.80E-05 

7.75E-05 

7.73E-05 

Emissions  measurements  were  performed  at  the  DPG  BangBox  test  facility  consisting  of  a  15.2m  diameter 
hemisphere.  Test  items  were  suspended  by  a  copper  wire  a  distance  of  1 .5m  above  the  chamber  floor. 
Tracer  gas  was  released  in  the  chamber  after  item  deployment  and  the  chamber  was  cleaned  between  tests. 
Table  14  lists  measured  emission  factors  for  selected  species.  The  measured  carbon  was  nearly  equal  to  the 
initial  level  of  carbon  based  on  the  compositions  in  Table  13  for  all  items  except  the  detonating  cord. 

The  DPG  tests  were  modeled  with  the  SCM  using  the  same  approached  used  to  model  the  ATC  tests 
described  in  Section  5.2.  Table  15  lists  the  ratio  of  calculated  mixing  ratio  to  the  average  value  calculated 
from  the  reported  [Carson,  2003]  average  corrected  emission  factor. 
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Table  14:  Representative  Emission  Factors  for  PPG  Tests  from  Carson  et  al.,  2003. 


Pollutant 

Average  Corrected  Emission  Factor  (Ib/lb  NEW) 

M031 

M456 

ML09 

ML15 

Particulates 

TSP 

3.43E-01 

3.80E-01 

1.61E+00 

5.12E+00 

PMio 

1.15E-01 

3.40E-01 

3.85E-01 

7.56E-01 

PM2.5 

1.31E-02 

1.1  IE-02 

1.36E-02 

8.09E-02 

Metals 

Aluminum 

1.14E-02 

1.63E-02 

5.06E-03 

7.70E-03 

Antimony 

6.54E-05 

7.56E-05 

9.37E-02 

3.00E-01 

Copper 

3.20E-03 

1.33E-03 

1.03E-03 

4.40E-03 

Lead 

2.20E-04 

7.26E-05 

3.49E-01 

8.50E-01 

Magnesium 

8.81E-03 

1.02E-02 

3.54E-03 

7.83E-03 

Zinc 

1.73E-03 

1.78E-03 

1.02E-03 

5.34E-03 

Energetics 

1,3,5-Trinitrobenzene 

1.79E-05 

4.89E-07 

3.02E-06 

2.52E-06 

RDX 

5.08E-07 

1.72E-06 

2.33E-02 

2.09E-02 

SVOCs 

Naphthalene 

7.61  E-06 

2.63E-05 

2.21  E-05 

1.79E-05 

Phenol 

7.73E-06 

5.50E-06 

2.47E-05 

2.64E-05 

VOCs 

Ethylene 

1.58E-04 

1.58E-03 

2.10E-03 

2.27E-03 

Acetylene 

1.03E-04 

1.18E-03 

7.20E-04 

5.58E-04 

Benzene 

7.53E-05 

3.33E-04 

4.51  E-04 

5.90E-04 

Toluene 

4.68E-05 

1.88E-04 

1.71E-04 

3.00E-04 

HCN 

ND 

ND 

2.29E-02 

6.34E-03 

Formaldehyde 

1.1  IE-04 

4.56E-04 

6.74E-03 

2.40E-03 

CEM 

C02 

1.47E+00 

5.10E+00 

2.40E-01 

4.35E-01 

CO 

2.50E-02 

4.80E-02 

1.71E-01 

1.84E-01 

NOx 

1.69E-02 

6.96E-03 

4.10E-02 

4.54E-02 
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Table  15:  Ratio  of  SCM  Predicted  to  Measured  Mixing  Ratios  for  Selected  Detonation  Products 


Species 

M031 

M456 

ML09 

ML15 

CO 

2.2 

2.3 

1.9 

2.3 

co2 

1.2 

1.4 

0.8 

1.3 

ch2o 

0.5 

0.4 

0.6 

0.4 

ch4 

1.9 

2.2 

1.8 

2.1 

HCN 

— 

— 

3.1 

2.8 

NOx 

0.5 

0.8 

0.7 

0.7 

C2H4 

2.7 

3.2 

2.8 

3.1 

c2h2 

0.2 

0.5 

0.3 

0.3 

CH3HCO 

0.6 

0.4 

0.5 

0.7 

Propene 

2.6 

2.3 

2.1 

2.1 

Ethanol 

0.3 

0.2 

0.2 

0.2 

Butanol 

0.3 

0.2 

0.3 

0.4 

Napthalene 

0.02 

0.03 

0.04 

0.03 

Pentane 

0.2 

0.1 

0.2 

0.3 

Hexane 

0.2 

0.2 

0.009 

0.1 

1  butene 

0.04 

0.04 

0.03 

0.03 

Isobutene 

0.03 

0.04 

0.04 

0.03 

Benzene 

0.04 

0.03 

0.02 

0.02 

Toluene 

0.01 

0.01 

0.009 

0.01 

m/p/o-xylene 

0.01 

0.01 

0.009 

0.02 

5.4  Field  Program  to  Identify  TRI  Chemicals 

Under  the  SERDP  project  CP- 1 197  entitled  "Field  Program  to  Identify  TRI  Chemicals  and  Determine 
Emission  Factors  from  DoD  munitions  activities  [Spicer  -  Battelle]  an  instrument  package  was  developed 
to  measure  a  target  list  of  Toxic  Release  Inventory  (TRI)  chemicals  required  as  part  of  federal  facilities 
reporting  requirements  under  Section  313  of  the  Environmental  Protection  Agency’s  Emergency  Planning 
and  Community  Right-to-Know  Act.  Following  testing  and  optimization  in  the  laboratory,  the  instrument 
package  was  to  be  deployed  in  a  mobile  laboratory  for  field  campaigns.  Two  field  campaigns  were 
planned.  The  first  measured  point  of  discharge  emissions.  The  second  field  campaign  measured  emissions 
point  of  impact  explosions.  Both  field  campaigns  were  conducted  at  the  Aberdeen  Test  Center  at  Aberdeen 
Proving  Ground,  MD  and  managed  by  USAEC. 

The  point  of  discharge  emissions  is  not  the  focus  of  the  SCM  and,  consequently,  will  not  be  discussed  here. 
However,  under  an  additional  task,  ARI  deployed  an  innovative  aerosol  mass  spectrometer  (AMS) 
instrument  to  support  those  field  studies  by  measuring  the  non-refractory  chemical  composition  of  fine 
particles  from  firing  point  muzzle  exhaust  plumes.  Muzzle  exhaust  plumes  contained  in  the  Aberdeen  Test 
Center’s  (ATC)  point-of-discharge  test  chamber  were  sampled  for  about  20  minutes  after  the  firing  of  the 
each  munition.  The  non-refractory  chemical  composition  of  the  evolving  muzzle  exhaust  plume  fine 
particles  size  distribution  was  characterized  during  these  measurements.  A  summary  of  these  results  and  a 
citation  for  the  task  report  we  presented  is  given  in  Appendix  A. 

The  second  field  campaign  for  point  of  impact  explosions,  however,  provided  additional  experimental  data 
which  could  be  used  in  test  and  validating  the  SCM.  Unfortunately,  unforeseen  delays  in  being  able  to 
conduct  the  field  tests  meant  that  data  wasn’t  available  until  the  last  few  months  of  effort  for  the  SCM 
development  and  validation  work.  Consequently,  the  results  presented  are  based  on  a  preliminary  analysis 
by  Battelle  of  the  field  data  and  may  change  with  further  analysis. 
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Table  16  lists  three  point  of  impact  munitions  for  which  field  data  from  CP- 1 197  was  available  for 
comparison  with  SCM  calculations.  The  designation,  detonation  medium,  and  number  of  tests  are  also 
shown. 


Table  16:  Munitions  tested  under  CP-1197  [C.  Spicer,  private  communication,  2003] 


Item 

Designation 

Medium 

No.  per  test 

No.  of  tests 

155  mm  howitzer 
rounds 

M107HE 

over  steel  plate 

1 

2 

4.2  in.  mortar 

M329A2HE 

over  soil 

1 

2 

4.2  in  mortar 

M329A2HE 

over  steel  plate 

1 

2 

Table  16  summarizes  the  comparison  between  SCM  predictions  and  measured  emissions  for  the  test  items 
in  Table  16  listing  the  ratio  of  calculated-to-measured  emission  factors  for  those  species  for  which  field 
data  was  available.  In  contrast  to  the  USAEC  exploding  ordinance  and  DPG  field  test  data,  SCM  over 
predicts  C02  relative  to  measured  emission  factors  by  a  larger  amount  considering  all  three  munitions 
together.  It  is  not  clear  as  to  whether  or  not  this  is  due  to  difficulties  in  integrating  over  the  plume  in  the 
field  tests,  differences  between  reported  and  actual  explosive  formulations,  or  an  over  prediction  of  C02  in 
detonation  plume  by  the  SCM  thermochemical  model. 


Table  17:  SCM  Predicted  to  Measured  Emission  Factors  for  CP-1197  Field  Test  Munitions. 


155-mm  howitzer 

5.2  in  mortar 

5.2  in  mortar 

over  steel  plate 

Over  steel  plate 

over  soil 

C02 

1.3 

1.3 

1.4 

CO 

2.2 

2.3 

2.6 

NOx 

0.6 

0.5 

0.5 

Propene 

0.7 

0.5 

0.6 

Formaldehyde 

0.8 

0.6 

0.7 

Benzene 

0.05 

0.3 

0.05 

Toluene 

0.01 

only  upper  limit 
reported 

0.01 

6.0  Summary 

This  report  has  documented  the  Source  Characterization  (SCM)  model  for  predicting  the  release,  chemical 
transformations,  deposition,  and  dispersion  of  emissions  from  FIE  detonation  of  munitions  and  explosive 
charges.  The  SCM  incorporates  a  number  of  detailed  microphysical  processes  that  are  important  in 
understanding  munition  emission  factors  and  the  potential  environmental  impact.  Of  particular  interest 
relative  to  previous  models  are  the  sub-models  for  crater  dynamics,  afterburning  chemistry,  and  the  impact 
of  detonation  plume  vortex  structure  of  particle  deposition.  Data  analysis  and  model  validation  studies  to- 
date  have  indicated  that  the  SCM  model  holds  much  promise  as  a  numerical  tool  for  understanding 
detonation  emissions  and  supporting  environmental  impact  assessments  at  DoD  test  and  training  sites. 
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9.0  Appendix  A  -  Aerosol  Mass  Spectrometric  Point-of-Fire  Measurements 


Aerosol  Mass  Spectrometer  Characterization  of  Point-of-Discharge 
Fine  Particle  Emissions 


Added  Task  Summary  for  Project  CP-1159 


T.  B.  Onasch,  H.  Boudries,  D.  R.  Worsnop  and  C.  E.  Kolb 
Center  for  Aerosol  and  Cloud  Chemistry 
Aerodyne  Research,  Inc 
45  Manning  Road 
Billerica,  MAO  1821 

August  2003 


An  Aerodyne  Aerosol  Mass  Spectrometer  (AMS)  was  used  to  measure  the  fine  particulate  mass  (PM1.0) 
emitted  from  point-of-discharge  emission  plumes  in  the  Aberdeen  Test  Center  (ATC)  enclosed  gun 
discharge  test  facility  from  September  25  to  October  2,  2002.  The  AMS  provided  real-time  measurements 
of  the  size-resolved  mass  distributions  for  fine  particles  between  40  and  1000  nm  in  diameter  and  the  total 
mass  loadings  of  organic  and  non-refractory  inorganic  species  (PM1.0).  These  tests  were  carried  out  as 
part  of  the  SERDP  project  CP- 1 197,  led  by  Battelle  Columbus  (C.  Spicer,  PI),  in  coordination  with  the 
Army  Environmental  Center  and  the  Aberdeen  Test  Center. 

Fine  particulate  emissions  from  the  muzzle  of  weapons  firing  five  different  munitions  were  measured. 
After  the  firing  of  each  shell,  the  emissions  were  isolated  in  a  steel  test  chamber  for  approximately  20 
minutes  for  sampling.  For  each  munitions  type  beyond  the  first,  an  initial  background  was  established 
within  the  chamber  and  then  two  shots  were  fired,  one  in  the  morning  and  one  in  the  afternoon.  Three 
replicates  of  the  first  munitions  type  were  fired.  The  initial  particulate  mass  loadings  for  each  shot  were 
very  high,  ranging  from  1000  to  18000  pg  m"3.  Figure  1  shows  the  total  fine  particulate  mass  loadings  for 
each  test.  According  to  the  Munitions  Items  Disposition  Action  System  (MIDAS)  database  the  chemical 
compositions  of  the  munitions  vary;  however,  the  majority  of  the  propellant  mass  is  nitrocellulose  and 
other  nitrogen-rich  organic  compounds.  The  munitions’  casings  and  projectiles  are  composed  mainly  of 
aluminum,  zinc,  steel,  nickel,  and  copper  alloys.  As  expected,  the  discharge  particulate  emissions 
measured  during  this  study  were  dominated  by  the  propellant  materials,  although  for  several  munitions 
there  was  apparently  measurable  particulate  mass  from  ablated  metal  alloy. 
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Figure  1.  The  total  fine  particulate  mass  loadings  as  measured  by  the  AMS  for  each  individual  point-of- 
discharge  test. 

The  main  volatile  and  semi-volatile  particulate  components  quantified  by  the  AMS  included  inorganic 
sulfates,  nitrates,  carbonates,  chlorides  and  the  total  organic  particulate  mass.  Specific  air  toxic  compounds 
identified  and  quantified  by  the  AMS  include  particle-bound  cyanides,  polycyclic  aromatic  hydrocarbons, 
and  several  heavy  metals.  The  AMS  oven  temperature  was  varied  to  help  identify  and  distinguish  between 
non-refractory  and  refractory  compounds.  The  AMS  flash  vaporizes  particulate  matter  from  a  heated  oven 
during  detection.  Nominal  oven  temperature  for  non-refractory  inorganic  and  organic  compounds  is  600°C. 
At  oven  temperatures  up  to  1000°C,  more  refractory  compounds,  such  as  the  sulfates  and  chlorides  of  zinc, 
lead,  and  antimony,  vaporize  more  readily  and  are  thus  easier  to  detect  and  identify. 

The  fine  particle  emissions  were  consistently  internally  mixed  particles  that  coagulated  and  dry  deposited 
with  time.  Particulate  mass  and  chemical  composition  results  show  a  strong  dependence  on  the  extent  of 
propellant  combustion.  Ample  evidence  exists  indicating  that  several  different  exhaust  plumes  underwent 
further  oxidation,  during  an  afterburning  event.  This  conclusion  is  supported  by  gas  phase  CO  and  CO2 
data.  Fine  particulate  emission  factors  (pg  PM/g  propellant)  are  reported  in  Table  1.  A  full  report 
presenting  the  data  and  detailed  analysis  has  been  prepared  and  distributed  [Onasch  et  al.,  2003]. 
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2471 

115 

Chloride 

4 

4 

298 

11 

7 

Organics 

225 

145 

2160 

811 

978 

Carbonates 

293 

87 

905 

322 

629 

Cyanide 

45 

8 

144 

80 

126 

PAH 

0 

0 

3 

1 

1 

Zinc 

90 

1234 

1131 

0 

1 

Lead 

15 

60 

12 

6 

3 

Antimony 

5 

23 

0 

1 

0 

Total  PM 

976 

2032 

5554 

4202 

2899 

Table  1  -  Fine  particulate  emission  factors  for  each  munitions  type  (pg  PM/g  propellant). 

Onasch,  T.B.,  H.  Boudries,  D.R.  Worsnop,  and  C.E.  Kolb,  Aerosol  Mass  Spectrometer  Characterization  of 
Point-of-Discharge  Fine  Particulate  Emissions,  AR1-RR-1382,  Aerodyne  Research,  Inc.,  Billerica,  2003. 
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10.0  Appendix  B:  Chemical  Kinetic  Mechanism 


ELEMENTS 

ATOMIC 

CONSIDERED 

WEIGHT 

1.  H 

1.00797 

2.  C 

12.0112 

3.  0 

15.9994 

4  .  N 

14.0067 

C 

P 

H 

H 

A 

A 

R 

SPECIES 

S 

G 

MOLECULAR 

CONSIDERED 

E 

E 

WEIGHT 

1  . 

C(S) 

S 

0 

12.01115 

2  . 

H2 

G 

0 

2.01594 

3. 

02 

G 

0 

31 . 99880 

4  . 

H20 

G 

0 

18.01534 

5. 

O 

G 

0 

15.99940 

6. 

H 

G 

0 

1.00797 

7  . 

OH 

G 

0 

17.00737 

8  . 

H02 

G 

0 

33.00677 

9. 

H202 

G 

0 

34.01474 

10. 

CH20 

G 

0 

30.02649 

11 . 

HCO 

G 

0 

29.01852 

12  . 

CO 

G 

0 

28.01055 

13. 

C02 

G 

0 

44.00995 

14  . 

N 

G 

0 

14.00670 

15. 

N2 

G 

0 

28.01340 

16. 

NO 

G 

0 

30.00610 

17  . 

N02 

G 

0 

46.00550 

18  . 

HNO 

G 

0 

31.01407 

19. 

HONO 

G 

0 

47.01347 

20. 

N20 

G 

0 

44.01280 

21 . 

N03 

G 

0 

62.00490 

22  . 

HN03 

G 

0 

63.01287 

23. 

CH4 

G 

0 

16.04303 

24  . 

CH3 

G 

0 

15.03506 

25. 

CH30 

G 

0 

31.03446 

26. 

CH2 

G 

0 

14.02709 

27  . 

CH 

G 

0 

13.01912 

28  . 

CH20H 

G 

0 

31.03446 

29. 

CH30H 

G 

0 

32.04243 

30. 

C2H6 

G 

0 

30.07012 

31 . 

C2H5 

G 

0 

29.06215 

32  . 

C2H4 

G 

0 

28.05418 

33. 

C2H3 

G 

0 

27.04621 

TEMPERATURE  ELEMENT  COUNT 


LOW 

HIGH 

H 

C 

O 

N 

300.0 

5000.0 

0 

1 

0 

0 

300.0 

5000.0 

2 

0 

0 

0 

300.0 

5000.0 

0 

0 

2 

0 

300.0 

5000.0 

2 

0 

1 

0 

300.0 

5000.0 

0 

0 

1 

0 

300.0 

5000.0 

1 

0 

0 

0 

300.0 

5000.0 

1 

0 

1 

0 

200.0 

3500.0 

1 

0 

2 

0 

300.0 

5000.0 

2 

0 

2 

0 

300.0 

5000.0 

2 

1 

1 

0 

300.0 

5000.0 

1 

1 

1 

0 

300.0 

5000.0 

0 

1 

1 

0 

300.0 

5000.0 

0 

1 

2 

0 

200.0 

6000.0 

0 

0 

0 

1 

300.0 

5000.0 

0 

0 

0 

2 

200.0 

6000.0 

0 

0 

1 

1 

200.0 

6000.0 

0 

0 

2 

1 

200.0 

6000.0 

1 

0 

1 

1 

300.0 

5000.0 

1 

0 

2 

1 

200.0 

6000.0 

0 

0 

1 

2 

300.0 

5000.0 

0 

0 

3 

1 

300.0 

5000.0 

1 

0 

3 

1 

300.0 

5000.0 

4 

1 

0 

0 

300.0 

5000.0 

3 

1 

0 

0 

300.0 

3000.0 

3 

1 

1 

0 

300.0 

5000.0 

2 

1 

0 

0 

300.0 

5000.0 

1 

1 

0 

0 

250.0 

4000.0 

3 

1 

1 

0 

300.0 

5000.0 

4 

1 

1 

0 

300.0 

1500.0 

6 

2 

0 

0 

300.0 

5000.0 

5 

2 

0 

0 

300.0 

5000.0 

4 

2 

0 

0 

300.0 

5000.0 

3 

2 

0 

0 
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34  . 

C2H2 

G 

0 

26.03824 

300.0 

5000.0 

2 

2 

0 

0 

35. 

C2H 

G 

0 

25.03027 

300.0 

5000.0 

1 

2 

0 

0 

36. 

CH2CO 

G 

0 

42.03764 

200.0 

3500.0 

2 

2 

1 

0 

37  . 

CH3CO 

G 

0 

43.04561 

300.0 

5000.0 

3 

2 

1 

0 

38  . 

CH3HCO 

G 

0 

44.05358 

300.0 

5000.0 

4 

2 

1 

0 

39. 

CN 

G 

0 

26.01785 

300.0 

5000.0 

0 

1 

0 

1 

40. 

HCN 

G 

0 

27.02582 

300.0 

5000.0 

1 

1 

0 

1 

41 . 

HNCO 

G 

0 

43.02522 

300.0 

5000.0 

1 

1 

1 

1 

42  . 

NH 

G 

0 

15.01467 

300.0 

5000.0 

1 

0 

0 

1 

43. 

NH2 

G 

0 

16.02264 

300.0 

5000.0 

2 

0 

0 

1 

44  . 

NH3 

G 

0 

17.03061 

300.0 

5000.0 

3 

0 

0 

1 

(k  =  A 

T**b  exp (-E/RT) ) 

REACTIONS  CONSIDERED 

A 

b  E 

1 . 

H2+M=H+H+M 

4.57E+19 

-1.4  104000.0 

H2 

Enhanced  by 

2 . 500E+00 

H20 

Enhanced  by 

1 .200E+01 

CO 

Enhanced  by 

1 . 900E+00 

C02 

Enhanced  by 

3 . 800E+00 

2  . 

0+H20=0H+0H 

2 . 97E+06 

2.0  13400.0 

3. 

0+H2=H+0H 

5 . 06E+04 

2.7  6290.0 

4  . 

0+0+M=02+M 

6. 17E+15 

o 

LO 

1 

H2 

Enhanced  by 

2 . 500E+00 

H20 

Enhanced  by 

1 .200E+01 

CO 

Enhanced  by 

1 . 9O0EtOO 

C02 

Enhanced  by 

3 . 800E+00 

5. 

H+02=0+0H 

1 . 94E+14 

.0  16440.0 

6. 

H+02 ( +M) =H02 ( +M) 

4.52E+13 

o 

o 

Low  pressure 

limit:  .67000Et20  -.14200E+01 

.  00000E+00 

TROE  centering:  .10000E+01  . 10000E-89 

.  10000E+91 

H2 

Enhanced  by 

2 . 500E+00 

H20 

Enhanced  by 

1 .200E+01 

CO 

Enhanced  by 

1 . 900E+00 

C02 

Enhanced  by 

3 . 800E+00 

7  . 

H+0+M=OH+M 

4.72E+18 

O 

o 

\ — 1 

1 

H2 

Enhanced  by 

2 . 500E+00 

H20 

Enhanced  by 

1 .200E+01 

CO 

Enhanced  by 

1 . 900E+00 

C02 

Enhanced  by 

3 . 800E+00 

8  . 

0H+H2=H20+H 

2 . 16E+08 

1.5  3430.0 

9. 

OH+H+M=H20+M 

2 . 2 1E+22 

-2.0  .0 

H2 

Enhanced  by 

2 . 500E+00 

H20 

Enhanced  by 

1 .200E+01 

CO 

Enhanced  by 

1 . 900E+00 

C02 

Enhanced  by 

3 . 800EtOO 

10. 

H02+0=02+0H 

1.75E+13 

.0  -397.0 

11 . 

H02+H=H2+02 

6. 62E+13 

.0  2130.0 

12  . 

H02+H=OH+OH 

1 . 69E+14 

.0  874.0 

13. 

H02+0H=H20+02 

1 . 90E+16 

o 

o 

\ — 1 

1 

14  . 

H02+H02=H202+02 

4 . 20E+14 

.0  11980.0 

Declared  duplicate  reaction. 

.  . 

15. 

H02+H02=H202+02 

1.30E+11 

.0  -1629.0 

Declared  duplicate  reaction. 

.  . 
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16. 

H202 ( +M) =OH+OH (+M) 

2 . 95E+14 

.0 

48460.0 

Low  pressure  limit: 

:  . 12000Etl8 

. OOOOOEtOO 

. 45500E+05 

TROE  centering: 

. 50000E+00 

. 10000E-89 

. 10000Et91 

17  . 

H202+0=0H+H02 

9. 64E+06 

2.0 

3970.0 

18  . 

H202+H=H20+0H 

1.00E+13 

.0 

3590.0 

19. 

H202+H=H02+H2 

4 . 82E+13 

.0 

7950.0 

20. 

H202+0H=H20+H02 
Declared  duplicate 

reaction . . . 

1.00E+12 

.0 

.0 

21 . 

H202+0H=H20+H02 
Declared  duplicate 

reaction . . . 

5.80E+14 

.0 

9557.0 

22  . 

CH20+M=HC0+H+M 

1 . 63E+33 

-4 . 1 

92550.0 

23. 

CH20tM=H2+C0+M 

8 . 25E+15 

.0 

69540.0 

24  . 

CH20+02=HC0+H02 

2.05E+13 

.0 

38920.0 

25. 

CH20+0=HC0+0H 

1 . 81E+13 

.0 

3078.0 

26. 

CH20tH=HC0+H2 

1 .26E+08 

1 . 6 

2163.0 

27  . 

CH20+0H=HC0+H20 

3.43E+09 

1.2 

-447.0 

28  . 

CH20+H02=HC0+H202 

1 . 99E+12 

.0 

11660.0 

29. 

HCO+M=H+COtM 

1 . 85E+17 

O 

\ — 1 

1 

17000.0 

H2 

Enhanced  by 

1 . 890E+00 

H20 

Enhanced  by 

1 .200E+01 

CO 

Enhanced  by 

1 . 900E+00 

C02 

Enhanced  by 

3 . 800E+00 

30. 

HC0+02=C0+H02 

7.58E+12 

.0 

406.0 

31 . 

HC0+0=C0+0H 

3.00E+13 

.0 

.0 

32  . 

HC0+0=C02+H 

3.00E+13 

.0 

.0 

33. 

HCO+H=CO+H2 

7 . 23E+13 

.0 

.0 

34  . 

HC0+0H=C0+H20 

3.00E+13 

.0 

.0 

35. 

HC0+H02=C02+0H+H 

3.00E+13 

.0 

.0 

36. 

CO+O ( +M) =C02 ( +M) 

1 . 80E+10 

.0 

2380.0 

Low  pressure  limit:  .13500E+25  -.27880E+01  .41910E+04 
TROE  centering:  .10000E+01  .10000E-89  .10000E+91 


N2 

Enhanced  by 

1 . 330E+00 

H2 

Enhanced  by 

2 . 500E+00 

H20 

Enhanced  by 

1 .200E+01 

CO 

Enhanced  by 

1 . 900E+00 

C02 

Enhanced  by 

3 . 800E+00 

37  . 

C0+02=C02+0 

2.53E+12 

.0 

47700.0 

38  . 

CO+OH=C02+H 

1.50E+07 

1.3 

-765.0 

39. 

C0+H02=C02+0H 

5.80E+13 

.0 

22930.0 

40. 

N+02=N0+0 

6 . 40Et09 

1.0 

6280.0 

41 . 

N+OH=NO+H 

3.80E+13 

.0 

.0 

42  . 

N+H02=NO+OH 

1.00E+13 

.0 

2000.0 

43. 

N+C02=N0+C0 

1 . 90E+11 

.0 

3400.0 

44  . 

N+NO=N2+0 

3.27E+12 

.3 

.0 

45. 

N+N02=NO+NO 

4.00E+12 

.0 

.0 

46. 

N+N02=N20+0 

5.00E+12 

.0 

.0 

47  . 

N+N02=N2+02 

1.00E+12 

.0 

.0 

48  . 

N+HNO=N20+H 

5 . OOEtlO 

.5 

3000.0 

49. 

N+N20=N2+N0 

1.00E+13 

.0 

19870.0 

50. 

NO+M=N+0+M 

9. 64Etl4 

.0 

148400.0 

N2  Enhanced  by  1.500E+00 

C02  Enhanced  by  2.500E+00 

51.  NO+O (+M) =N02 (+M)  1.30E+15  -.8  .0 

Low  pressure  limit:  .47200E+25  -.28700E+01  .15510E+04 

TROE  centering:  .95700E+00  .10000E-89  .83320E+04 

52.  NO+H (+M) =HNO (+M)  1.52E+15  -.4  .0 

Low  pressure  limit:  .89600E+20  -.13200E+01  .73520E+03 
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TROE  centering: 

.  82000Et00 

.  10000E-89 

.  10000E+91 

53. 

NO+OH ( +M) =HONO (+M) 

1 . 99E+12 

-.1  -721.0 

Low  pressure  limit: 

.  50800E+24 

-  .  25100E+01  - 

.  67560Et02 

TROE  centering: 

.  62000Et00 

.  10000E-89 

.  10000E+91 

H20  Enhanced 

by  5.000E+00 

54  . 

N0+H02=N02+0H 

2 . 11E+12 

.0  -479.0 

55. 

NO+HCO=HNO+CO 

7 . 23E+12 

.0  .0 

56. 

N02+0=02+N0 

3.91E+12 

.0  -238.0 

57  . 

N02+0 ( +M) =N03 ( +M) 

1.33E+13 

.0  .0 

Low  pressure  limit: 

.  14900E+29 

-  .  40800E+01 

.  24670E+04 

TROE  centering: 

.  82 600E+00 

.  10000E-89 

.  31910E+04 

58  . 

N02+H=N0+0H 

1.32E+14 

.0  361.6 

59. 

N02+0H (+M) =HN03 (+M) 

2.41E+13 

.0  .0 

Low  pressure  limit: 

.  64200Et33 

-  .  54900E+01 

.  23500E+04 

TROE  centering: 

.  83700E+00 

.  10000E-89 

.  16570E+04 

60. 

N02+CH20=H0N0+HC0 

8 . 02E+02 

2.8  13730.0 

61 . 

N02+HC0=C0tH0N0 

1 .24E+23 

-3.3  2354.0 

62  . 

N02+HC0=H+C02+N0 

8.39E+15 

-.8  1927.0 

63. 

N02+C0=C02+N0 

9.03E+13 

.0  33780.0 

64  . 

N02tN02=N03+N0 

9. 64Et09 

.7  20920.0 

65. 

N02+N02=2N0+02 

1 . 63E+12 

.0  26120.0 

66. 

HN0+02=N0+H02 

1.00E+13 

.0  25000.0 

67  . 

HN0+0=0H+N0 

1 . 81E+13 

.0  .0 

CO 

'vO 

HNO+H=H2+NO 

1 . 81E+13 

.0  993.5 

69. 

HN0+0H=H20tN0 

1.00E+13 

.0  993.5 

70. 

HN0+HC0=CH20+N0 

6.02E+11 

.0  1987.0 

71 . 

HN0+N0=N20+0H 

2 . OOEL12 

.0  26000.0 

72  . 

HN0+N02=H0N0+N0 

6.02E+11 

.0  1987.0 

73. 

HN0+HN0=H20+N20 

8.51E+08 

.0  3080.0 

74  . 

H0N0+0=0H+N02 

1.20E+13 

.0  5961.0 

75. 

H0N0+H=H2+N02 

1 .20E+13 

.0  7352.0 

76. 

H0N0+0H=H20+N02 

1 .26E+10 

1.0  135.1 

77  . 

N20 (+M) =N2+0 ( +M) 

7 . 91E+10 

.0  56020.0 

Low  pressure  limit: 

.  91300E+15 

. 00000E+00 

.  57 690E+05 

H20 

Enhanced  by 

7 . 500E+00 

NO 

Enhanced  by 

2 . 000E+00 

CO 

Enhanced  by 

2 . 000E+00 

C02 

Enhanced  by 

3 . 000E+00 

78  . 

N20+0=02+N2 

1.00E+14 

.0 

28000.0 

79. 

N20+0=2N0 

1.00E+14 

.0 

28000.0 

o 

CO 

N20+H=N2+0H 
Declared  duplicate 

reaction . . . 

2.53E+10 

.0 

4550.0 

81 . 

N20+H=N2+0H 
Declared  duplicate 

reaction . . . 

2 . 23E+14 

.0 

16750.0 

82  . 

N20+0H=H02+N2 

2.00E+12 

.0 

40000.0 

83. 

N20+C0=N2+C02 

5.01E+13 

.0 

44000.0 

84  . 

CH4+02=CH3tH02 

3.97E+13 

.0 

56687.4 

LO 

CO 

CH4+H=CH3+H2 

1.32E+04 

3.0 

8009.1 

'X> 

CO 

CH4+0=CH3+0H 

7 . 22E+08 

1 . 6 

8454 . 6 

87  . 

CH4+0H=CH3tH20 

1.57E+07 

1 . 8 

2772.0 

CO 

CO 

CH4+H02=CH3+H202 

9.03E+12 

.0 

24631.2 

O'! 

CO 

CH3+M=CH2+H+M 

1.02E+16 

.0 

90288.0 

90. 

CH3+02=CH30+0 

1.32E+14 

.0 

31284.0 

91 . 

CH4 ( +M) =CH3+H (+M) 

2.40E+16 

.0 

104544.0 

Low  pressure  limit: 

:  . 47000EL48 

-  .  82000E+01 

. 11722E+06 

TROE  centering: 

. 69000EL00 

. 90000E+02 

. 22100E+04 

92 . 

CH3+H=CH2+H2 

6.02E+13 

.0 

15048.0 
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93. 

CH3+0=CH20+H 

8.43E+13 

.0 

.0 

94  . 

CH3+0H=CH2+H20 

7 . 22E+12 

.0 

2772.0 

95. 

CH3+H02=CH30+0H 

1 . 81E+13 

.0 

.0 

96. 

CH30+H=CH20+H2 

1 . 81E+13 

.0 

.0 

97  . 

CH2+02=C02+H+H 

1 . 60E+12 

.0 

1000.0 

98  . 

CH2+02=CH20+0 

5.00E+13 

.0 

9000.0 

99. 

CH2+02=C02+H2 

6.90E+11 

.0 

500.0 

100. 

CH2+02=C0+H20 

1 . 90E+10 

.0 

-1000.0 

101 . 

CH2+02=C0+0H+H 

8 . 60E+10 

.0 

-500.0 

102  . 

CH2+02=HC0+0H 

4.30E+10 

.0 

-500.0 

103. 

CH2+H=CH+H2 

3.01E+13 

.0 

.0 

104  . 

CH2+0=C0+H+H 

1.51E+13 

.0 

.0 

105. 

CH2+0=C0+H2 

1.51E+13 

.0 

.0 

106. 

CH2+0H=CH20+H 

3.01E+13 

.0 

.0 

107  . 

CH2+HCO=CH3+CO 

1 . 81E+13 

.0 

.0 

108  . 

CH+0=C0+H 

3.97E+13 

.0 

.0 

109. 

CH3+CH3 ( +M) =C2H6 (+M) 

3.61E+13 

.0 

.0 

Low  pressure  limit: 

. 12700E+42 

-  .  70000E+01 

. 27522E+04 

TROE  centering: 

. 62000E+00 

. 73000E+02 

. 11800E+04 

110. 

C2H6+02=C2H5+H02 

6.02E+13 

.0 

51678.0 

Ill . 

C2H6+0=C2H5+0H 

9 . 99E+08 

1.5 

5781 . 6 

112  . 

C2H6+0H=C2H5+H20 

7 . 22E+06 

2.0 

861.3 

113. 

C2H6tH02=C2H5tH202 

1.32E+13 

.0 

20395.0 

114  . 

C2H6+CH3=C2H5+CH4 

1.51E-07 

6.0 

6025.1 

115. 

C2H5tH2=C2H6+H 

3.07E+00 

3.6 

8415.0 

116. 

C2H5+02=C2H4+H02 

1.02E+10 

.0 

-2178.0 

117  . 

C2H5+0=CH20+CH3 

1 . 13E+13 

.0 

.0 

118  . 

C2H5+H=CH3+CH3 

3.61E+13 

.0 

.0 

119. 

C2H5+CH3=CH4+C2H4 

1 . 14E+12 

.0 

.0 

120. 

C2H5+C2H5=C2H4+C2H6 

1.44E+12 

.0 

.0 

121 . 

C2H4+M=C2H3+H+M 

2.59E+17 

.0 

96228.0 

122  . 

C2H4+M=C2H2+H2tM 

3.49E+16 

.0 

71280.0 

123. 

C2H4+0=CH3+HC0 

4.74E+06 

1 . 9 

178.2 

124  . 

C2H4+0=CH2+CH20 

8 . 13E+06 

1 . 9 

178.2 

125. 

C2H4+0=CH2C0+H2 

6.77E+05 

1 . 9 

178.2 

126. 

C2H4+H=C2H3+H2 

5.42E+14 

.0 

14850.0 

127  . 

C2H4+0H=C2H3+H20 

2.05E+13 

.0 

5920.2 

128  . 

C2H4+CH3=C2H3+CH4 

4 . 15E+12 

.0 

11088.0 

129. 

C2H3+0=C2H2+0H 

1.00E+13 

.0 

.0 

130. 

C2H3+0=CH3+C0 

1.00E+13 

.0 

.0 

131 . 

C2H3+0=CH2+HC0 

1.00E+13 

.0 

.0 

132  . 

C2H3+H=C2H2+H2 

1 .20E+13 

.0 

.0 

133. 

C2H2+0H=C2H+H20 

3.01E+13 

.0 

12870.0 

134  . 

C2H2+OH=CH2CO+H 

3.01E+13 

.0 

12870.0 

135. 

C2H+0=CH+C0 

1.02E+13 

.0 

.0 

136. 

C2H+H2=C2H2+H 

1.08E+13 

.0 

2158.2 

137  . 

CH2CO+H=CH3+CO 

1 . 81E+13 

.0 

3366.0 

138  . 

CH2C0+0H=CH20+HC0 

1.02E+13 

.0 

.0 

139. 

NH2+02=HN0+0H 

1.78E+12 

.0 

14900.0 

140. 

NH2+0=HN0+H 

6. 63E+14 

-.5 

.0 

141 . 

NH2+0=NH+0H 

6.75E+12 

.0 

.0 

142  . 

NH2+H=NH+H2 

6. 92E+13 

.0 

3650.0 

143. 

NH2+0H=NHtH20 

4 . 00E+06 

2.0 

1000.0 

144  . 

NH2+N=N2+2H 

7 . 20E+13 

.0 

.0 

145. 

NH2+N0=N2+H20 

Declared  duplicate  reaction... 

1.30E+16 

-1.3 

.0 

146. 

NH2+N0=N2+H20 

•k-k'k-k'k-k'k-k 

-.6 

.0 
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147  . 

NH2+N0=N20tH2 

5.00E+13 

.0 

24640.0 

148  . 

NH2+NO=HNO+NH 

1.00E+13 

.0 

40000.0 

149. 

NH2+N02=N20+H20 

3.28E+18 

-2.2 

.0 

150. 

NH3+M=NH2+H+M 

2 .20E+16 

.0 

93470.0 

151 . 

NH3+0=NH2+0H 

9.40E+06 

1 . 9 

6460.0 

152  . 

NH3+H=NH2+H2 

6 . 40E+05 

2.4 

10170.0 

153. 

NH3+0H=NH2+H20 

2.04E+06 

2.0 

566.0 

154  . 

NH3+H02=NH2+H202 

3.00E+11 

.0 

22000.0 

155. 

NH2+H02=NH3+02 

1.00E+13 

.0 

.0 

156. 

NH2+NH2=NH3+NH 

5.00E+13 

.0 

10000.0 

157  . 

HCN ( +M) =H+CN ( +M) 

8.30E+17 

-.9 

123800.0 

Low  pressure  limit: 

.  35700E+27 

-  .  26000E+01 

. 12490E+06 

TROE  centering: 

. 95700E+00 

.  10000E-89 

. 83320E+04 

158  . 

HCN+0=CN+0H 

2.70E+09 

1 . 6 

29200.0 

159. 

HCN+H=H2+CN 

3.73E+14 

.0 

24850.0 

160. 

HCN+0H=H20+CN 

3.90E+06 

1 . 8 

10290.0 

161 . 

HCN+OH=NH2+CO 

7 . 83E-04 

4.0 

4000.0 

162  . 

HCN+HC0=CH20+CN 

2.41E+12 

.0 

34180.0 

NOTE:  A  units  mole-cm-sec-K,  E  units  cal/mole 
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